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PREFACE 


This  publication  is  regarded  by  the  author  as  supple¬ 
mentary  training  material  for  students  who  are  familiar 
with,  or  are  studying,  elementary  theory  of  sampling  includ¬ 
ing  stratification,  cluster  sampling,  ratio  and  regression 
estimation,  sampling  with  probability  proportional  to  size, 
and  multiple-stage  sampling.  After  studying  sampling  methods 
one  at  a  time,  it  is  important  to  get  a  unified  view  of  the 
several  methods  and  the  conditions  under  which  they  have 
about  the  same  or  different  variances. 

In  sampling  various  populations  we  quite  often  find  two 
or  more  techniques  that  are  roughly  equal  in  efficiency  and 
reduce  sampling  variance  about  as  much  as  possible.  Admin¬ 
istrative  feasibility,  costs,  and  freedom  from  potential 
bias  are  important  criteria  for  selecting  a  sampling  plan 
and  become  primary  criteria  when  the  choice  fs  among  plans 
having  small  differences  in  sampling  variance. 

Ability  to  prejudge  accurately  the  efficiency  of  alter¬ 
native  sample  designs  with  reference  to  various  survey 
objectives  and  populations  is  important.  Such  ability  comes 
from  experience  and  detailed  study  of  alternative  techniques 
of  sampling  a  population  and  of  making  estimates.  Quite 
often  only  two  or  three  alternatives  are  compared  in  an 
analysis  because  of  limitations  of  data  or  only  a  few  alter¬ 
natives  are  of  interest.  In  this  publication  many  alternative 
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sampling  and  estimation  plans  are  applied  to  a  small  popula¬ 
tion  of  apple  trees  and  the  results  are  recorded  in  tables 
for  comparative  purposes.  The  focus  of  attention  is  on  the 
magnitude  of  the  differences  in  efficiency  in  relation  to 
the  patterns  of  variation  that  exist. 

For  some  readers,  parts  of  the  presentation  are  probably 
too  detailed.  However,  it  is  important  to  understand  fully 
the  alternatives  and  to  put  mathematical  expressions  for 
estimators  and  their  variances  in  forms  that  are  most  meaning¬ 
ful  for  comparative  purposes.  Exercises  are  distributed 
through  the  text. 

Chapter  I  makes  use  of  graphical,  or  geometrical,  inter¬ 
pretations  in  the  comparison  of  four  alternative  ways  of  using 
an  auxiliary  variable.  There  is  a  brief  presentation  of  the 
relevant  theory  for  each  plan  which  is  followed  by  a  dis¬ 
cussion  of  the  plans  including  a  numerical  example.  Sampling 
with  probability  proportional  to  size  in  comparison  to  other 
methods  is  of  special  interest.  For  comparison,  a  part  of 
each  variance  formula  is  written  as  the  sum  of  squares  of 
deviations  from  a  line. 

Chapter  II  expands  the  comparisons  made  in  Chapter  I  to 
include  interactions  in  efficiency.  For  example,  the  compar¬ 
ative  efficiency  of  sampling  units  of  various  size  is  related 
to  the  method  of  estimation  and  to  stratification.  Chapter  III 
provides  some  further  comparisons,  but  the  emphasis  is  on  how 


11 


theory  and  ingenuity  solved  an  important  problem  in  the  samp¬ 
ling  of  fruit  trees.  Some  comparisons  involving  two-stage 
sampling  using  apple  trees  as  an  example  are  included  in 
Chapter  IV. 

This  volume  was  written  because  it  was  a  pleasure  and 
because  I  always  learn  something  from  making  comparisons 
like  those  contained  herein. 


Earl  E.  Houseman 
Statistician 
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CHAPTER  I 


SIMPLE  USES  OF  AN  AUXILIARY  VARIABLE 


1.1  INTRODUCTION 

Proficiency  in  the  use  of  auxiliary  information  to  reduce 
sampling  variance  is  an  important  goal  in  the  formulation  of 
a  sampling  plan.  In  this  chapter  we  will  compare  four  alterna¬ 
tive  methods  of  using  an  auxiliary  variable  in  the  design  of  a 
sample  or  in  the  estimator  and  one  without  using  an  auxiliary 
variable,  giving  a  total  of  five  alternative  methods.  The 
methods  discussed  are  commonly  found  in  textbooks  on  sampling. 
It  is  important  to  know  whether  an  auxiliary  variable  is  worth 
using  and  how  to  use  it  most  effectively.  Achievement  of 
greater  efficiency  in  the  use  of  an  auxiliary  variable  is 
usually  inexpensive  compared  to  increasing  sample  size,  but 
incorrect  use  could  cause  an  increase  rather  than  a  decrease 
in  sampling  error. 

For  each  of  the  five  alternatives  there  is  an  estimator 
and  the  variance  of  each  estimator  can  be  expressed  in  a  form 
that  is  suitable  for  interpretation  of  the  sampling  variance 
as  a  function  of  deviations  of  points  from  a  line.  The  emphasi 
in  this  chapter  is  on  simple  dot  charts  as  a  useful  aid  to 
understanding  or  judging  the  comparative  effectiveness  of  alter 
native  methods  in  different  situations.  Special  attention  will 
be  given  to  sampling  with  probability  proportional  to  size  and 
how  it  compares  with  other  ways  of  using  an  auxiliary  variable 


including  stratification  and  optimum  allocation.  After  a 
review  of  notation,  definitions,  and  theory,  a  numerical 
example  will  be  presented  which  makes  use  of  some  data  col¬ 
lected  in  a  research  project  to  develop  techniques  for 
estimating  apple  production. 

Consider  a  population  of  N  sampling  units  and  let  Y^,..., 
represent  the  unknown  values  of  Y  and  let  X^,...,  X^  repre¬ 
sent  the  known  values  of  an  auxiliary  variable  X.  A  sample 
is  to  be  selected  and  the  values  of  Y  for  the  n  su's  (sampling 
units)  in  the  sample,  namely  y^,...,  yn,  are  to  be  obtained. 
The  corresponding  values  of  X  for  the  su's  in  the  sample  are 

x,  , .  .  .  ,  x  .  We  assume  that  the  objective  is  to  estimate  the 
I  n  J 

N 

IY. 

population  mean,  Y  =  — ^  .  Also,  in  the  interest  of  keeping 

the  notation  as  simple  as  possible,  let  Y  and  X  represent  the 

N  N 

population  totals.  That  is,  Y  =  ZY^  and  X  =  £X^.  This  gives 
"Y",  for  example,  a  dual  meaning  as  in  "the  characteristic  Y" 
or  as  the  total  for  the  population.  However,  the  meaning 
should  be  clear  from  the  context. 

A  resume  of  the  theory  for  each  of  the  five  alternatives, 
which  will  be  called  plans,  is  presented  after  a  brief  review 
of  sampling  with  equal  and  unequal  probabilities  of  selection. 
1.1.1  EQUAL  PROBABILITIES  OF  SELECTION 

A  sample  obtained  by  selecting  one  su  at  a  time,  at 
random  with  equal  probability  and  without  replacement,  is 
called  a  simple  random  sample.  When  the  variance  of  Y  in  the 
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population  is  defined  as 


2 

N 

£  (Y. 

v  l 

-  2 
-  Y)z 

a  = 

N 

the 

variance  of 

the 

mean,  y, 

of  a 

simpl 
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given  by 

V  (y) 
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the  variance 

of 

Y  is  defined  as 

N 

s2  = 

l  (Y. 

v  l 

-  Y) 

N-l 

and 

the  variance 

of 

y  is 

V(y) 

N-n 

N 

£ 

n 

(1.1) 

random  sample  of  n 

(1.2) 


(1.3) 


(1.4) 


2 

In  the  discussion  that  follows,  S  will  be  used  as  the 
definition  of  the  variance  of  y. 

The  mean,  y,  of  a  simple  random  sample  is  an  unbiased 

n  -  2 
2  2^i-y) 

estimate  of  Y  and  the  variance,  s  =  — -  ,  among  su's  in 

2 

the  sample  is  an  unbiased  estimate  of  S  .  Incidentally,  the 
writer  from  a  practical  point  of  view  advises  use  of  the  word 
’’unbiased”  with  some  caution.  In  the  mathematical  theory,  the 
meaning  of  ’’unbiased”  is  usually  clear,  but  in  practice  ’’un¬ 
biased  estimate”  is  often  misleading  to  persons  who  are 
interested  in  estimates  from  a  survey  and  are  unaware  of  the 
restricted  meaning  of  the  term.— ^ 


1/  See  sections  4.4  and  4.5  of  Expected  Value  of  a  Sample  Esti¬ 
mate,  Statistical  Reporting  Service,  USDA,  September  1974. 
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Exercise  1.1  Show  that  either  definition  of  the  variance 
among  the  N  values  of  Y  leads  to  the  same  answer  for  the  vari¬ 
ance  of  y.  That  is ,  show  that  equations  1.2  and  1.4  are  the 
same . 

1.1.2  UNEQUAL  PROBABILITIES  OF  SELECTION 

Some  sampling  plans  specify  that  sampling  units  be 

selected  with  pps  (probability  proportional  to  size).  For 

simplicity,  sampling  with  replacement  is  assumed. 

It  is  often  very  important  to  make  a  clear  distinction 

t  h 

between  the  probability  of  selecting  the  i  su  of  a  popula¬ 
tion  when  a  particular  random  draw  is  made  and  the  probability 
th 

of  the  l  su  being  included  in  a  sample.  To  help  make  the 
distinction  clear,  the  letter  "P"  or  "p"  will  be  used  to  repre¬ 
sent  selection  probability  and  "f"  will  represent  inclusion 
probability,  that  is,  the  probability  of  any  given  su  being  in 
the  sample.  When  simple  random  sampling  is  applied,  each  su 
has  a  probability  equal  to  ^  of  being  in  the  sample.  That  is, 
the  inclusion  probability,  f,  is  equal  to  for  simple  random 
sampling. 

With  regard  to  sampling  with  pps  and  replacement,  let  P^, 

P2>...,P^  be  the  set  of  selection  probabilities  for  the  N  su's 

N 

in  the  population.  It  is  specified  that  EPi=l.  Thus,  "select¬ 
ing  a  sample  with  probabilities  proportional  to  X^"  means  that 
X^^  N 

P^  =  where  X=£X^.  Since  the  sampling  is  with  replacement, 
the  selection  probabilities  remain  constant  from  one  random 
draw  to  another. 
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The  unbiased  estimator  of  Y  for  a  sample  of  n  is 

1  1  n  yi 

?  =  W  (n’  z  57 


(1.5) 


In  the  estimator,  i  is  an  index  of  the  n  random  draws  because 
the  same  su  might  be  selected  more  than  once.  To  illustrate, 
if  on  the  4th  draw  su  number  15  in  the  population  is  selected, 
y ^  and  p4  are  equal  to  Y and  And  if  the  15th  su  is 

selected  again  on  the  12th  draw,  y ^  an<l  P^2  are  eclual  t0  Y^,. 
and  P^<-.  In  practice,  techniques  for  avoiding  the  selection 
of  the  su  more  than  once  are  usually  introduced  but  such  tech¬ 
niques  are  for  later  consideration. 

Yi 

Each  of  the  n  values  of  —  in  Eq.  1.5  is  an  unbiased  esti- 

P  • 

1  n  yi 

mate  of  the  population  total.  Thus ,  (— )  E  —  is  a  simple 

1  i 

average  of  n  independent,  unbiased  estimates  of  Y,  and  (^) 
appears  in  Eq.  1.5  so  y  will  be  an  estimator  of  Y  instead  of 
the  population  total. 


The  variance  of  y,  Eq.  1.5,  is 

V (y)  =  — 

77  n 


where 


and 


-tNY*  ^  -i  ry 
°  -  (iyjEP, (pi  -Y)Z  =  ( iy)aZ 

1  l  NZ  L 


(1.6) 

(1.7) 


N 

Y  =  ZY. 

l 


Is  Eq.  1.6  reasonable?  Study  the  estimator.  For  any 
Yi 

given  value  of  i,  —  in  repeated  sampling  is  a  random  variable 


which  has  an  expected  value  equal  to  the  population  total,  Y. 
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y. 


By  definition,  the  variance  of  — -  is 

^i 


N 


°t  ■  IpiCpr  ■  Y) 


where  i  is  the  index  to  the  N  su's  in  the  population.  In  the 


i  n  y  • 

estimator,  —  l  —  is  the  average  of  n  independent  estimates; 
n  Pi  2 

therefore,  the  variance  of  this  average  is  —  .  And,  since  we 

2 

are  interested  in  estimating  Y  rather  than  Y,  a  must  be  divided 


by  N  as  shown  in  Eq.  1.7. 


2/ 


1.2  RESUME  OF  THEORY  FOR  FIVE  PLANS 

2 

As  discussed  above,  we  will  use  S  ,  Eq.  1.3,  as  the  defini¬ 
tion  of  the  population  variance  for  simple  random  sampling  with 

2 

replacement  and  a  ,  Eq.  1.7,  is  the  definition  of  population 

variance  for  sampling  with  pps  and  replacement.  Notice  that, 

1  2 

when  the  P.  all  equal  a  defined  in  1.7  becomes  1.1. 

For  convenient  reference,  the  estimators  and  their  variances 
for  the  five  plans  to  be  discussed  are  listed  in  Table  1.1,  page 
29.  The  variances  are  expressed  as  population  values  (parameters) 
rather  than  as  sample  estimates  of  variance.  Each  variance 
formula  is  written  in  a  form  which  shows  a  sum  of  squares  of 
deviations  of  points  from  a  line  (or  lines).  Also,  an  alternative 


2/  A  good  reference  is:  Cochran,  W.G.,  Sampling  Techniques: 
Stratified  Random  Sampling,  Chapter  5;  Ratio  Estimates, 
Chapter  6;  Regression  Estimates,  Chapter  7;  and  for  sampling 
with  probability  proportional  to  size  see  Sections  9.9,  9.10, 
9.11,  and  9.12  of  Chapter  9. 
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expression  for  the  variance  of  the  estimator  for  each  plan  is 
shown.  For  simplicity,  an  assumption  is  made  that  the  sampling 
fractions  are  small  when  the  sampling  is  without  replacement. 
Thus,  the  fpc  (finite  population  correction)  factor  has  been 
omitted  from  the  variance  formulas.  The  fpc,  namely  pp,  can 
always  be  included  if  needed.  Notice  in  Table  1.1  that,  for 
a  constant  size  of  sample,  it  is  only  the  sums  of  squares  that 
differ  among  the  plans. 

A  dot  chart  that  shows  one  point  for  each  pair  of  values 
of  X^  and  provides  simple,  graphical  interpretations  of  the 
sums  of  squares  in  the  variance  formulas  for  the  five  plans. 
Each  variance  formula  for  the  first  four  plans  involves  the 
deviations  of  from  a  line  through  the  point  (X,Y).  The 
fifth  plan  involves  line  segments.  How  do  the  lines  for  the 
five  plans  differ  and  how  can  one  judge  the  sampling  variance 
for  one  plan  compared  to  another  by  looking  at  a  dot  chart? 
1.2.1  PLAN  1  -  MEAN  ESTIMATOR 


In  the  first  three  plans,  simple  random  sampling  is 

assumed.  These  three  plans  differ  only  with  regard  to  the 

method  of  estimating  Y.  The  first  plan  is  to  use  the  sample 
n 

Eyi 

average  y  =  — —  as  an  estimator  of  Y.  As  a  symbol  for  an 
estimator  we  will  use  y,  and  a  subscript  will  be  used  to 
distinguish  the  different  estimators.  Thus,  the  first  esti¬ 
mator  and  its  variance  are 


y- 


y 


n 


(1.8) 


7 


(1.9) 


vCyxD 


n 


s(Yi-Y) 

5Fl 


The  formula  for  the  variance  of  y^  contains  the  expres- 
N  -  2 

sion  z(Yi-Y)  .  As  shown  in  Figure  1.1,  the  vertical  distance 
between  a  point  (X^Y^  and  a  horizontal  line  through  (X,Y)  is 

N  -  2 

equal  to  (Y^-Y).  Hence  z(Y^-Y)  may  be  interpreted  as  the  sum 
of  squares  of  the  deviations  of  Y  from  a  horizontal  line  through 
(X,Y).  The  closer  the  points  are  to  this  horizontal  line,  the 

A 

smaller  the  variance  of  y^. 

In  the  general  context  of  regression  estimation,  Plan  1  is 
a  special  case.  Cochran,  in  Chapter  7,  Sampling  Techniques, 

A 

discusses  regression  estimation  where  y  in  the  following  equation 
is  the  regression  estimator: 

y  =  y  +  b  (X-x)  (1.10) 


The  value  of  the  regression  coefficient,  b,  might  be  preassigned 
or  it  might  be  computed  from  the  sample  data.  It  it  is  pre¬ 
assigned,  b  is  a  constant  when  one  considers  the  expected 
value  of  y.  If  b  is  constant,  it  is  clear  from  the  theory  of 

expected  values  that  E(y)  =  Y  because  the  expected  value  of  y 

-  -  3/ 

is  Y  and  the  expected  value  of  the  second  term,  b(X-x)  is  zero—  . 
Thus,  the  expected  value  of  y  is  Y  regardless  of  the  value  that 
is  preassigned  to  b.  There  are  cases  where  a  preassigned  value 
of  b  equal  to  1  is  of  interest  but  that  is  not  pertinent  to  the 


3 /  E[b(X-x)]  =  E(bX)  -  E(bx)  =  bX  -  bE(x)  =  0  because  E(x)  =  X  . 
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present  discussion.  The  point  of  interest  is  that  Plan  1  may¬ 
be  regarded  as  a  special  case  of  regression  estimation  where 
b  is  given  a  preassigned  value  equal  to  zero.  In  Plans  2  and 
3,  the  value  of  b  is  computed  from  the  sample. 

1.2.2  PLAN  2  -  RATIO  ESTIMATOR 

When  we  let  b  equal  ^  ,  the  right  side  of  Eq.  1.10 


becomes  X  ^  which  is  the  estimator  for  Plan  2.  Thus, 
x 

y,  -  x  i 


(l.u) 


This  estimator  is  called  a  ratio  estimator  since  it  is  the  ratio 
of  two  random  variables  y  and  x.  For  simple  random  sampling  the 

A 

variance  of  y2  is  often  written  as  follows: 


where 


v(y2)  -  ;p  -  (j-KSy  +  r2sx  '  2RSw] 


XY- 


N  .  2 
7  m.-Yr 

q^  =  x _ 

I  N-l 


Cl. 12) 


N 

2  S(X.-X) 

q  =  x 
^X  N-l 


SXY 


N 

Z(Xi-X)(Yi-Y) 


and 


R  = 


N 

ZY. 

i 

w~ 

EX. 

1 


Y 

X 
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The  variance  formula  for  Plan  2,  Table  1.1,  shows  that 
the  deviations,  (Y^-RX^),  are  squared  and  summed. 

Exercise  1.2  With  reference  to  the  variance  of  Eq.  1.12} 

2  2  2 

show  that  Sy  +  R  S-^  -  2RS^y  =  j- — — - 

Consider  a  line  through  the  origin  and  the  point  (X,Y) , 

Y 

see  Figure  1.1.  The  slope  of  this  line  is  R  =  —  .  The  vertical 

X 

distance  between  this  line  and  a  point  (X^,Y^)  is  (Y^-RX^). 

2 

Therefore,  the  sum  of  squares,  l(Y^-RX^)  ,  in  the  variance 
formula  for  is  the  sum  of  squares  of  the  deviations  of  the 

points  from  the  line  through  the  origin  and  (X,Y),  The 

A  A. 

only  difference  between  the  variances  of  y^  and  y^  is  the  dif- 

-  2  ? 

ference  between  Z(Y^-Y)  and  E(Y^-RX^)  .  The  points  for  the 

assumed  population  in  Figure  1.1  are  somewhat  closer  to  the 
line  through  the  origin  and  (X,Y)  than  to  a  horizontal  line 
through  (X,Y).  Therefore,  one  would  expect  y^  to  have  a 
smaller  sampling  variance  than  y^. 

Exercise  1.3  Verify  that  Y^-RX^  is  the  vertical  distance 
between  a  point  (X^,Y^)  and  a  straight  line  that  passes  through 
the  origin  and  (X,Y). 

1.2.3  PLAN  3  -  REGRESSION  ESTIMATOR 

A  * 

The  estimator,  y^,  in  Plan  3  is  called  a  regression  esti¬ 
mator.  It  makes  use  of  a  line  that  is  derived  by  applying  the 
least  squares  method  in  fitting  a  line  to  the  sample  values  of 
x  and  y.  The  equation  for  the  least  squares  line  (fitted  to 
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the  sample  data)  may  be  written  as  follows: 

yi  =  y  +  b (xi~  x)  (1.13) 

E  (x.-x) (y.-y) 

where  b  =  - * — ,  i  =  1,...,  n 

Z  (xi-x)Z 

and  y^  is  the  point  on  the  line  where  x  is  equal  to  x^.  The 
estimator  of  Y  is  obtained  by  substituting  X  for  x^  in  1.13 
which  gives 

y3  =  y  +  b(X-x)  (1.14) 

To  understand  the  variance  formula  for  y^,  suppose  a 
least  squares  line  is  determined  for  the  population  of  points 
shown  in  Figure  1.1.  It  is 

Y\  =  Y  +  B(Xi-X)  (1.15) 

E(X.-X) (Y.-Y) 

where  B  =  - * -  ,  i  =  1,...,  N. 

s(xi-x)z 

and  Y^  is  the  point  on  the  line  where  X  is  equal  to  X^.  This 
line  has  been  determined  so  the  sum  of  the  squares  of  the 

N  *  2 

deviations  of  Y^  from  it  is  a  minimum.  That  is,  z(Y^-Y^) 

is  less  than  the  sum  of  the  squares  of  the  deviations  from  any 
other  straight  line.  The  sum  of  squares  of  the  deviations  of  Y^ 
from  the  least-squares  regression  line  can  be  written  as 
follows : 

N  ?  N  7 

S(Yi-Yi)Z  =  E{Yi-[Y+B(Xi-X)]}^  (1.16) 
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The  expression  on  the  right  side  of  1.16  appears  in  the  variance 
formula  for  y^  in  Table  1.1  which  is 


2  N 

st  £{Y.  -  [Y+B (X .  -  X)  ]  }' 

V(y3)  -  -  (i)  — ! - 


N-l 


(1.17) 


Exercise  1.4  Show  that  the  right  side  of  1.16  reduces  to 

2  -2 

(1-r  ) E  (Yi~Y)  where  r  is  the  coefficient  of  correlation  between 
X  and  Y. 

1.2.4  DISCUSSION  OF  PLANS  1,  2,  and  3 

The  variances  of  y^,  y2>  and  y^  have  been  related  to  the 
sums  of  squares  of  deviations  from  three  lines  respectively: 

(1)  a  horizontal  line  through  (X,Y),  (2)  a  ratio  line  (that  is, 
a  line  through  the  origin  and  (X,Y),  and  (3)  a  regression  line 
(which  is  a  line  determined  by  the  method  of  least  squares). 
Since  the  sum  of  squares  of  deviations  from  the  regression  line 

A 

is  least,  the  variance  of  y^  will  generally  be  less  than  the 
variances  for  y^  and  y2»  The  comparative  variances  can  be 
judged  from  visual  examination  of  how  close  the  points  are  to 
each  of  the  three  lines. 

The  variance  of  y2  is  not  always  less  than  the  variance 

of  y^.  Moreover,  the  correlation  coefficient  is  not  a  reliable 

measure  of  how  the  variances  of  y^  and  y2  compare.  According  to 

2  2 

Eq.  1.12,  2RS^y  must  be  larger  than  R  S^  or  the  variance  of  y2 


will  be  larger  than  the  variance  of  y^.  In  other  words,  use 
of  an  auxiliary  variable  in  a  ratio  estimator  could  result  in 
an  increase  rather  than  a  decrease  in  variance. 
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The  variance  formulas  discussed  above  are  population 
variances  (parameters)  which  must  be  estimated  from  the  sample. 
For  all  three  plans,  formula  for  estimating  the  sampling  vari¬ 
ances  are  of  the  same  format  as  the  population  variance  formula. 
The  only  difference  is  that  the  sum  of  squares  is  computed  from 
sample  data  instead  of  data  for  the  entire  population.  The 
variance  formulas  for  Plans  2  and  3  are  large  sample  approxi¬ 
mations,  which  are  commonly  used  in  practice.  (See  Cochran's 
book  sections  6.4  and  7.4.) 

In  a  survey  involving  many  variables  and  tabulations  by 
various  classifications,  the  first  two  estimators  (plans)  are 
commonly  used.  Although  the  variance  of  y^  is,  to  some  degree, 

A  A 

generally  less  than  the  variance  of  y^  or  y^ ,  its  use  is  gen¬ 
erally  limited  to  special  situations  where  low  error  is  very 
important  and  the  variance  of  y^  is  appreciably  less  than  the 
variance  of  y^  or  y  For  example,  it  might  be  used  to  esti¬ 
mate  the  production  of  a  particular  commodity  or  when  it  is 
very  important  to  make  estimates  with  a  high  degree  of  accuracy 
for  a  few  selected  characteristics. 

All  three  of  the  estimators  may  be  used  with  sampling 
plans  other  than  simple  random  sampling;  for  example,  ratio 
estimators  and  stratified  random  sampling  are  quite  common. 

Exercise  1.5  For  the  special  case  where  the  regression 
line  is  the  same  as  the  ratio  line ,  show  that  the  variance 
of  y^  is  equal  to  the  variance  of  y Can  V(y^)  ever  be 
larger  than  V(y~)? 
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Exercise  1.6  Compare  Elans  1}  23  and  3  with  regard  to 
the  following  three  dot  charts  representing  three  different 
relations  between  X  and  Y. 


Y 


Y 


Case  2 


Case  3 


X 


For  each  of  the  three  cases  rank  the  three  plans  from  largest 
to  smallest  sampling  variance . 

1.2.5  PLAN  4  -  SAMPLING  WITH  PPS 

Plans  2  and  3  used  the  auxiliary  variable  in  estimation 
and  not  in  the  design  or  selection  of  a  sample.  Plan  4  is  to 
select  a  sample  of  n  elements  with  replacement  and  to  use 
probabilities  of  selection  proportional  to  X^.  By  substituting 
x.  X. 

^ -  for  p.  in  Eq.  1.5  and  — for  P.  in  1.7,  the  following  ex- 

A  1  A  1 

pressions  are  obtained  for  the  estimator  and  its  variance: 


-  1  n  yi 
y,  =  X(-)  E  — 

(1.18) 

i 

and 

v(y4)  -  ^  -  (i)(i)x(|-)(Yi-Rxi)2 

(1.19) 

The  formula  for  the  variance  of  y4  shows  that  (Y^-RX^) 
are  the  deviations  which  are  squared.  Thus,  the  line  involved 
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in  Plan  4  is  the  same  as  the  line  for  the  ratio  estimator. 

2 

Notice  that  the  squares  of  the  deviations,  (Y^-RX^)  ,  are 

weighted  by  owing  to  the  unequal  probability  of  selection. 

1 

For  the  ratio  estimator,  the  squares  of  the  deviations  were 

weighted  equally.  Incidentally,  the  appropriate  formula  for 

estimating  the  variance  of  y^  from  sample  data  is  not  of  the 

same  form  (and  will  not  reduce  to  the  same  form)  as  Eq.  1.19. 

In  practice  one  often  finds  that  the  variance  of  the 

deviations,  (Y^-RX^) ,  increases  as  X  increases.  That  is,  the 

values  of  Y  are  usually  more  widely  scattered  for  large  values 

of  X  than  for  small  values  of  X.  If  the  relation  between  X  and 

Y  is  like  the  dot  chart  in  Figure  1.2,  Plan  4  will  have  a  lower 

sampling  variance  than  the  first  three  plans.  A  line  through 

(X,Y)  and  the  origin  fits  the  data  about  as  well  as  any  line. 

But,  y^  would  have  the  least  sampling  variance  because,  as  shown 

2 

in  the  formula  for  its  variance,  the  largest  values  of  (Y^-RX^) 
receive  the  smallest  weights  in  the  sum  of  squares.  Judging  the 
effectiveness  of  Plan  4  is  more  than  a  matter  of  observing  how 

well  the  data  fit  a  line  through  (X,Y)  and  the  origin.  In  fact , 

it  is  easy  to  misjudge  the  effectiveness  of  sampling  with  pps . 

We  will  return  to  this  point  after  presentation  of  Plan  5. 

2 

Exercise  1.7  Start  with  a  as  defined  in  1.7  and  show 

N  -  X 

1  IN  x  2  A  i  Y 

that  it  reduces  to  (^)  £  ^r-(Y^-RX^)  when  P^  =  and  R  =  ^  • 

1.2.6  PLAN  5  -  STRATIFIED  SAMPLING 

This  plan  makes  use  of  the  variable  X  as  a  basis  for 
stratification.  Suppose  the  sampling  units  in  the  population 
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have  been  listed  in  order  from  smallest  to  largest  values  of  X. 

The  list  is  then  divided  into  L  strata.  Let 

=  the  population  number  of  su's  in  stratum  h, 

n^  =  the  sample  number  of  su’s, 

n. 

f^  =  jq—  =  the  sampling  fraction, 
h 

Y,  •  and  X,  .  =  the  values  of  Y  and  X  for  the  i**1 
hi  hi 

su  in  stratum  h, 

2 

SYh  =  the  variance  of  Y  within  stratum  h, 

Y^  =  the  average  value  of  Y  in  stratum  h,  and 
X^  =  the  average  value  of  X  in  stratum  h. 

We  are  primarily  interested  in  proportional  allocation  of  the 
sample  to  strata  for  comparison  with  Plans  1,  2,  and  3,  and  in 
optimum  allocation  for  comparison  with  Plan  4. 

With  proportional  allocation  the  sampling  fractions,  f^, 
are  all  equal  and  it  is  appropriate  to  use  the  unweighted  sample 
mean  as  an  estimator  of  Y.  Hence, 


Yr  =  y 


(1. 20) 


Assuming  simple  random  sampling  within  strata  and  that  the 
fpc's  are  negligible, 


V(y5)  = 

S2 

_5 

n 

where 

-  — 
N 

lNhSYh 

and 

= 

^Yhi-V2 

i  h 

V1 

(1.21) 
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With  reference  to  a  dot  chart  for  showing  deviations  that 
are  squared  in  the  variance  formulas,  instead  of  one  line,  we 
now  have  a  series  of  line  segments,  one  for  each  stratum,  as 
shown  in  Figure  1.3.  Each  line  segment  is  a  horizontal  line 
through  the  stratum  mean.  The  sampling  variance,  S^,  is  an 
average  of  the  squares  of  deviations  from  these  horizontal  line 
segments.  If  the  points  are  close  to  the  line  segments,  the 
sampling  variance  will  be  small  for  stratified  random  sampling. 

Consider  what  happens  to  the  sum  of  squares  for  stratified 
random  sampling  as  the  number  of  strata  increases,  that  is, 
as  the  difference  between  the  largest  and  smallest  value  of 
X  for  each  stratum  decreases.  If  the  relation  between  X  and  Y 
over  the  whole  population  is  approximately  linear,  the  sum  of 
squares  of  the  deviations  from  the  line  segments  will  become 
approximately  equal  to  the  sum  of  squares  of  the  deviation  from 
a  regression  line  as  in  Plan  3.  Under  those  conditions  Plans  3 
and  5  would  have  approximately  the  same  sampling  variance.  If 
the  relation  between  X  and  Y  is  not  linear,  the  sampling  variance 
for  Plan  5  might  be  less  than  the  sampling  variance  for  Plan  3, 
depending  on  the  width  of  the  stratum  intervals,  the  degree  of 
nonlinearity,  and  how  close  the  points  are  to  a  curved  line. 

Suppose  the  ratio  line  (that  is,  a  straight  line  through 
(X,Y)  and  the  origin)  fits  the  points  about  as  well  as  any  line. 
In  this  case,  the  sampling  variances  for  Plans  2,  3  and  5 
(assuming  the  stratum  intervals  are  small)  would  be  approximately 
equal . 
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Plan  4  must  be  judged  with  regard  to  how  well  the  prob¬ 
abilities  of  selection  fit  the  situation  as  well  as  the  close¬ 
ness  of  the  points  to  the  ratio  line.  It  is  helpful  to  compare 
it  with  using  the  auxiliary  variable  for  stratification  and 
optimum  allocation  of  the  sample  to  strata.  We  know  that  the 
optimum  size  of  sample  from  stratum  h  is  proportional  to  N^Sy^. 

Or,  in  terms  of  sampling  fractions,  the  optimum  sampling 
fraction,  f^,  is  proportional  to  Sy^. 

In  stratified  sampling,  the  optimum  sampling  fractions  are 
proportional  to  X^  when  Sy^  is  proportional  to  X^.  In  this  case, 
the  selection  probabilities  in  sampling  with  pps  would  be  approxi¬ 
mately  in  proportion  to  X^  provided  the  stratum  intervals  are 
small.  In  other  words,  when  Sy^  is  proportional  to  X^  the 
optimum  sampling  fractions  in  stratified  sampling  are  in  close 
agreement  with  the  selection  probabilities  in  sampling  with 
pps.  It  is  very  important  to  recognize  that  the  situation  most 
favorable  for  sampling  with  pps  occurs  when  (1)  the  data  follow 
the  ratio  line,  and  (2)  the  conditional  standard  deviation  of 
Y  is  proportional  to  X.  ("Conditional  standard  deviation" 
refers  to  the  standard  deviation  of  Y  for  a  given  value  of  X.) 

The  dot  chart.  Figure  1.2,  meets  those  conditions.  Notice  that 
the  vertical  distance  between  the  two  dotted  lines  is  proportional 
to  X;  hence,  the  conditional  standard  deviation  of  Y  is,  at  least 
roughly,  proportional  to  X. 

Recognition  of  a  relation  like  the  one  in  Figure  1.2  as  a 
good  case  for  sampling  with  pps  provides  guidance  when  making  a 
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choice  among  alternatives  including  the  possibility  of  making 
a  transformation  of  X  that  would  provide  a  better  measure  of 
size.  Sometimes  a  simple  transformation  like  X^  =  X^  +  c > 
where  C  is  a  constant,  will  provide  a  measure  of  size,  X', 
such  that  the  conditional  standard  deviation  of  Y  will  be  in 
proportion  to  X'.  In  some  cases,  a  simple  transformation  can 
change  sampling  with  pps,  compared  to  Plan  1,  from  a  substantial 
increase  in  sampling  variance  to  an  important  reduction.  With 
pps  sampling  it  is  important  that  the  maximum  values  of  Y  approach 
zero  as  X  approaches  zero. 

One  might  feel  that  sampling  with  probability  proportional 
to  X  does  not  fully  remove,  from  the  sampling  variance,  variation 
among  strata  when  X  is  the  criterion  for  stratification.  Look 
at  the  pps  estimator.  It  is  variation  in  the  stratum  ratios, 

1  Y,. 

R,  =  Tf—  Z  tt— -  rather  than  variation  in  Y,  that  needs  to  be 
h  Nh  i  Xh.  h 

considered.  It  will  be  easier  to  discuss  this  point  in  the 
next  chapter  when  stratification  in  combination  with  different 
methods  of  estimation  is  considered. 

Some  numerical  results  as  well  as  dot  charts  will  be 
presented  later  in  this  chapter  and  in  Chapter  II. 

Exercise  1.8  Refer  to  Figure  1.2  and  verify  from  theorems 
pertaining  to  similar  triangles  that  the  range  in  values  of  Y 
is  proportional  to  X.  In  this  case3  as  a  rough  approximation , 
we  may  regard  the  standard  deviation  of  Y  as  being  in  proportion 
to  X.  Is  it  possible  in  sampling  with  probability  proportional 


19 


to  X  to  have  a  lower  sampling  variance  than  sampling  with 
stratification  by  X  and  optimum  allocation?  When? 

Exercise  1.9  (a)  Refer  to  Figure  1.1  and  rank  Plans  13 

23  3S  and  5  from  least  variance  to  highest.  Ans .  3S  53  23  1 
with  3  and  5  being  close  depending  on  the  number  of  strata. 

(b)  It  appears  that  the  variance  for  Plan  4  would  be 
much  larger  than  the  variance  for  stratified  random  sampling 
with  optimum  allocation.  Why?  Look  at  the  conditional 
standard  deviation  of  Y. 

(c)  Since  the  range  in  the  optimum  sampling  fractions 
for  stratified  random  sampling  is  small3  would  you  agree  that 
Plan  4  would  have  a  much  larger  variance  than  Plan  1? 

(d)  Consider  the  simple  transformation  XT  =  X^  +  C  where 
C  is  a  constant .  Is  there  a  value  of  C  such  that  X'  would  be 
an  effective  measure  of  size. 

Exercise  1.10  Refer  to  Exercise  1.6  and  for  each  case 
rank  all  five  plans  with  regard  to  sampling  variance . 

Exercise  1.11  Prepare  a  dot  chart  showing  a  relation 
between  X  and  Y  such  that  stratified  random  sampling  with 
allocation  proportional  to  N^,  Plan  53  will  have  a  smaller 
sampling  variance  than  the  regression  e stimator 3  Plan  3. 

Exercise  1.12  Prepare  a  dot  chart  such  that  the  variance 
for  Plan  5  with  proportional  allocation  will  be  approximately 
equal  to  the  variance  for  Plan  1  and  (at  the  same  time)  the 
variance  for  Plan  5 3  with  optimum  allocation  will  be  much 
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less  than  the  variance  for  Plan  1.  This  would  he  a  case  where 
gain  from  stratification  would  he  entirely  attrihutahle  to 
varying  sampling  fractions  rather  than  stratification  to  remove 
variation  associated  with  differences  among  stratum  means. 

1.2.7  SUMMARY 

If  there  is  no  relation  between  X  and  Y,  including  a 
relation  between  X  and  the  conditional  standard  deviation  of 
Y,  information  about  X  offers  no  possibilities  for  reducing 
sampling  variance;  in  fact,  the  sampling  variance  could  be  in¬ 
creased  by  using  X.  If  there  is  a  relation,  some  alternative 
ways  to  take  advantage^  of  it  have  been  shown.  Clearly,  the 
most  effective  way  of  using  an  auxiliary  variable  depends  on 
what  the  relation  is  like. 

In  the  sampling  and  estimation  specifications  for  a 
particular  survey,  an  auxiliary  variable  would  generally  be 
used  in  only  one  way.  For  example,  attempting  to  use  a  relation¬ 
ship  between  X  and  Y  as  a  basis  for  stratification  and  also  in 
estimation  is  generally  not  advisable.  Try  to  fully  utilize 
the  potential  contribution  of  an  auxiliary  variable  in  one 
way.  Whether  an  auxiliary  variable  is  used  in  stratification 
or  in  estimation  might  depend  on  the  nature  of  other  auxiliary 
variables  that  are  available.  For  example,  some  kinds  of 
auxiliary  variables  are  readily  useful  in  stratification  but 
not  estimation.  Consider  using  quantitative  measures  in  esti¬ 
mation  or  in  sampling  with  pps  and  using  nonquant itat ive  measures 
in  stratification.  This  point  will  receive  further  attention. 
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1.3  NUMERICAL  EXAMPLE 


Although  our  interest  is  in  the  practical  application  of 
sampling  theory,  a  major  objective  in  the  presentation  of 
numerical  illustrations  in  this  and  later  chapters  is  to 
improve  one's  comprehension  of  patterns  of  variation  that 
exist  and  to  develop  one's  skill  at  judging  the  effectiveness 
of  alternative  sampling  and  estimation  methods  in  specific 
situations.  It  is  informative  to  apply  several  alternatives 
to  the  same  population  even  though  some  of  the  alternatives 
are  not  practically  feasible. 

The  data  for  the  following  example  were  taken  from  a 
research  project  to  develop  techniques  for  sampling  apple  trees 
to  forecast  and  estimate  apple  production.  The  primary  purpose 
was  to  make  an  intensive  investigation  of  ways  of  sampling  a 
tree  rather  than  how  to  select  a  sample  of  trees.  As  a  part 
of  this  project,  the  branches  on  six  apple  trees  were  mapped. 
Included  among  the  measurements  that  were  taken  are  the  cross- 
sectional  area  of  each  branch  and  the  number  of  apples  on  each 
branch.  There  was  a  total  of  28  primary  branches  on  the  six 
trees.  A  primary  branch,  which  is  a  branch  from  the  tree  trunk, 
probably  would  not  be  used  as  a  sampling  unit  in  practice.  How¬ 
ever,  data  for  these  28  primary  branches  are  useful  as  a 
numerical  example  of  alternative  ways  of  using  an  auxiliary 
variable.  Also  the  results  will  be  useful  in  later  discussions 
and  comparisons  of  methods  of  sampling  within  trees. 
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For  purposes  of  this  numerical  example,  the  28  primary- 
branches  is  the  population  of  sampling  units.  We  assume  the 
purpose  of  sampling  is  to  estimate  the  total  number  of  apples 
on  the  six  trees.  The  auxiliary  variable  X  is  the  csa  (cross- 
sectional  area)  of  a  branch.  The  fruit  counts,  Y,  and  the 
csa’s,  X,  for  the  28  limbs  are  presented  in  Table  1.2.  Let 
us  compare  the  five  plans  outlined  above  by  referring  to  a  dot 
chart.  Figure  1.4  shows  the  points  (X^,Y^)  and  three  lines: 

(1)  the  horizontal  line  for  Plan  1,  (2)  a  ratio  line  through 
the  origin  and  (X,Y)  which  pertains  to  Plans  2  and  4,  and 
(3)  the  least  squares  regression  line  for  Plan  3.  To  order 
the  sampling  variances  from  smallest  to  largest,  one  would 
undoubtedly  rank  the  first  three  plans  in  the  order  3,  2,  and 
1,  with  1  having  a  much  larger  variance  than  the  other  two. 

Since  the  scatter  of  the  points  increases  as  the  csa  increases, 
one  might  expect  Plan  4  to  be  better  than  Plan  2,  but  Plan  4 
is  somewhat  difficult  to  judge.  In  Chapter  II,  similar  com¬ 
parisons  of  the  plans  will  be  made  using  terminal  branches 
(and  hence  more  points)  as  sampling  units. 

The  total  number  of  sampling  units,  28,  is  too  small  to 
provide  a  good  example  of  stratified  random  sampling  in  com¬ 
parison  to  the  other  four  plans.  However,  for  purposes  of 
illustration,  a  comparison  will  be  made.  Since  28  is  divisible 
by  4,  it  is  convenient  to  divide  the  branches,  after  being 
ordered  by  csa,  into  four  strata  of  7  branches  each  as  presented 
in  Table  1.2. 
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The  stratum  boundaries  are  indicated  by  vertical  dotted 
lines  in  Figure  1.4.  It  is  evident  that  line  segments,  for 
the  stratified  random  sampling  as  specified  in  the  preceding 
paragraph,  do  not  fit  the  data  as  well  as  the  regression  line. 
Plan  3.  Although  the  sampling  variance  for  Plan  5  is  clearly 
much  less  than  the  variance  for  Plan  1 ,  it  is  undoubtedly  greater 

than  the  variance  for  Plan  3.  Its  rank  compared  to  Plans  2  and 
4  is  uncertain. 

We  will  now  compare  the  judgments  formed  from  looking  at 
Figure  1.4  with  numerical  results.  The  relative  variances  of 
the  five  estimators,  assuming  n  =  1  (that  is,  a  sample  of  one 
branch),  are  presented  in  Table  1.3.  Relative  variances  are  the 
variances  divided  by  Y2.  Although  it  is  not  possible  to  select 
a  stratified  random  sample  of  one  branch,  it  is  appropriate  to 
let  for  purposes  of  comparing  Plan  5  with  the  other  plans. 

In  this  example,  the  relationship  between  X  and  Y  is  such 
that  all  four  Plans  2,  3,  4,  and  5  provide  large  reductions  in 
sampling  variance.  Stratification,  as  applied,  reduced  the 
sampling  variance  by  more  than  80  percent  compared  to  Plan  1 
but  not  as  much  as  Plans  2,  3,  and  4  because  it  did  not  utilize 
as  fully  the  information  provided  by  X.  If  it  were  feasible  to 
divide  the  population  into  more  strata,  perhaps  8  or  10  instead 
of  4,  the  relative  variance  for  Plan  5  would  have  been  less  than 
0.307  and  perhaps  nearly  as  low  as  the  variance  for  the  regression 
estimator,  Plan  3.  However,  from  the  results  that  we  have  seen, 
at  appears  that  the  auxiliary  variable  X  can  be  used  to  reduce 
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the  sampling  variance  from  1.117  to  about  0.200.  Some  of  the 
practical  considerations  in  the  choice  of  a  plan  will  be  dis¬ 
cussed  later.  In  the  next  section  our  understanding  of  sampling 
with  pps  will  be  extended  by  comparing  it  to  stratification  with 
optimum  allocation. 

1.3.1  VARYING  THE  SAMPLING  FRACTION  WITH  SIZE  OF  SAMPLING  UNIT 
From  Figure  1.4  it  is  clear  that  the  variance  of  the 
number  of  apples  increases  with  the  size  of  branch.  The  stan¬ 
dard  deviation  within  strata  and  the  average  csa  per  branch  are 
presented  in  Table  1.4. 

Since  the  largest  Sy^  is  about  10  times  larger  than  the 
smallest,  the  largest  sampling  fraction  (with  stratified 
sampling  and  optimum  allocation)  would  be  about  10  times 
larger  than  the  smallest.  This  range  of  variation  in  sampling 
fractions  is  large  enough  to  expect  optimum  allocation,  compared 
to  proportional,  to  give  a  substantial  reduction  in  variance. 

The  relative  variance  for  optimum  is  0.211  compared  to  0.301 
for  proportional. 

With  reference  to  sampling  with  pps,  notice  that  the 
conditional  standard  deviation  of  Y  is  roughly  in  proportion 
to  X.  This  is  indicated  by  the  fact  that  the  ratio  o  f  Syh  to 
X^,  Table  1.4,  is  nearly  constant.  Also,  the  points  in  Figure  1.4 
follow,  approximately,  a  line  through  the  origin  and  (X,Y).  There¬ 
fore,  it  is  reasonable  to  find  that  0.211,  the  variance  for 
stratified  sampling  with  optimum  allocation,  is  close  to  0.194, 
the  variance  for  sampling  with  pps. 
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Since  is  approximately  in  proportion  to  X^,  csa  is  a 

good  measure  of  size.  However,  it  is  informative  to  compare 
the  five  plans  when  circumference  is  used  or  a  measure  of  size 
of  branch.  To  examine  the  relation  between  number  of  apples 
and  circumference,  see  Figure  1.5.  Notice  that  the  least  squares 
line  (Plan  3)  departs  farther  from  the  origin  than  did  the  least 
squares  line  for  csa.  Figure  1.4.  This  is  reflected  in  the 
variances  which  are  presented  in  Table  1.5.  The  relative  vari¬ 
ance,  0.256,  for  Plan  3  is  considerably  less  than  the  relative 
variances  for  Plans  2  and  4.  Also  notice  that  circumference 
is  less  effective  than  csa  for  all  three  Plans  2,  3,  and  4. 

Exercise  1.13  Refer  to  Table  1.2  and  compute  the  four 
values  of  taking  the  circumference  as  the  auxiliary  variable . 
Compare  these  values  of  X^  with  the  values  of  given  in 

Table  1.4.  What  does  this  comparison  indicate  regarding  the 
use  of  circumference  as  a  measure  of  size  in  pps  sampling? 

Notice  that  csa  is  a  mathematical  transformation  of  circum¬ 
ference.  The  question  might  be  asked,  MIs  there  a  better 
transformation?”  This  question  will  be  given  further  attention 
in  the  next  chapter.  For  the  research  study,  a  csa  measurement 
was  made  by  wrapping  a  tape  around  the  base  of  a  branch.  The 
tapes  had  been  calibrated  to  give  a  direct  reading  of  the  csa 
assuming  the  branch  is  circular.  Figure  1.4  suggests  that  csa 
is  a  good  measure  of  size  for  sampling  with  pps,  but  broader 
experience  is  needed.  In  a  later  illustration  it  will  become 
evident  that  sampling  with  pps  is  a  good  practical  method  of 
selecting  a  sample  of  branches. 
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Exercise  1.14  By  careful  planning  one  can  compute  sub- 


2  Y- 

totals  and  totals  of  ZY.  ,  ZX .  ,  ZY .  ,  ZX.Y.,  and  Z— that 

J  l  l  l  li  X. 

l 

provide  intermediate  results  from  which  the  variances  for 

several  alternative  plans  are  easily  obtained.  For  purposes 

2 

of  computation  show  that  the  values  of  S  for  the  five  plans 
may  be  written  as  follows : 

2  1  2  (ZYi)2 

51  ■  -  -rr-3 

ZY . 

52  =  (xfirHz Y?  -  2RZX.Y.  +  R2ZX?]  where  R  = 

2  N- 1J  L  l  li  l J  £X^ 

2  2  2 

=  S^(l-r  )  where  r  is  the  correlation  coefficient 


(— ")  j — i.  -  y 

l 


c,2 


Since  there  are  7  branches  in  each  stratum  the  expression 
2 

in  Table  1.1  for  reduces  to 

ZY2 

[ ZY2  -  — y-^-1  where  Y^  is  the  total  of  Y  for 

stratum  h. 

From  Table  1.2  the  following  intermediate  results  are 
obtained : 


zYi  ' 

7,199 

zXi 

=  157.76 

v2 

zYi  ' 

3,844,283 

zx2i 

=  1,329.98 

ZX.Y.  = 

67,633. 47 

Y2 

1  1 

eX7 

=  392,247.3 
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The  stratum  totals ,  Y^j  are 

Y1  =  202,  Y2  =  923,  Y 3  =  1,594,  and  Y 4  =  4,480 
2  2  2  2  2 

Compute  the  values  of  S^,  S2>  S“,  a^,  and  . 


Answer: 


1 

,2 

>2 

.2 

*3 

2 

f4 

.2 


73,828 
16, 339 
12,292 
12 , 826 
20,274 
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Table  1.2 — Data  for  Primary  Branches  on  Six  Apple  Trees 

(Arrayed  by  csa) 


Stratum 

Branch^ 

2 / 
csa-7 

3/;No.  of 
uir‘  '.apples 

Stratum 

Branch^/ 

1  2/ 

3/:tJo.  of 
; apples 

1 

1-4 

.87 

3.3 

5 

3 

6-3 

4.84 

7.8 

183 

1-5 

1.03 

3.6 

34 

1-2 

5.09 

8.0 

40 

5-6 

1.34 

4.1 

4 

2-4 

5.75 

8.5 

396 

1-3 

1.83 

4.8 

59 

6-2 

5.89 

8.6 

250 

5-4 

1.83 

4.8 

18 

2-3 

6.16 

8.8 

157 

5-5 

1.83 

4.8 

17 

5-1 

6.16 

8.8 

179 

6-4 

1.99 

5.0 

65 

4-2 

7.18 

9.5 

389 

2 

2-5 

2.68 

5.8 

89 

4 

2-2 

8.94 

10.6 

333 

4-4 

2.86 

6.0 

238 

4-1 

9.28 

10.8 

696 

4-5 

2.86 

6.0 

81 

2-1 

9.63 

11.0 

473 

4-3 

3.57 

6.7 

254 

3-1 

11.60 

12.1 

762 

1-1 

3.68 

6.8 

76 

3-3 

12.84 

12.7 

517 

5-3 

4.48 

7.5 

97 

3-2 

13.45 

13.0 

622 

5-2 

4.72 

7.7 

88 

6-1 

15.38 

13.9 

1,077 

TOTAL 

157.76 

221.0 

7,199 

1/  Tree  (first  digit)  and  branch  within  a  tree  (second  digit). 
2/  Cross-sectional  area  of  branch  in  square  inches. 

3/  Circumference  of  branch  in  inches. 


Table  1.3 — Relative  Variances  of  Estimators 


A 

Plan 

Relative  variance  of  y 

1 

1 . 117 

2 

0.247 

3 

0.186 

4 

0.194 

5 

0.307 
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Table  1.4--Mean  csa  and  Standard  Deviation  of  Y  by  Strata 


Stratum 

Mean  csa, 

Standard 

Deviation , 

SYh 

SYh 

*h 

1 

1.53 

24.8 

16.2 

2 

3.55 

78.4 

22.1 

3 

5.87 

129 

22.0 

4 

11.59 

240 

20.7 

Table  1 . 5- -Relative  Variances  When  the  Auxiliary  Vari¬ 
able  is  Circumference 


Plan 


Relative  variance  of  y 


1  1.117 

2  0.559 


3 


0.256 


4 


0.438 
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Figure  1.1 — Deviations  in  Variance  Formulas  for 
Plans  1,  2,  and  3 
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Figure  1.2 — Deviations  in  Variance  Formulas  for 
Plans  2  and  4 
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*(X1,Y^ 


Stratum  1  Stratum  2  Stratum  3  Stratum  4 


Stratum  5 


X 


Figure 

for 


1.3 — Deviations  in  Variance  Formula 
Plan  5,  Stratified  Random  Sampling 
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Figure  1.4 — Relation  between  Number  of  Apples  and  CSA 
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Figure  1.5 — Relation  between  Number  of  Apples  and  Circumference 


CHAPTER  II 


FURTHER  OBSERVATIONS  ON  USES  OF  AN  AUXILIARY  VARIABLE 

2.1  INTRODUCTION 

The  effects  on  sampling  variance  of  various  factors  in 
sample  design  and  estimation  are  not  independent.  For  example, 
the  difference  in  the  sampling  variance  between  a  mean  esti¬ 
mator  and  a  ratio  estimator  might  vary  with  the  definition  of 
the  sampling  unit  or  with  the  criteria  used  for  stratification. 
In  this  chapter  some  numerical  examples  that  display  such  inter¬ 
actions  will  be  given.  The  objective  is  to  further  develop  a 
perception  of  patterns  (or  components)  of  variation  and  ability 
to  judge  how  alternative  methods  rank  with  regard  to  sampling 
variance.  As  you  study  and  acquire  experience  in  sampling  try 
to  visualize  the  pattern  of  variation  in  a  population  to  be 
sampled  and  test  your  skill  at  prejudging  the  effectiveness  of 
alternative  sampling  plans. 

The  data  for  the  examples  in  this  chapter  are  taken  from 
the  research  project  on  methods  of  estimating  apple  production 
which  was  referred  to  in  Chapter  I.  The  sampling  alternatives 
that  are  considered  require  a  map  of  each  tree  that  is  sampled. 
That  is,  a  map  of  a  tree  which  defines  the  sampling  units 
(branches)  is  the  sampling  frame.  Methods  of  probability 
sampling  are  available  which  do  not  require  preparing  a  com¬ 
plete  map  of  a  tree.  This  will  be  discussed  in  Chapter  III. 
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As  background,  refer  to  Figure  2.1  which  is  a  map  of  one 
of  the  six  apple  trees  used  for  the  numerical  example  presented 
in  Chapter  I.  The  map  shows  the  scheme  that  was  used  for  iden¬ 
tifying  branches.  For  example,  3-1-4  refers  to  third-stage 
branch  number  4  from  second-stage  branch  number  1  and  first- 
stage  branch  number  3.  Branches  from  the  tree  trunk  were 
mapped  until  ’'terminal”  branches  were  reached.  "Terminal 
branch”  refers  to  the  last  stage  of  branching  where  the  mapping 
of  branches  was  terminated.  The  csa's  (cross  sectional  areas) 
of  the  terminal  branches  ranged  from  about  3/4  to  2  square  inches 
which  seemed  to  be  about  the  smallest  practical  size  of  branch 
to  consider  as  a  sampling  unit.  There  were  28  primary  branches 
and  135  terminal  branches  on  the  six  trees.  The  average  number 
of  apples  on  a  terminal  branch  was  about  50. 

When  following  a  tree  trunk  to  primary  branches,  to  second- 
stage  branches,  etc.,  small  branches  are  sometimes  found  which 
are  not  large  enough  to  be  classified  as  terminal  branches. 

For  example,  six  apples  were  found  on  small  branches  on  primary 
branch  number  2  before  the  4  second-stage  branches  2-1,  2-2,  2-3, 
and  2-4  were  reached.  Apples  on  such  branches  have  been  called 
"path”  fruit,  meaning  fruit  on  the  path  of  a  terminal  branch. 

Path  fruit  present  some  special  problems  which  will  be  discussed 
in  Chapter  III.  The  amount  of  path  fruit  is  relatively  small 
and  will  be  ignored  in  this  chapter. 

For  each  of  the  first  four  plans  that  were  discussed  in 
Chapter  I  primary  and  terminal  branches  will  be  compared  as 
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sampling  units.  Then,  using  terminal  branches,  the  first  four 
plans  will  then  be  applied  within  strata  (trees)  for  comparison 
with  each  of  the  four  plans  when  there  is  no  stratification. 

2.2  COMPARISON  OF  PRIMARY  AND  TERMINAL  BRANCHES  AS  SAMPLING 

UNITS 

The  number  of  applies  on  each  of  the  28  primary  branches 
and  the  csa  of  each  branch  were  presented  in  Table  1.2.  Data 
for  the  135  terminal  branches  are  presented  in  Table  2.1.  The 
number  of  apples  on  primary  branches  included  path  fruit  whereas 
the  numbers  on  terminal  branches  do  not.  The  difference  is  pre¬ 
sumed  to  be  negligible  for  purposes  of  an  exercise  in  variance 
comparisons.  Figures  1.4  and  2.2  are  the  dot  charts  for  primary 
and  terminal  limbs  respectively. 

Table  2.2  presents  relative  variances  for  terminal  and 
primary  branches.  The  relative  variances  for  primary  branches 
are  taken  from  Table  1.3  in  Chapter  I,  and  relative  variances 
for  terminal  branches  were  computed  using  the  same  variance 
formulas . 

When  interpreting  variances  it  is  essential  that  the 
dimensions  of  the  variances  be  clear.  What  variation  does 
a  particular  variance  measure  and  in  what  units  is  the  variance 
expressed?  Are  the  relative  variances  in  Table  2.2  comparable? 
Let  us  examine  the  formula  for  the  relative  variance  (RV)  of 
y^ ,  which  is 

RelVar  (yp  =  ^(i)  (S^)  (2.1) 
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where 


S 


2 

1 


£  (Yi  -  Y)  2 
N-l 


2 

A  quantity  like  is  sometimes  called  "unit  variance"  as  it 
is  a  measure  of  variation  among  individual  sampling  units.  The 


quantity  — may  be  called  "unit-relative  variance"  which  is 

Y2 


the  square  of  the  coefficient  of  variation  among  individual 

units.  In  Eq.  2.1,  when  n  =  1  the  relative  variance  of  y^, 

is  the  unit  -  re lat ive  variance.  A  similar  interpretation  of 

the  variance  formula  for  the  other  estimators  holds.  Thus, 

2 

S2  is  the  unit  variance  that  pertains  to  the  ratio  estimator, 


The  variances  presented  in  Table  2.2  are  unit- relative 
variances  which  may  be  regarded  as  sampling  variances  for 
samples  of  one  branch.  Usually  sampling  variances  for  alterna¬ 
tive  plans  are  compared  under  one  of  two  conditions:  equal 
sampling  fractions  or  equal  costs.  In  this  chapter  the  com¬ 
parisons  will  be  under  an  assumption  of  equal  sampling  fractions. 
The  sampling  fractions  are  and  yjEj-,  respectively,  for  one 
primary  branch  and  one  terminal  branch. 

To  achieve  comparability,  the  variances  for  primary  branches 
will  be  converted  to  the  equivalent  of  one  terminal  branch.  That 
is,  we  want  to  find  the  variances  for  primary  branches  that 
correspond  to  a  sampling  fraction  of  y-yy-.  There  is  an  average 
of  — yg-  =  4.82  terminal  branches  per  primary  branch. 
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Ignoring  the  fpc  (finite  population  correction),  the 

a  117 

relative  variance  of  the  first  estimator,  y,  ,  is  -  for  a 
sample  of  n^  primary  branches  and  is  --  for  a  sample  of  n 

terminal  branches.  (The  numbers,  1.17  and  0.660,  are  from 
Table  2.2.)  Since  the  sampling  fractions  are  the  same  for 
terminal  and  primary  branches  when  n  =  4.82  n',  we  will  sub¬ 
stitute  ^  ^  for  n'.  Thus, 

1.17  =  (4.  82)  (1.17)  =  5.639 
n'  n  n 

Therefore,  5.639  compares  with  .660  when  the  sampling  fractions 
are  equal.  The  variance,  5.639,  might  be  described  as  the 
relative  variance  among  primary  branches  expressed  on  the  basis 
of  one  terminal  branch. 

The  conversion  factor,  4.82,  also  applies  to  the  other 
estimators.  Thus,  all  of  the  unit  variances  for  primary  branches 
must  be  multiplied  by  4.82  to  convert  them  to  the  equivalent  of 
one  terminal  branch.  This  leads  to  Table  2.3,  which  reflects 
differences  in  sampling  efficiency  under  the  condition  that  the 
sampling  fraction  is  the  same  for  primary  and  terminal  branches 
and  for  all  four  plans. 

The  variances  in  Table  2.3  are  also  meaningful  in  terms  of 
sampling  fractions  that  would  be  required  when  all  four  esti¬ 
mators  have  the  same  variance.  Such  sampling  fractions  would  be 
proportional  to  the  variances  in  Table  2.3,  assuming  the  fpc’s 
are  negligible.  As  an  example,  using  primary  branches  as 
sampling  units,  the  variance  of  y2  will  be  the  same  as  the 
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variance  of  y^  when  the  sampling  fraction  of  Plan  2  is  21  per- 
1  191 

cent ,  5-  -g~3§  =  .21,  of  the  sampling  fraction  for  Plan  1.  As 

another  example,  for  the  sampling  variance  of  y^  using  primary 
branches  to  be  the  same  as  the  sampling  variance  of  y^  using 
terminal  branches,  the  sampling  fraction  would  need  to  be  2.8, 


0.897  ..  , 

0 — times  larger. 

Exercise  2.1  (a)  Find  the  relative  variance  of  y^  for  a 

random  sample  of  five  terminal  branches .  Plan  4  is  sampling 
with  pps  and  replacement.  Ans .  0.064. 

(b)  Assume  simple  random  sampling  of  primary  branches  and 
find  the  number  of  primary  branches  so  that  the  relative  vari¬ 
ance  of  (y1)  =  0.064.  The  answer 3  ignoring  the  fpc3  is  18.3. 

There  were  only  28  primary  branches  in  the  population  so  the 

N  -  n 

fpc  should  be  taken  into  account .  Include  the  fpc3  ^  ,  in  the 
variance  formula  for  y^  and  recompute  the  sample  size  that  is 
needed.  Ans.  11  primary  branches. 

(c)  With  reference  to  (a)  and  (b)3  135y^  and  28y^  are 
estimators  of  the  population  total  number  of  apples.  Will  the 
relative  variances  of  these  two  estimators  of  the  total  be 
equal  when  the  sample  sizes  are  5  terminal  branches  with 

Plan  4  and  11  primary  branches  with  Plan  1? 

(d)  It  was  stated  above  that3  when  the  fpc  is  negligible3 
the  variances  in  Table  2.3  are  proportional  to  the  sampling 
fractions  needed  to  have  the  same  relative  variances  of  the 
estimates  for  all  of  the  alternatives.  The  answers  to  (a) 


and  (b )  were 


135 


and 


18.  3 


2"g—  when  the  fpc  was  ignored.  Verify 
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that  these  sampling  fractions  are  proportional  to  the  corres¬ 
ponding  variances  presented  in  Table  2.3. 

Table  2.3  shows  two  major  differences  in  efficiency: 

(1)  Plan  1  vs  the  other  three  plans  and  (2)  primary  vs  terminal 
branches  as  sampling  units.  Table  2.4  presents  the  relative 
variances  for  Plans  2,  3  and  4  as  a  proportion  of  the  variances 
for  Plan  1.  Notice  that  the  proportions  of  the  variation  among 
primary  branches  which  was  accounted  for  by  variation  in  the 
size  (csa)  of  branches  was  much  higher  than  the  proportions  for 
terminal  branches.  Variation  in  the  size  of  terminal  branches 
was  partially  controlled  by  the  specifications  and  process  for 
determining  a  terminal  branch.  For  primary  branches  the  corre¬ 
lation  between  X  and  Y  was  .91.  It  was  .69  for  terminal  branche 

In  Table  2.5  the  variances  for  terminal  branches  are  ex¬ 
pressed  as  a  proportion  of  the  variances  for  primary  branches. 
Here  we  see  that  the  largest  reduction  in  variance  is  under 
Plan  1.  However,  even  after  variance  associated  with  variation 
in  the  csa  has  been  taken  into  account  in  the  estimator  or  pro¬ 
cess  of  selection  (Plans  2,  3  and  4),  the  sampling  variances 
for  terminal  branches  are  about  one  third  of  the  sampling  vari¬ 
ances  for  primary  branches.  This  is  a  manifestation  of  intra¬ 
class  correlation- - the  general  tendency  for  things  that  are 
close  together  in  time  or  space  to  be  alike.  If  there  was  no 
intra-class  correlation,  the  sampling  variances  for  Plans  2,  3 
and  4  would  have  been  about  the  same  for  primary  and  terminal 
branches.  With  Plan  1  the  difference  in  variance  between  pri¬ 
mary  and  terminal  branches  is  attributable  to  the  difference 
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in  correlation  between  csa  and  number  of  apples  as  well  as 
intra-class  correlation. 

The  interaction  shown  in  Table  2.3  between  the  variances 
for  the  four  plans  and  the  two  kinds  of  sampling  units  seems 
typical.  The  situation  might  be  viewed  in  this  way.  When 
the  sampling  units  are  large  and  auxiliary  information  is  not 
used  in  the  sample  design  or  in  estimation,  the  sampling  vari¬ 
ance  is  large  and  there  is  a  large  potential  for  reducing 
sampling  variance.  An  auxiliary  variable  that  is  effective 
in  reducing  sampling  variance  will  probably  be  relatively 
more  effective  when  the  sampling  units  are  large.  This  was 
displayed  in  Table  2.4.  Or,  when  an  effective  auxiliary  vari¬ 
able  is  used,  the  relative  difference  in  sampling  variance 
between  large  and  small  sampling  units  will  probably  be  less 
as  displayed  in  Table  2.5. 

The  same  phenomenon  has  been  observed  in  various  other 
situations.  In  area  sampling,  for  example,  if  geographic 
stratification  is  effective,  it  will  tend  to  be  relatively 
more  effective  when  the  area  sampling  units  are  large  than 
when  they  are  small.  This  is  not  a  justification  for  large 
sampling  units.  The  implication  is  that  matters  of  sample 
design  and  estimation  are  more  critical  when  the  sampling 
units  are  large  and  vary  widely  in  size. 

There  is  a  limit  to  the  reduction  in  variance  that  can 
be  achieved  through  sample  design  and  estimation  techniques. 
That  is,  assuming  a  fixed  sampling  fraction,  one  might  imagine 


a  practical  minimum  variance  as  a  goal  to  be  achieved  by  design. 
There  might  be  a  number  of  alternatives  which  will  approach  that 
goal.  Table  2.3  shows  three  alternatives  with  relative  vari¬ 
ances  between  0.3  and  0.4. 

Exercise  2.2  Table  2.63  which  will  be  discussed  later 3 
shows  the  number  of  apples  on  each  of  the  six  trees  in  the 

EfYj-Y)2 

column  headed  Y^.  The  variances  among  trees3  - -  ,  is 

•  ish 

464,295,  where  is  the  number  of  apples  on  the  i  tree  and 
N  is  the  number  of  trees.  Verify  that  the  relative  variance 
of  Y.  is  0.344.  This  is  the  relative  variance  of  y1  when  a 
tree  is  the  sampling  unit  and  the  size  of  the  sample  is  one 
tree.  Convert  this  variance3  0.344,  to  the  equivalent  of  one 
terminal  branch.  Ans .  1.1A.  Compare  the  answer  with  the 
variances  in  Table  2.2  for  Elan  1. 

Exercise  2.2  Assume  that  a  simple  random  sample  of 
terminal  branches  on  the  six  trees  is  to  be  selected  and  that 
Ny  is  the  estimator  of  the  total  number  of  apples  on  the  six 
trees.  Ignoring  the  fpc3  how  many  terminal  branches  need  to 
be  selected  so  the  variance  of  Ny  is  equal  to  the  variance  of 
an  estimate  based  on  a  random  selection  of  one  tree  and  a  count 
of  all  apples  on  the  tree ?  Assume  that  6y  is  the  estimator  for 
the  sample  of  one  tree  where  y  is  the  number  of  apples  on  the 
sample  tree.  Refer  to  exercise  2.2  for  the  variance  among  trees 
and  to  Table  2.2  for  the  variance  among  terminal  branches. 

Ans.  The  variance  of  an  estimate  from  a  sample  of  2  t erminal 
branches  is  equal  to  the  variance  of  an  estimate  from  a  sample 

I 
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of  one  tree.  There  were  22.5  terminal  branches  per  tree  so 
2  terminal  branches  is  less  than  one-tenth  of  one  tree.  This 
result  is  typical  of  the  low  sampling  efficiency  of  a  large 
sampling  unit.  Moreover3  it  is  very  difficult  to  make  an 
accurate  count  of  all  apples  on  a  tree. 

2.3  STRATIFICATION  BY  TREES 

Table  2.6  presents  variances,  covariances,  and  other 
information  for  each  of  the  six  trees.  These  data  pertain  to 
terminal  branches.  They  will  be  used  to  determine  the  vari¬ 
ances  for  five  different  estimators  based  on  stratified  random 
sampling  with  trees  as  strata  and  a  constant  sampling  fraction. 
Stratified  sampling  with  pps  within  trees  will  also  be  con¬ 
sidered  which  gives  a  total  of  six  alternatives.  For  these 
six  alternatives,  designated  as  plans  6  through  11,  we  want 
to  find  sampling  variances  that  are  comparable  with  the  vari¬ 
ances  presented  in  Table  2.2  for  nonstrat if ied  sampling  of 
terminal  branches. 

It  is  advantageous  to  become  sufficiently  familiar  with 
sampling  theory  to  avoid  searching  textbooks  for  a  formula 
and  checking  it  to  be  sure  it  is  applicable.  A  formula  as 
found  in  a  textbook  might  be  appropriate  but  need  adaptation. 

By  recalling  a  few  things  from  the  theory  of  random  variables, 
correct  variance  formulas  can  be  readily  derived  for  finding 
the  sampling  variances  for  the  sampling  and  estimation  plans 
that  follow. 
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For  comparative  purposes,  relative  variances  are  recorded 
in  Table  2.7  for  four  plans  that  have  been  discussed  and  for 
six  additional  plans  that  will  be  discussed  in  the  next  section. 
2.3.1  PLAN  6 --MEAN  ESTIMATOR 

In  Plan  6  the  sample  is  allocated  to  trees  (strata)  in 
proportion  to  the  number  of  terminal  branches  on  the  trees. 

You  may  notice  that  Plan  6  is  the  same  as  Plan  5  except  that 
the  strata  are  trees  instead  of  size- of-branch  classes.  The 
estimator  of  the  population  mean,  Y,  is 


^6 


N 


1  - 


FT  *1  + 


(2.2) 


where  y^ 


is  the  sample  average  for  stratum  h 
(i.e.,  the  average  number  of  apples  per 
terminal  branch  on  tree  h) , 


h  is  the  index  for  strata  (trees), 
i  is  the  index  for  sampling  units  within  stratum  h 
(branches  on  tree), 

N^  is  the  total  number  of  sampling  units  (terminal 
branches)  in  stratum  h, 

N  =  ENh  is  the  total  number  of  sampling  units  in 
the  population,  and 

n^  is  the  number  of  sampling  units  in  the  sample  from 
stratum  h  (number  of  branches  in  the  sample  from 
tree  h) . 
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Exercise  2.4  Since  the  sample  is  stratified  and  allo¬ 


cated  to  strata  in  proportion  to  N^,  the  estimator  is  a  simple 
average  of  all  values  of  y,^  in  the  sample.  Show  that  this  is 
true . 

The  estimator,  y^ ,  was  written  as  shown  in  Eq.  (2.2) 

because,  to  find  its  variance,  we  need  to  consider  it  as  a 

Nh 

function  of  the  stratum  means.  The  weights  ,  are  constant. 

A 

Therefore,  the  variance  of  y^  depends  on  the  variance  of  the 

stratum  means.  The  sample  from  one  stratum  is  independent  of 

the  sample  from  another  stratum.  Therefore,  the  stratum  means , 

Nh 

y^,  are  independent  random  variables,  and  the  terms,  y^,  in 
y^  are  independent  random  variables.  We  know  from  the  theory 
of  random  variables  that  the  variance  of  the  sum  of  independent 
random  variables  is  the  sum  of  the  variances  of  the  random  vari¬ 
ables.  This  gives  the  basis  for  writing  the  variance  of  y^  as 
follows : 

V(y6)  =  x[V(Jl  yh)]  =  vcjlyp  ♦...+  V(JJiyL) 

h 

We  also  know  that  the  variance  of  a  constant  times  a  variable 
equals  the  square  of  the  constant  times  the  variance  of  the 
variable.  Hence, 

V(y6)  =  E  [  (ip)  2V  (y,  )  ]  =(jr)2V(y  )  +  ...  +  (JjVv(yL)  (2.3) 

h 


Next,  we  need  an  expression  for  the  variance  of  y^.  Since  the 
sample  within  each  stratum  is  a  simple  random  sample,  the  vari¬ 
ance  of  y^ 


is  as  follows: 


N.-n,  SYh 

vcyj  -  C^V1)-  Yh 


n. 


(2.4) 


where 


Yh 


phrV 

1 


The  subscript  Y  in  is  included  to  show  that  the  variance 

refers  to  the  variable  Y.  Later  we  will  need  to  take  the 

2 

variance  of  X  into  account  and  will  use  to  represent  the 

variance  of  X  within  stratum  h  and  S^y^  to  represent  the  co- 

variance  of  X  and  Y  within  stratum  h. 

For  simplicity  and  convenience  assume  that  the  sampling 

nh  Nh'nh 

fractions,  f^  =  ,  are  small  so  the  fpc’s,  — |q —  ,  may  be 

n  n 

ignored.  Thus,  dropping  the  fpc  and  substituting  the  variance 

of  y^  in  Eq.  (2.3)  gives: 


N,  ,  sl.  N,  ,  sj,  N.  ,  S2 

V(yJ  =  -H  =  C/)2  -li  +  ...♦  (rji)2  -H 

w6  VN  nn  n,  J  nT 


(2.5) 


Since  the  sampling  specifications  called  for  a  constant  sampling 


n 


fraction,  is  constant  from  stratum  to  stratum  which  means 


that 


n. 


N, 


_  n 
Nl  -  N 


where 


En^  =  n  and  EN^  =  N 


n. 

Substituting  ^  for  —  in  Eq.  (2.5)  and  simplifying  the  ex- 


N, 


pression  we  obtain 

A  -I  Nl  0  “\  ^  1  O 

'  k  2  IT  4  '  ^  S' 


n  lN  °Y1 


NI  2 

+  +  _ "  Q  ^  "1 

...  N  ^YLJ 


(2.6) 
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Exercise  2.5  Perform  the  algebra  that  is  necessary  to 
go  from  Eq.  (2.5)  to  Eq.  (2.6). 

For  Plan  6,  let 


S6  ■  EWh4  ■  W1SY1 


WLSYL 


(2.7) 


N, 

where  \  "  jj- 


2 

Since  Sy^  will  be  replaced  by  corresponding  variances  that  are 

2  2 

involved  later  in  Plans  7  and  8,  let  =  Sy^  so  the  notation 

will  reflect  the  number  of  the  plan  or  estimator.  Then  Eq.  (2.7) 
becomes 


S 


2 

6 


EWhS6h 


(2.8) 


and  Eq.  (2.6)  simplifies  to  the  following  form 

y(y6)  ■  &s6  (2-9) 

2 

where  is  a  weighted  average  of  the  within  stratum  variances. 

2 

The  values  of  are  recorded  in  Table  2.6  in  the  column 
2  2 

headed  and  the  value  of  is  1367  which  is  recorded  in 
the  line  labeled  '’Separate.”  The  reason  for  calling  this  line 
"Separate”  will  be  explained  later. 

For  purposes  of  comparing  variances  for  alternative  plans 
the  choice  of  a  sample  size  is  arbitrary.  Previously,  the 
sampling  variances  for  alternative  plans  were  compared  assuming 
n  =  1.  Even  though  it  is  impossible  to  select  a  stratified 
random  sample  of  only  one  unit,  it  is  possible  to  let  n  =  1 
in  Eq.  2.9  and  regard  the  variance  of  y^  as  the  sampling 
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variance  for  a  hypothetical  sample  of  one  unit.  As  with  simple 
random  sampling,  a  stratified  random  sample  of  n  units  would 
have  a  sampling  variance  equal  to  jjj-  times  the  sampling  vari¬ 
ance  for  a  hypothetical  stratified  random  sample  of  one  unit-- 
provided  n  is  large  enough  so  the  n^,  which  must  be  integers, 
are  approximately  in  proportion  to  N^.  Remember,  these  numer¬ 
ical  examples  are  being  worked  as  though  the  sampling  fraction, 
f^,  is  constant  and  small. 

Exercise  2.6  Calculate  the  variance  of  y^  assuming 

2 

n  =  1.  In  other  words 3  find  the  value  of  .  Also 3  calculate 
the  relative  variance  of  y^  when  n  =  1.  Your  answer  should 
agree  with  the  relative  variance  of  y  ^  which  is  recorded  in 
To.hle  2.7. 

Exercise  2.7  Since  y^  is  an  estimate  of  Y,  Ny^  is  an 
estimate  of  the  population  total.  Find  the  standard  error  of 
Ny^  for  n  =  1.  Ans.  4991. 


2.3.2  PLAN  7- -RATIO  ESTIMATORS  BY  STRATA 


Yi 


Plan  7  is  the  same  as  Plan  6  except  that  (X,  — ) , 


instead  of  y^,  is  used  in  Eq.  2.2  as  an  estimator  of  the 
stratum  mean,  Y^.  Thus, 


N, 


y,  N,  y,  nt  yT 

Y 7  ~  1  N~^h  =  n-  T~ )  +...+  CXL  t~)  (2.10) 


N. 


The  derivation  of  the  relative  variance  of  y^  follows  the 
derivation  in  Plan  6.  Simply  replace  the  variance  of  y^  in 
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-  yh 

Eq.  (2.3)  with  the  variance  of  (X^  7—). 

xh 


The  variance  of 


(X,  — ) ,  ignoring  the  fpc,  is 


-  yh  12 

V(X,  — )  =  (— )S* 

K  h  -  n,  7h 

xh  h 


where 


2  2  2  2 

S7,  =  +  Rr  S*  -  2R,  Svv, 

7h  Yh  h  Xh  h  XYh 


2  2 

Notice  that  Sy^  is  the  same  as  S£  in  Table  1.1  except  that 
2 

Sy^  is  a  variance  within  stratum  h  rather  than  a  variance  over 

2  2 

the  whole  population.  Substituting  Sy^  for  in  Eqs .  (2.5), 
(2.6),  and  (2.7)  leads  to  the  following  results: 


vCy7)  -  (i)s2? 


(2.11) 


where  S?  =  EW,S~, 

7  h  7h 

2  2 

The  values  of  Sy^  and  Sy  are  presented  in  Table  2.6. 

Exercise  2.8  The  estimator  yy,  Eq.  (2.10  )t  was  expressed 
in  a  form  to  show  its  similarity  to  y^.  Is  there  a  modifica¬ 
tion  of  Eq.  (2.10)  that  would  he  better  for  computing  the  value 

A  A 

of  yy  from  sample  data?  How  would  you  compute  the  value  of  yy 
from  a  sample? 

Exercise  2.9  From  the  data  presented  in  Table  2.6 ,  find 
the  relative  variance  of  y*j  for  n  =  1.  The  answer ,  0.279,  is 
in  Table  2.7.  How  would  you  explain  why  the  sampling  variance 
for  y ~j  is  less  than  the  sampling  variance  for  y^  ? 
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The  sampling  variance  for  Plan  2  (no  stratification  and 
the  ratio  estimator)  was  0.382  compared  with  0.279  for  Plan  7 
(stratification  and  separate  ratio  estimators  by  strata).  The 
geometrical  interpretation  of  the  sum  of  squares  for  Plan  7 
compared  with  Plan  6  is  analogous  to  Plan  2  compared  with 
Plan  1.  Horizontal  lines  for  Plan  6  (one  for  each  tree)  are 
replaced  by  lines  through  the  origin  and  the  stratum  means  of 
X  and  Y.  With  the  ratio  estimator,  y^,  the  effect  of  strati¬ 
fication  depends  on  how  much  the  ratio  lines  differ  among 
strata.  More  will  be  said  later  about  stratification  and 
ratio  estimators. 

Exercise  2.10  Notice  with  reference  to  Eq.  (2.10)  that 
Nh^h  is  the  population  total  of  X  for  stratum  h.  Let 
X^  =  N^X^  and  substitute  X^  in  Eq.  (2.10)  which  gives 


In  .  1  [X 
-  N  LX1 
n 


With  y y  in  this  forms  write  a  formula  for  the  variance  of  y^. 

2.3.3  PLAN  8 --REGRESSION  ESTIMATORS  BY  STRATA 

Plan  8  is  like  Plans  6  and  7  except  that  the  regression 
estimator  (see  Plan  3,  Chapter  I)  is  used  stratum  by  stratum. 
Thus,  instead  of  Eq.  (2.2)  or  (2.10)  we  have 

-  E<?r  Cyh  +  V*h  -  V]} 

N1  -  -  nt 

~  W~  1  +  bi(xrxi)l  +*..+  [yL+bLCXL-xL)]  (2.12) 
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In  the  derivation  of  the  variance  of  yg,  the  variance  of 
^h+^h^h~xh^  replaces  the  variance  of  y^  in  Eq .  (2.3). 

A 

This  leads  to  an  equation  for  the  variance  of  yg  which  is 
similar  to  the  variances  of  y^  and  y^.  Thus, 

VCV  =  S8 

where  S8  *  EWhS8h 

and  S2h  =  S2h(l-r2) 


where  r,  is  the  correlation  between  X  and  Y  within  stratum  h. 
h 


Exercise  2.11 

Find 

the 

re lative 

variance 

Of  y8 

for  n  =  1 

Compare  your  result 

with 

the 

relative 

variance 

for  y 

g  that  is 

recorded  in  Table  2. 

7 . 

2.3.4  DISCUSSION  OF  PLANS  6,  7,  and  8 

Compare  the  estimators,  y^,  y^,  and  yg,  and  their  vari- 

A  A  A 

ances  with  y^,  y  ^  and  y^,  and  their  variances,  Table  2.7.  In 
essence  each  stratum  in  Plans  6,  7,  and  8  is  treated  as  a 
separate  population  and  the  estimators  and  their  variances  with¬ 
in  each  of  the  strata  are  combined  using  appropriate  weights. 
Geometric  interpretations  of  the  sampling  variances  with  refer¬ 
ence  to  sums  of  squares  is  analogous  to  the  interpretations  given 
in  Chapter  I  for  Plans  1,  2,  and  3.  There  is  one  line  for  each 
stratum  and  each  of  the  estimators,  y^,  y^,  and  yg. 
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Figure  2.3  presents  a  dot  chart  for  each  of  the  six  trees. 

A 

For  each  tree,  the  solid  line  is  the  ratio  line  involved  in 
and  the  broken  line  is  the  ratio  line  for  y no  stratification. 

A  A 

As  recorded  in  Table  2.7,  the  relative  variances  of  y  ^  and  y^ 
are  0.279  and  0.382  respectively.  This  indicates  the  degree  to 
which  the  6  ratio  lines  fit  the  data  better  than  the  single 

A 

line.  Figures  analogous  to  Fig.  2.3  could  be  prepared  for  y^ 
compared  with  y^,  for  yg  compared  with  y^,  for  yg  compared  with 

y  7 ,  o  t  c . 

Separate  stratum  estimators  like  y  ^  and  y-g  are  seldom  used 
in  practice.  However,  Plans  7  and  8  were  included  for  compara¬ 
tive  purposes  and  further  understanding  of  possible  alternatives. 
There  will  be  additional  discussion  of  these  plans  after  Plans 
9,  10,  and  11  have  been  presented. 

2.3.5  PLAN  9  -  -  COMBINED  RATIO  ESTIMATOR 

Instead  of  making  a  ratio  estimate  for  each  stratum  and 
combining  the  separate  stratum  estimates,  the  data  from  the 
strata  are  combined  before  computing  a  ratio.  Likewise,  in 
Plan  10,  results  for  individual  strata  are  combined  and  used  to 
determine  a  "combined  regression  estimator."  This  explains  the 
two  titles  "Separate"  and  "Combined"  in  Table  2.6.  The  "Separate" 
line  contains  averages  of  within  stratum  variances  for  Plans  7, 

8,  and  11  which  use  separate  stratum  estimators.  The  entries 
in  the  "Combined"  line  pertain  to  the  combined  stratum  estima¬ 
tors  in  Plans  9  and  10.  The  distinction  between  separate  and 
combined  is  not  applicable  to  the  mean  estimator,  Plan  6. 
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2 

S6  is  shown  in  both  lines  of  the  table. 


The  "combined  ratio  estimator"  is 


(2.14) 


where 


Nh 

*s  =  ZN“  ^h 


x 


s 


The  letter  "s"  in  y  and  x  is  used  to  indicate  that  y  and  x 

s  s  J  s  s 

are  means  that  pertain  to  a  stratified  random  sample. 

To  find  the  variance  of  y^ ,  it  is  convenient  to  remember 
that  the  large  sample  approximation  of  the  relative  variance 
(RelVar)  of  the  ratio  of  any  two  random  variables  u  and  v  is 

RelVar(^)  =  RelVar(u)  +  RelVar(v)  -  2RelCov(u,v) 

Therefore,  since  y  and  x  are  random  variables  we  have 

J  s  s 

ys  .  _  .  . 

RelVar  (-3— )  =RelVar (y  ) +RelVar  (xs ) - 2RelCov (yg , x  )  (2.15) 

x 

s 

Exercise  2.12  Verify  that  the  relative  variance  of  y g 

ys 

is  equal  to  the  relative  variance  of  the  ratio ,  —  . 

*s 

With  reference  to  Eq.  2.15,  notice  that  yg  is  the  same  as 

/\  /v 

y^.  We  found  for  Plan  6,  Exercise  2.6,  that  the  RelVar  of  y^, 
and  therefore  of  y  ,  was  0.512  for  n  =  1.  The  RelVar  of  xg  is 
determined  in  the  same  way.  According  to  Table  2.6,  the  average 
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-  2 

within  stratum  variance  of  X  is  0.2566  and  X  =  2.0240.  Thus, 
RelVar(xs)  for  n  =  1  is  =  0.127. 

Exercise  2.13  Find  the  average  within  stratum  covariance 
of  X  and  Y  in  Table  2.6.  Then  find  the  value  of  RelCov (ys , xg) 
for  n  =  1.  Ans .  0.166. 

ys 

According  to  Eq.  (2.15)  the  RelVar  of  —  is 

*s 

0.512  +  0.127  -  2(0.166)  =  0.307 

Therefore,  the  RelVar  of  yg  is  0.307.  This  answer  is  recorded 
in  Table  2.7. 

Exercise  2.14  Start  with  Eq.  2.15  and  show  that  the 
variance  of  yg  is  given  by 

V(yg)  =  V(ys)  +  R2V(xs)  -  2R[Cov(ys,xs)] 

Y  .  *  -  2 

where  R  =  —  .  Suggestion:  Notice  that  V(yQ)=Y  [RelVar  (yQ) ] , 

X  y  y 

-  2 

then  multiply  the  right  hand  side  of  Eq.  (2.15)  by  Y  .  For  n  =  1 
the  value  of  V(yg),  V(xg),  Cov(ys,xs)  and  R  are  given  in  the 
"Combined"  line  of  Table  2.6.  Using  these  values  compute  the 

A 

value  of  V(yg) .  The  answer  is  817,  which  is  also  in  the  "Com¬ 
bined"  line. 

Exercise  2.15  Beginning  with  the  variance  of  yg  as  expressed 
algebraically  in  Exercise  2.14 ,  show  that 

v(y9)  -  i 
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where 


S9  ^WhS9h 


and 


9h 


=  S 


Yh 


R  4 


-  2RS 


XYh 


Suggestion :  Sinoe 

Eq.  (2.9)  to  obtain 
Formulas  for  V(xs) 

S9h  wzth  S7h ’ 


yg  and  y^  are  the  same  you  may  refer  to 
the  appropriate  formula  for  V  (yc ) . 
and  Cov(xs,ys)  are  analogous .  Compare 


Exercise  2.16  Continuing  from  the  formula  for  V(yg) 
which  is  given  in  Exercise  2.153  show  that 


v(y9)  =  i  rfMh 


^Yhi-Rxhi 

1 

N,“  - 1 


2 

-] 


The  formula  for  the  variance  of  yg ,  which  is  given  in 
Exercise  2.16,  shows  that  the  deviations  which  are  squared 
are  deviations  from  a  line  through  the  origin  and  (X,Y)  where 
X  and  Y  are  the  overall  means  of  X  and  Y.  This  line  for  the 
combined  ratio  estimator,  yg ,  is  the  same  as  the  line  pertaining 
to  y2»  the  ratio  estimator  without  stratification.  Thus,  if 

A  /V 

yg  has  a  lower  variance  than  y ^  it  is  attributable  to  the  effect 
of  stratification  which  assures  proportional  representation  in 
the  sample  by  strata.  That  is,  there  is  proportional  repre¬ 
sentation  by  strata  of  the  deviations  of  Y^  from  the  combined 
ratio  line.  In  Plan  7  there  was  proportional  representation 
and  separate  ratio  lines  by  strata.  RelVar  of  yg  was  0.307 
compared  to  0.382  for  y2  and  0.279  for  y^  (see  Table  2.7). 
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2.3.6  PLAN  10- -COMBINED  REGRESSION  ESTIMATOR 

As  in  the  case  of  the  combined  ratio  estimator,  data 
for  strata  may  be  combined  and  a  single  (or  combined)  regres¬ 
sion  used  instead  of  separate  regressions.  The  estimator,  y^ , 
for  the  combined  regression  looks  like  y^  but  it  is  an  average 
within  stratum  regression  that  is  determined  from  combined 
within  stratum  variances  and  covariances.  Since  the  sampling 

fraction  is  constant,  the  appropriate  weights  for  combining  the 

N1  NL 

within  stratum  variances  and  covariances  are  , .  .  .  ,  which 
are  the  same  weights  used  previously  for  combining  variances. 
The  combined  within  stratum  variances  of  Y  and  X,  1367  and 
.2566,  and  the  combined  within  stratum  covariance,  12.23,  are 
shown  in  the  "Combined”  line  of  Table  2.6.  These  numbers  are 
needed  for  computing  the  "Combined"  regression  coefficient, 
the  "Combined"  correlation  coefficient,  and  the  variance  of 

A 

y^g  for  n  =  1,  which  are  also  shown  in  the  "Combined"  line  of 
Table  2.6.  The  corresponding  numbers  for  sampling  without 
stratification  are  shown  in  the  last  line  of  Table  2.6. 
Algebraically 


y 


10 


=  y  +  bc(X-xc) 


(2.16) 


where 


b  = 
s 


XY 


?X 


SXY  "  £WhSXYh 

SX  ■  f  hSXh 
N, 

wh  - 


and 
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Notice  that  lower  case  letters,  x  and  y,  are  used  in  the 


definition  of  bg  to  indicate  that  it  is  computed  from  sample 

values.  In  Table  2.6,  the  value  of  B  is  shown  which  is  the 

s 

population  value  that  bs  is  an  estimate  of.  The  bar  in  the 


of  within  stratum  covariances  and  variances.  (Previously, 


2  2 

we  had  used  Sy,  S^.,  and  S^.y  to  represent  the  overall  variances 


and  covariances  without  stratification.)  The  subscript  MsM 
is  used  as  a  code  indicating  that  stratified  random  sampling 
and  combined-stratum  estimation  are  involved.  To  recapitulate 
b^  is  a  least  squares  estimate  of  the  regression  coefficient 
within  stratum  h,  bg  is  an  estimate  of  the  combined  regression 
coefficient  in  the  combined  regression  estimator,  and  b  is  the 
least-squares  regression  coefficient  computed  from  a  simple 
random  sample  without  stratification. 

The  variance  of  y1Q  is 


(2.17) 


where 


S 


XY 


and 


r 


s 


The  variance  of  y1Q  involves  squares  of  the  deviations  of 


from  a  line  with  a  slope  equal  to  Bs  that  passes  through  (X,Y) 
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Remember  the  assumptions  and  that  the  variance  formula  is  a 

large  sample  approximation.  Further  discussion  of  the  basis 

for  the  formula  for  the  variance  of  y^Q  will  be  omitted.  For 

4  / 

more  detail  the  reader  is  referred  to  Cochran.— 


Exercise  2.17  Verify  that  the  regression  coefficient  for 
the  combined  within  stratum  regression  is  47 .7  and  that  the 
combined  within  stratum  correlation  coefficient  is  0.653.  Then 
verify  that  the  relative  variance  of  y-^Q  is  0.294  for  n  =  1. 

2.3.7  PLAN  11- -SAMPLING  WITH  PPS  WITHIN  STRATA 


As  in  Plan  4,  sampling  with  replacement  is  assumed  for 

simplicity.  That  is,  in  stratum  h  the  probability  of  the  / 

sampling  unit  being  selected  on  any  given  draw  is  proportional 

to  X...  The  estimator  is 
hi 

yn  =  £  tfr5  yiih  (2.18] 


where 


-  A, /hi 

^llh  ^n,  . 

h  l  hi 


Notice  that  y-Q^  is  like  y^,  the  difference  being  that  y^^  is 
an  estimate  of  the  stratum  mean  whereas  y^  is  an  estimate 
of  the  population  mean  Y.  Also,  notice  that  the  estimator, 
y^,  can  be  obtained  by  substituting  y-Qh  f°r  in  2.2. 


4/  Cochran,  W.  G.,  Sampling  Techniques,  Second  Edition, 
Chapter  7.  John  Wiley  §  Sons,  Inc.,  1963. 
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It  follows  that  the  variance  of  y^  can  be  obtained  by  substi¬ 
tuting  the  variance  of  y-Qh  f°r  the  variance  of  y^  in  Eq.  (2.3). 

A  A 

Owing  to  the  similarity  of  y-Q^  and  y^  the  formula  for 
the  variance  of  y^  is  applicable.  Simply  add  a  subscript  h 
in  the  formula  for  the  variance  of  y^ ,  which  gives 


where 


and 


v<Yllh>  ■  CHTH-4UPhiCpj=i 

h  1  hi 


V 


x,. 

hi 


hi  X, 


Y,  =  ZY,  . 
h  .  hi 


X,  =  EX,  . 
h  .  hi 
1 


Substitution  of  V(y-Q^)  i-n  Eq.  (2.3)  gives 

vCyn)  =  EC^)2[(i-)(4)EPhi(^i 

h  h  l  hi 


Yh)2] 


(2.19) 


As  in  the  derivation  of  2.6,  assume  that  n^  is  propor¬ 
tional  to  N^,  which  means  that  n^  =  Substituting 

(^•)N^  for  n^  in  Eq.  (2.19)  leads  to  the  following  which  ex¬ 
presses  the  variance  of  y^  in  a  form  like  that  used  for  the 
other  estimators: 


*  ^Sll 


(2.20) 


where 


sii  ■  fhsiih 


and 


llh 


i  ZP  (-2ii 
7T  .  hi1?,  • 
N,  l  hi 
h 


V 
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2 

Like  the  variances  for  the  other  estimators,  S-^  is  a  weighted 

average  of  the  appropriate  within  stratum  variances.  The 

2 

within  stratum  variances,  S^^*  are  presented  in  Table  2.6  and 
2 

S-^,  the  variance  of  for  n  =  1,  is  recorded  in  the  ’’separate" 

1  ine . 


Exercis e  2.18 
of  S^lh  for  h  =  1. 
in  Table  2.6.  Is 
above  for  finding 
Exercise  2.19 
presented  in  Table 
check  your  answer 
Table  2.  7. 


Using  the  data  in  Table  2.13  find  the  value 

Check  your  answer  with  the  value  recorded 

there  a  better  expression  than  the  one  given 

2 

the  values  of  S-^-^? 

From  the  data  by  individual  trees  that  are 

2.63  find  the  RelVar  of  y^  for  n  =  1  and 

2 

with  the  value  of  that  is  recorded  in 


A  geometrical  interpretation  of  th 
matter  of  making  an  interpretation  for 
the  average  situation  over  all  strata, 
reference  is  made  to  the  discussion  of 

2.3.8  SUMMARY  AND  DISCUSSION 


e  variance  of  y^ 
each  stratum  and 
For  this  purpose 
y^  in  Chapter  I. 


j 


is  a 
udging 


Sampling  variances  for  10  out  of  11  plans  are  presented 
in  Table  2.7  for  terminal  branches  as  sampling  units.  Plan  5 
was  not  applied  to  terminal  branches.  All  of  the  plans  have 
an  important  practical  shortcoming.  It  is  necessary  to  define, 
label,  and  list  all  terminal  branches  on  a  tree  before  it  is 
sampled.  Some  ways  of  avoiding  this  will  be  discussed  in 
Chapter  III,  However,  Chapter  II  was  intended  as  an  exercise 
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in  the  use  of  theory  to  find  the  variances  for  alternative 
sampling  and  estimation  plans  and  as  a  study  of  the  differ¬ 
ences  in  the  variances  for  several  alternatives. 

As  you  gain  experience  through  evaluations  of  sampling 
plans  you  will  become  increasingly  aware  of  prevailing  patterns 
of  variation.  You  will  observe  manifestations  of  the  general 
tendency  for  things  to  be  stratified  in  space  or  time,  or  the 
tendency  for  things  that  are  close  together  in  space  or  time 
to  be  alike.  There  are  exceptions.  For  example,  in  a  field 
where  the  plant  population  is  very  dense  there  might  be  a 
negative  intra-plot  correlation  among  plants  within  very  small 
plots  owing  to  competition  between  adjacent  plants. 

From  the  results  in  Table  2 . 7  we  find  that  the  two  plans 
with  the  largest  variance  are  Plans  1  and  6.  Neither  plan 
makes  use  of  size  of  branch  as  an  auxiliary  variable.  The 
reductions  in  variance  from  use  of  csa  as  an  auxiliary  variable 
are  substantial.  This  strongly  suggests  exploration  of  practi¬ 
cal  ways  of  using  csa  as  a  measure  of  branch  size  unless  it  is 
possible  and  feasible,  when  determining  terminal  branches,  to 
restrict  the  sizes  within  narrow  limits. 

The  variances  for  "separate"  ratio  and  regression  esti¬ 
mators  are  moderately  less  than  the  corresponding  variances 
for  the  "combined"  ratio  and  regression  estimators.  A  small 
difference  in  favor  of  "separate"  estimators  is  indicated  by 
general  experience  and  mathematical  considerations.  However, 
"combined"  ratio  or  regression  estimators  are  generally  used 
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in  practice  because:  (1)  they  are  more  convenient;  (2)  in  some 
situations,  bias  in  the  '’separate"  estimators  is  appreciable 
relative  to  the  standard  error;  and  (3)  the  variance  formulas, 
which  are  large  sample  approximations,  are  better  approxima¬ 
tions  for  the  "combined"  estimators.  Separate  ratio  or 
regression  estimators  might  be  preferable  to  combined  esti¬ 
mators  when  the  number  of  strata  is  very  small  and  the  ratios 
or  regressions  differ  widely  among  the  strata.— ^ 

With  stratified  random  sampling,  sampling  variance  is  a 
function  of  variation  within  strata.  It  is  generally  better 
to  judge  the  impact  of  stratification  by  considering  within 
stratum  variation  than  by  the  differences  among  strata.  Making 
a  choice  between  two  alternative  methods  of  stratification 
solely  on  the  basis  of  differences  among  strata  could  in  some 
cases  be  misleading.  For  example,  the  variance  among  the  means 
of  four  strata  could  be  much  larger  than  the  variance  among  the 
means  of  30  strata.  That  does  not  necessarily  mean  that  the 
sampling  variance  for  the  four  strata  will  be  the  least.  Also, 
the  effect  of  stratification  and  of  optimum  allocation  among 
strata  are  not  independent  of  the  method  of  estimation. 

In  the  preceding  discussion,  stratification  was  considered 
as  a  matter  of  reducing  sampling  variance.  In  the  design  of 
a  sample  for  a  survey,  attention  needs  to  be  given  to  the 
domains  (subpopulations)  for  which  estimates  are  important. 

5_/  See  Cochran,  Sampling  Techniques,  for  a  discussion  of  the 
properties  of  the  separate  and  combined  estimators. 
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lhis  might  be  a  primary  determiner  of  the  stratification  and 
allocation  of  the  sample.  The  sampling  variances  of  domain 
estimates  depend,  among  other  things,  on  how  close  the  bound¬ 
aries  of  strata  for  sampling  purposes  correspond  to  the  domains 
for  which  estimates  are  required. 

2.4  FURTHER  COMPARISON  OF  SAMPLING  WITH  PPS  TO  STRATIFIED 

SAMPLING  WITH  OPTIMUM  ALLOCATION 

In  Chapter  1,  sampling  with  PPS  was  compared  to  stratified 
random  sampling  with  optimum  allocation  to  strata.  The  data 
for  terminal  branches  provide  a  better  set  of  data  for  study 
of  sampling  with  pps  including  the  possibility  of  a  transfor¬ 
mation  of  X  or  Y  to  reduce  sampling  variance.  For  this  purpose, 
five  size-of-branch  strata  based  on  csa  will  be  used.  Table  2.8 
shows  the  definition  of  these  strata  and  presents  key  informa¬ 
tion  about  the  five  strata.  (Reference  is  made  to  Sections  1.2.5, 
1.2.6,  and  1.3.1  on  sampling  with  pps  and  stratified  random 
sampling  with  optimum  allocation.)  For  stratified  random 
sampling  using  the  mean  estimator,  Eq.  2.2,  the  optimum  sampling 
fraction  for  stratum  h  is  given  by 


f'  = 
±h 


hSYh 


)n 


In  previous  comparisons  we  assumed  n 
ested  in  the  values  of 


£h  = 


Yh 


ZNhSYh 


1.  Hence,  we  are  inter- 


(2.21) 
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which  are  presented  in  Table  2.9.  The  values  of  are  the 

sampling  fractions  within  strata  for  a  hypothetical  sample  cf 

X. 

one  branch  and  are  comparable  to  the  probabilities  =  j—  for 

t  h 

selecting  one  branch  with  pps ,  where  X.  is  the  csa  of  the  1 
N  1 

branch  and  X  =  EX^.  Let  equal  the  average  value  of  P^  for 
the  branches  in  stratum  h.  The  values  of  P^  are  presented  in 
Table  2.9  for  comparison  with  values  of  f^. 

Exercise  2.20  Verify  the  values  of  f^  and  P^  in  Table  2.9 
for  one  or  two  of  the  strata.  The  data  presented  in  Table  2.8 
are  sufficient  for  this  purpose. 

The  values  of  f^  and  P^  agree  quite  well,  which  means  the 
probability  of  a  branch  being  in  a  sample  is  roughly  the  same 
for  both  methods,  except  for  variation  of  P^  within  a  stratum. 
The  next  question  is  how  well  do  the  lines  that  are  involved 
fit  the  data.  Turn  to  Figure  2.4.  The  points  appear  to  fit 
horizontal  lines  (which  are  not  shown)  for  stratified  random 
sampling  approximately  as  well  as  the  line  through  the  origin 
and  (X,Y)  for  pps. 

The  variance  for  sampling  with  pps,  Plan  4,  has  already 
been  obtained.  According  to  Table  2.2  it  is  0.319.  For  the 
stratified  sampling  with  optimum  allocation,  which  will  be 
called  Plan  12,  the  estimator  of  Y  is 

y12  =  CR)ENhyh  (2‘22) 
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and  the  variance  of  y^7  is  given  by 


where 


v^i2)  ■  |(si2^ 


c 2  _  rENh  SYh,2 

b12  1  N  J 


(2.23) 


Exercise  2.21  Using  the  data  in  Table  2.8}  find  the 
RelVar  of  y ^  fOI>  n  =  1.  Ans.  0.273. 

Exercise  2.22  From  Eqs.  (2.5)  and  (2.21)  derive  alge¬ 
braically  the  variance  of  y ^  which  is  given  by  2.23. 

In  this  example,  the  RelVar  for  stratified  random  sam¬ 
pling  with  optimum  allocation  was  0.273  compared  with  0.319 
for  sampling  with  pps .  The  difference  in  variance  is  attrib¬ 
utable  to  the  difference  in  probabilities  of  selection  and  to 
how  well  the  lines  that  are  implicitly  involved  fit  the  data. 

A  question  posed  earlier  was  whether  some  transformation 
of  X  might  provide  a  better  measure  of  size  for  pps  sampling. 

A  simple  transformation  would  be  Xl  =  +  C,  where  C  is  a 

constant  and  X?  is  the  transformed  variable.  The  least  squares 
regression  line,  see  Figure  2.4,  crosses  the  horizontal  axis 
at  0.46.  Since  the  least  squares  line  fits  the  data  "better" 
than  any  other  straight  line  the  transformation  X^  -  0.46  is 
suggested.  With  this  transformation  the  least  squares  and 
pps  lines  become  the  same.  But,  with  the  transformation 
X^  -  0.46  the  maximum  values  of  do  not  approach  zero  as 
X7  =  X^  -  0.46  approaches  zero  which  indicates  that  the  trans¬ 
formation  is  not  a  good  one.  However,  it  is  informative  to 
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make  the  transformation.  Instead  of  p:  =  ,  we  now  have 

1  A 

xr  -  0.46 

“  e (X. -0 . 46)  *  ^he  average  values  of  P7  within  strata, 

which  are  labeled  P^,  are  presented  in  Table  2.9  for  compari¬ 
son  with  the  values  of  P^  and  the  optimum  sampling  fractions, 

fh’ 

When  X^  is  used  as  the  auxiliary  variable,  the  relative 
variance  is  1.403  compared  to  0.319  when  X^  is  the  auxiliary 
variable.  The  relative  variance  for  Plan  1  was  only  0.660. 
Thus,  the  transformed  variable  Xr  results  in  an  increase  in 
variance  compared  to  simple  random  sampling  with  equal  proba¬ 
bilities  of  selection.  Before  transformation,  the  selection 
probabilities  for  sampling  with  pps  were  .0032  and  .0140  for 
the  smallest  and  largest  branches  and  were  .0012  and  .0171 
after  transformation.  Thus,  the  range  in  the  selection  proba¬ 
bilities  were  greatly  increased  by  the  transformation  from  a 
factor  of  4.3  =  to  a  factor  of  14.2  =  .  The 

transformation  does  not  effect  the  optimum  sampling  fractions, 

££- 

Exercise  2.23  Verify  two  of  the  values  of  P^  that  are 
presented  in  Table  2.9. 

Exercise  2.24  Verify  that  when  the  transformation 

Xi  =  Xi  ‘  °-46  is  made  the  pps  line  and  the  least  squares 
regression  line  become  the  same. 
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With  reference  to  Figure  2.4,  consider  two  lines  through 
the  origin  that  represent  maximum  and  minimum  values  of  Y. 

(See  Figure  1.2  in  Chapter  I  which  was  portrayed  as  a  good 
case  for  pps).  As  an  approximation,  theory  suggests  that  a 
good  measure  of  size  is  one  that  is  proportional  to  the  differ¬ 
ence  between  two  lines  representing  the  maximum  and  minimum 
values  of  Y.  A  look  at  Figure  2.4  with  this  in  mind  suggests 
that  a  transformation  such  as  -  0.46  is  not  a  good 

possibility  for  reducing  variance. 

Exercise  2.25  Refer  to  Table  2.8  and  for  each  stratum 
divide  Sy^  by  X^  and  examine  Sy^  as  a  proportion  of  X^.  Does 
this  indicate  that  a  transformation  of  X  would  be  advisable? 

Ans .  No.  A  transformation  of  X  is  not  indicated  and  one  should 
accept  csa  as  a  measure  of  size  unless  there  is  evidence  to  the 
contrary  from  other  sources . 

Table  2.10  provides  a  comparison  of  variances  for  four 
plans  when  X  and  X'  =  X  -  0.46  are  the  auxiliary  variables. 

Exercise  2.26  Explain  why  the  transformation  of  X  to 
has  no  effect  on  the  variances  for  Plans  5  and  12. 

Exercise  2.27  From  the  data  presented  in  the  last  line  of 
Table  2.6,  compute  the  RelVar  of  the  ratio  estimator  (Plan  2) 
after  the  transformation,  that  is,  when  XT  =  X^  -  0.46  is  used 
as  an  auxiliary  variable  instead  of  X^.  Notice  that  a  trans¬ 
formation  of  this  kind  does  not  affect  the  variance  or  covariance , 
but  the  value  of  R  is  changed. 
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We  have  established  that  the  csa  of  a  branch  is  a  good 
measure  of  size  with  regard  to  probabilities  of  selection. 
Judging  from  Figure  2.4,  the  least  squares  line  and  the  ratio 
line  differ  enough  to  raise  a  question  of  transforming  Y  in¬ 
stead  of  X.  It  would  be  possible  to  add  a  constant  to  Y  so 
the  least  squares  regression  line  for  X  and  Y'  (where  Y'  =  Y+C) 
would  pass  through  the  origin  and  be  the  same  as  the  ratio  line 
involving  X  and  Y  .  Such  a  transformation  would  not  change  the 
selection  probabilities  and  it  appears  that  an  appreciable 
reduction  in  variance  might  be  obtained. 

Exercise  2.28  Find  the  value  of  C  in  Y'  =  Y+C  such  that 
the  regression  line  for  X  and  Y'  will  pass  through  the  origin. 
Ans .  C  =  2  4 . 

Consider  the  result  from  Exercise  2.28  and  the  transfor¬ 
mation  Y'=  Y  +  24.  The  estimator  of  Y,  without  stratification, 
would  be 

1  .  n  y  •  +24 

Y  =  „  (X)(Z  - )  -  24 

i  i 

The  RelVar  of  y  is  0.294  which  is  about  8  percent  less  than 
0.319,  the  RelVar  for  Plan  4.  If  it  is  necessary  to  estimate 
C  from  the  sample,  the  variance  of  C  must  be  taken  into  account 
and  there  is  no  gain  from  the  transformation.  However,  if 
there  is  good  prior  information  about  the  value  of  C,  the 
possibility  of  the  transformation  Y+C  might  be  worth  considering. 

When  considering  sampling  with  pps ,  look  for  a  measure  of 
size  that  is  close  to  being  proportional  to  Y^ .  If  Y^  is 
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and  the 


Y. 

exactly  proportional  to  X^,  the  ratio  — is  constant 
sampling  variance  is  zero.  Also,  the  situation  is  a  good  one 
for  sampling  with  pps  when  the  standard  deviation  of  Y  for  a 
fixed  value  of  X  is  proportional  to  the  value  of  X. 

In  stratified  random  sampling  some  statisticians,  in  the 
absence  of  a  better  basis,  allocate  a  sample  to  strata  accord¬ 
ing  to  estimates  of  the  proportion  of  the  total  that  each 
stratum  accounts  for.  For  example,  stratum  5  accounts  for 
23  percent,  the  aPPles  so  23  percent  of  the  sample 

would  be  allocated  to  stratum  5.  Only  8  percent  of  the  sample 
would  be  allocated  to  stratum  1  even  though  it  contains  21 
percent  of  the  branches.  In  practice,  prior  data  often  pro¬ 
vide  an  estimate  of  the  proportion  of  the  population  total  that 

ft 

each  stratum  accounts  for.  Thus,  the  size  of  sample,  n^ ,  for 
stratum  h  would  be 


IT  ft 

nu  =  Pu  n 
h  h 


f ! 


(2.24) 


where  is  the  proportion  that  stratum  h  accounts  for  and  n 
is  the  total  size  of  the  sample.  The  sampling  fraction  for 
stratum  h  is 


ft 


ft 


n 


and  for  n  =  1 


tf 
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"  rh 

Thus,  f^  =  compares  with  the  sampling  fractions  (or  selec- 
n 

tion  probabilities)  for  the  methods  discussed  previously. 

ft 

Values  of  f^  for  the  numerical  example  on  apples  are  in 
T  ab 1 e  2.9. 


ft 

Exercise  2.29  Verify  two  of  the  values  of  f ^ . 

Allocating  a  sample  to  strata  in  proportion  to  prior 
estimates  of  the  amounts  that  the  strata  account  for  is  a 
good  plan  where  the  coefficients  of  variations  of  Y  are  nearly 
constant  among  strata.  In  Table  2.8  notice  that  the  coeffi¬ 
cient  of  variations  tend  to  decrease  as  the  branch  size 
increases.  This  phenomena  appears  very  frequently.  The  fact 
that  the  first  stratum  has  the  highest  coefficient  of  variation 
means  that  the  sample  for  the  first  stratum  should  be  larger 

ff 

than  n^,  given  by  Eq.  2.24.  This  is  also  indicated  by  the 

ft 

comparison  of  f^  and  f^.  As  a  ’’rule  of  thumb,"  some  statis¬ 
ticians  have  adopted  a  practice  of  allocating  a  sample  according 
to  (2.24)  and  then  doubling  the  size  of  the  sample  for  the  first 
stratum  and  increasing  the  sample  for  the  second  stratum  by 
50  percent.  Small  departures  from  an  optimum  allocation  have 
a  negligible  impact  on  variance.  Moreover,  in  practice  an 
exact  optimum  allocation  cannot  be  achieved  because  exact 
values  of  within  stratum  variances  are  unknown. 
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Table  2.1- -Number  of  Apples  and  Cross  Section  Areas 
(Terminal  Branches) 


Table 

2 . 2 - -Relative  Variances  of  Numbers 
Primary  and  Terminal  Branches 

of 

Apples  Among 

Relative 

Variances  Among 

Plan 

Estimator 

Primary  Branches 

Terminal  Branches 

1 

Yi 

1.17 

0 . 660 

2 

y  2 

0 . 247 

0.382 

3 

ys 

0.186 

0.350 

4 

y^ 

0.194 

0.319 

Table 

2 . 3 - - Relat ive  Variances  of  Numbers  of  Apples 
in  Terms  of  One  Terminal  Branch 

Expressed 

Relative  Variances 

Among 

Plan 

Estimator 

Primary  Branches  Terminal  Branches 

1 

y  i 

5.639 

0.660 

2 

y2 

1.191 

0 . 382 

3 

ys 

0.897 

0 . 350 

4 

y4 

0.935 

0.319 
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Table  2 . 4 - -Relative  Variances  Expressed  as  a  Proportion  of 
the  Variances  for  Plan  1 


Plan  Primary  Branches  Terminal  Branches 


1 

1.00 

1.00 

2 

0.21 

0.58 

3 

0.16 

0.53 

4 

0.17 

0.48 

Table  2 . 5 - - Relat ive  Variances  for  Terminal  Branches  as  a 

Proportion  of  the  Variances  for  Primary  Branches 

Plan  Primary  Branches  Terminal  Branches 

1  1.00  0.12 

2  1.00  0.32 

3  1.00  0.39 

4  1.00  0.34 
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Table  2.6--Data  for  Six  Trees 
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Table  2 . 7 -- Summary  of  Sampling  Plans- -Estimators  and  Their  Relative 
Variances  for  a  Sample  of  One  Terminal  Branch 


Method  of 


1/ 


Plan  Sampling  - 


2/ 

Estimator  — 


Relative 
Variance 
for  n  =  1 


A 


mean 


Yi  =  y 


0.660 


A 


ratio 


y2-^ 


0.382 


regression 


y3  =  y  +  b  (X  -  x) 


0.350 


p.p.s 


y4 


y4  = 


0.319 


mean 


y6  =  y  =  N  ENhyh 


0.512 


separate  ratios 


J.  _  v  il 

y7  '  N  EXh  ~ 


0.279 


Nn 


separate  regressions  yft  =  s  rr-  [y,  +b,  (X,  -x,  )  ]  0.256 


combined  ratio 


8  N  L/h  ^hv  h  h 

-  V  ^ 

y9  =  y  - 

X 

S 


0.307 


10 


combined  regression  y, „  =  yc  +  b  (X  -  x  ) 


0.294 


11 


p.p.s 


yn-  In  Vxi 


0.241 


If  A  Simple  random  sampling  without  stratification 

B  Sampling  with  replacement  and  p.p.s  without  stratification 
C  Simple  random  sampling  within  strata  (trees) 

D  Sampling  with  p.p.s  within  strata  (trees) 

2/  h  is  the  index  to  strata 

s  as  in  y  refers  to  stratification.  Thus  ys  is  the  mean  of  a  stratified  random 
sample  ancl  y  is  the  mean  of  a  simple  random  sample. 
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Table  2.8--Bata  for  Si^e -of -Branch  Strata 


Stratum 

Branch  Size 

No .  of 
Branches 

*h 

No .  of 
Apples 

?h 

SYh 

1 

.60-1.00 

28 

0.816 

586 

20.93 

16.79 

0.80 

2 

1.01-1.40 

45 

1.161 

1500 

33.33 

17.22 

0.52 

3 

1.41-1.80 

28 

1.543 

1765 

63.04 

31.42 

0.50 

4 

1.81-2.20 

21 

1.935 

1488 

70.86 

39.53 

0.56 

5 

2.21  + 

13 

2.550 

1634 

125.69 

52.72 

0.42 

Total 

or  Average 

135 

1.423 

6973 

51.65 

Stratum 

Table  2.9 

fh 

--Sampling  Fractions 

-  1 1 

P  P  £ 

n  h  rh 

1 

.  0046 

.0042 

.0027 

.  0030 

2 

.  0047 

.  0060 

.  0054 

.  0048 

3 

.  0086 

.  0080 

.0083 

.0090 

4 

.0109 

.  0101 

.0113 

.  0102 

5 

.  0145 

.  0133 

.0161 

.0180 

Table  2.10-- 

Effect  of 

Transformation 

Relative  Sampling  Variance 

Estimator 

Plan 

Before 

Transformation 

After 

Transformation 

Ratio(no  stratification) 

2 

0.382 

0.350 

Regress  ion (no  stratification) 

3 

0.350 

0.350 

PPS  (no  stratification) 

4 

0.319 

1.403 

Mean  (stratification  by  size) 

12 

0.273 

0.273 
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Figure  2.1- -Map  of  apple  tree  no.  3  showing  branch 
identification  and  number  of  apples 
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03 

c/) 


Number  of  apples 
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Figure  2. 2- -Number  of  apples  by  csa — terminal  branches 


Number  of  apples  Number  of  apples  Number  of  apples 


Figure  2.3- 


-Dot  charts  and  ratio  lines  by  trees 
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Figure  2 .4- -Size-of -branch  strata,  least  squares  regression 
line,  and  line  for  pps. 


CHAPTER  III 


RANDOM- PATH  SAMPLING  OF  FRUIT  TREES 

3.1  INTRODUCTION 

The  methods  discussed  in  Chapter  II  required  a  map  of 
each  tree  to  be  sampled.  A  map  of  a  tree  provides  a  good 
sampling  frame,  but  drawing  a  map  and  measuring  the  csa's  of 
all  branches  is  too  time  consuming.  In  the  research  for  prac¬ 
tical  ways  of  probability  sampling,  photographs  of  trees  taken 
when  the  trees  had  no  leaves  have  been  studied  for  possible 
use  as  sampling  frames,  including  the  estimation  of  branch 
sizes  for  sampling  with  pps .  Photography  has  also  been  con¬ 
sidered,  in  the  context  of  double  sampling,  as  a  means  of 
counting  and  estimating  the  number  of  fruit  on  the  tree. 

In  this  chapter  the  random-path  method  proposed  by  Jessen 
for  sampling  fruit  on  a  tree  will  be  illustrated  using  one  of 
the  six  apple  trees  in  the  analysis  in  Chapter  II.  Two  random 
path  methods  will  be  compared  with  two  methods  that  were  dis¬ 
cussed  in  Chapter  II.  The  comparisons  will  be  made  as  though 
only  one  terminal  branch  from  a  tree  is  to  be  selected  and 
used  to  estimate  the  total  number  of  fruit  on  the  tree. 


6/  Jessen,  Raymond  J. ,  Determining  the  Fruit  Count  on  a  Tree 
by  Randomized  Branch  Sampling,  Biometrics,  March  1955. 


3.2  FOUR  METHODS  OF  SAMPLING  A  TREE 


The  four  sampling  methods  and  estimators,  as  described  in 
this  section,  include  only  apples  that  are  on  terminal  branches. 

A  small  proportion  of  the  apples  are  not  on  terminal  branches. 
Methods  of  including  these  apples  will  be  discussed  later. 

(1)  The  first  method  is  included  primarily  as  a  base  for 
comparison.  It  is  the  same  as  plan  1  discussed  in  Chapter  II. 

After  all  terminal  branches  on  a  tree  have  been  identified  and 
numbered,  one  branch  is  selected  at  random  with  equal  probability. 
Apples  on  the  sample  branch  are  then  counted.  The  estimator  is 
Ny^  where  N  is  the  number  of  branches  on  the  tree  and  y^  is  the 

number  of  apples  on  the  sample  branch.  Since  y^  refers  to  a 

sample  value,  one  would  expect  i  to  be  the  index  for  branches 
in  a  sample .  But,  we  are  considering  a  sample  of  only  one  branch 
and  it  is  convenient  to  let  i  be  the  index  to  branches  in  the 
population.  Thus,  if  the  5th  branch  in  the  population,  ,  Y^, 

...,  Y^ ,  happens  to  be  selected,  y^  =  Y^.  The  first  method  will 
be  referred  to  as  DS-EP,  which  means  direct  selection  from  a  list 
of  all  branches  with  equal  probabilities. 

(2)  The  second  method  is  a  random-path  technique.  Beginning 
from  the  bottom  of  the  tree,  the  primary  branches  are  all  identified 
and  one  of  the  primary  branches  is  selected  at  random  with  equal 
probabilities.  The  selected  primary  branch  is  then  examined  to 
identify  all  second-stage  branches  from  it.  Then,  one  second-stage 
branch  is  selected  at  random  with  equal  probabilities.  The  process 
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is  discontinued  when  a  terminal  branch  has  been  selected.  The 

estimator  is  —  where  p.  is  the  probability  of  selecting  the 
Pi  i 

particular  terminal  branch  that  happens  to  be  selected.  As  a 
short  title  for  this  method  RP-EP  will  be  used  where  RP  represents 
random  path  and  EP  means  equal  probability  of  selection  at  each 
stage  of  branching. 

(3)  Like  the  first  method,  the  third  requires  a  complete 
identification  of  all  terminal  branches  prior  to  selection.  The 
csa  (cross  sectional  area)  of  each  terminal  branch  is  measured  and 
one  branch  is  selected  with  pps ,  probability  proportional  to  csa. 

y. 

The  estimator  is  X — ,  where  x.  is  the  csa  of  the  selected  branch 

i 

and  X  is  the  sum  of  the  csa's  of  all  terminal  branches  on  the  tree. 
DS-PPS  is  the  short  title  for  this  method,  which  is  the  same  as 
plan  2  in  Chapter  II. 

(4)  The  fourth  method  is  a  random-path  method  which  differs 
from  method  two  in  the  probability  of  selection.  At  each  stage 
of  branching  the  csa’s  of  the  branches  at  that  stage  are  measured 


Xi 

and  one  branch  is  selected  with  pps.  The  estimator,  — ,  is  like 

^i 

the  estimator  for  the  second  method  but  the  values  of  P^  are 
different.  This  method  is  titled  RP-PPS. 

3.3  BRANCH  IDENTIFICATION  AND  DESCRIPTION  OF  DATA 

Data  for  tree  No.  3,  which  was  represented  in  Figure  2.1, 
are  presented  in  Table  3.1  in  a  way  that  shows  the  stage  of 
branching.  There  were  only  three  primary  branches.  Their  csa's 
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11.60,  13.45,  and  12.84,  are  presented  in  the  column  titled  csa 
under  1st  stage.  The  sum  of  these  csa's  is  37.89,  Thus,  if  a 
primar^  branch  is  selected  with  pps ,  the  first  primary  branch 
would  have  a  selection  probability  equal  to  ^  ^  .  For  further 

illustration  of  the  recording  system,  notice  that  the  second  digit 
of  the  identification  number  shows  four  second-stage  branches  from 
the  second  primary  branch.  The  csa’s  of  these  four  branches  and 
their  sum  are  recorded  under  2nd  stage.  This  scheme  of  branch 
identification  and  recording  is  continued  until  a  terminal  branch 
is  reached.  When  this  occurs,  the  number  of  apples  on  a  terminal 
branch  is  recorded  to  the  right  of  its  csa.  Thus,  Table  3.1  shows, 
for  example,  that  branch  2-3  was  a  terminal  branch  with  a  csa  equal 
to  1.99  square  inches  and  that  it  had  124  apples  on  it.  The  numbers' 
in  parenthesis  are  numbers  of  "path"  apples  which  will  be  discussed 
later . 

Terminal  branches  were  defined  as  branches  having  a  csa 
between  3/4  and  2  square  inches.  Adherence  to  exact  size  is  not 
possible.  For  example,  the  first  terminal  branch  1-1-1  has  a  csa 
equal  to  2.68  and  is  large  enough  so  an  additional  stage  of  branch¬ 
ing  was  probably  considered.  If  from  1-1-1  there  were  two  branches 
of  about  equal  size,  those  two  branches  would  have  been  terminal 
branches.  Probably  1-1-1  divided  into  several  branches  that  were 
too  small  to  be  considered  as  terminal  branches.  As  another  case, 
suppose  at  the  last  stage  of  branching  there  is  a  branch  with  a 
csa  of  1.5  square  inches  which  is  dying  and  clearly  has  no  fruit. 
This  branch  could  be  shown  on  the  map  but  marked  for  exclusion  and 
not  counted  as  a  branch  for  sampling  purposes. 
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Along  a  path  from  the  base  of  a  tree  to  a  terminal  branch 
there  are  some  branches  which  are  much  too  small  to  qualify  as 
terminal  branches.  Fruit  on  small  branches  along  a  path  to  a 
terminal  branch  have  been  called  path  fruit  and  must  be  ac¬ 
counted  for  in  some  way.  For  example,  on  the  path  from  the 
base  of  1-2  to  the  bases  of  1-2-1  and  1-2-2  there  were  three 
apples.  These  three  apples  are  recorded  in  Figure  2.1  and  in 
Table  3.1  next  to  the  csa  of  branch  1-2.  The  counts  of  apples 
along  the  paths  are  shown  in  parenthesis  in  Table  3.1  to  dis¬ 
tinguish  such  apples  from  apples  on  the  terminal  branches. 

There  are  various  ways  of  dealing  with  the  path  fruit;  but, 
first  let  us  examine  the  probability  of  selecting  any  given 
terminal  branch  with  regard  to  each  of  the  four  methods. 

3.4  PROBABILITY  OF  SELECTION  AND  ESTIMATION 

With  the  DS-EP  method  each  one  of  the  26  terminal  branches 

1  "t  hi 

has  a  selection  probability  equal  to  jg-*  With  DS-PPS  the  i 

X. 

terminal  branch  has  a  probability  of  selection  equal  to  y—  where 
is  its  csa  and  X=EX^.  Calculating  the  probabilities  for  the 
random-path  methods  is  more  involved.  For  example,  consider 
terminal  branch  1-2-1-2  and  RP-EP.  In  Table  3.1  notice  that 
the  numbers  of  branches  at  each  stage  on  the  path  to  terminal 
branch  1-2-1-2  are  3,  5,  2,  and  2.  Thus,  the  probability  of 
selecting  branch  1-2-1-2  is  (y)  (y-)  (y)  (y)  =  which  is 

the  product  of  the  probabilities  of  selection  at  the  four 
stages.  With  RP-PPS  the  product  of  corresponding  probabilities 
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11.60 


5.61,  ,4.13,  ,2.32 


at  the  four  stages  is  (37790")  (74-794)  (5795-)  (3T80')  =  -04862 

An  estimator  for  each  of  the  four  methods  was  presented 
above.  However,  all  of  the  four  estimators  can  be  written  in 
the  same  form, 

y,- 


Yi  =  p7 


(3.1) 


where  Y.  is  the  estimator  and  i  =  1,  2 . N  is  the  index 

t  h 

to  terminal  branches  on  the  tree.  If  the  1  branch  happens 

yi 

to  be  selected,  then  —  is  the  estimate  of  the  total  number 

Pi 

of  apples  on  the  tree  (except  that  path  apples  are  not  in- 

t  h 

eluded)  where  y^  is  the  number  of  apples  on  the  1  terminal 
branch  and  p^  is  the  probability  of  selecting  it.  The  value 
of  p^  depends  on  the  method  of  sampling.  In  fact  there  are 
four  sets  of  probabilities,  one  for  each  method.  These  four 
sets  of  values  are  presented  in  Table  3.2  in  the  columns  headed 
Pi ,  ?2’  an<4  ?4«  The  columns  headed  Y^,  ,  Y^,  and  Y^  con¬ 

tain  estimates  of  the  total  number  of  apples  on  the  tree 
depending  upon  the  terminal  branch  that  happens  to  be  selected. 
A  discussion  of  these  estimates  follows. 

Exercise  3.1  Refer  to  Table  3.1  and  for  each  of  the  four 
sampling  methods  compute  the  selection  probabilities  for  termi¬ 
nal  branches  2-3  and  3-2-4.  Compare  your  answers  with  the 
probabilities  presented  in  Table  3.2. 

To  include  the  path  apples  there  are  at  least  two  possi¬ 
bilities  that  might  be  considered  with  regard  to  the  DS-EP  and 
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DS-PPS  methods:  (1)  Count  all  path  apples  and  add  the  count  to 
the  estimate  of  the  total  number  of  apples  on  terminal  branches 
that  is  obtained  from  a  sample  of  terminal  branches;  (2)  Define 
sampling  units  that  are  sections  of  path  between  the  base  of 
the  tree  and  the  terminal  branches.  Then  select  a  sample  of 
such  sections  to  estimate  the  path  fruit. 

With  the  random-path  methods,  it  is  necessary  to  count  the 
apples  on  each  section  of  the  path  along  the  path  to  a  terminal 
branch.  Also,  it  is  necessary  to  determine  the  probability 
that  each  section  of  the  path  had  of  being  in  the  sample. 

Since  the  RP-PPS  method  is  of  primary  interest,  it  will 
be  used  to  illustrate  how  the  path  fruit  can  be  accounted  for 
in  the  estimation  process.  Consider  the  three  apples  (see 
Table  3.1)  on  the  path  section  1-2  which  is  the  section  between 
the  base  of  1-2  and  the  third  stage  branches.  To  find  the 
probability  of  these  three  being  in  the  sample,  consider  re¬ 
peated  application  of  a  random-path  sampling  method.  These 
three  apples  will  be  in  the  sample  whenever  this  path  section 
is  traversed.  Therefore,  under  the  RP-PPS  method,  the  proba¬ 
bility  of  this  path  section  being  in  the  sample  is 
( \j[  g-Q )  ( p4  ]  9  4  ^  =  .1150  which  gives  -  =  26.1  apples  to 

be  included  in  the  estimate  whenever  this  path  section  is 
traversed. 

It  is  important  to  observe  that  a  random  path  always  ends 
with  one  and  only  one  terminal  branch.  There  are  three  termi¬ 
nal  branches,  1-2-1-1,  1-2-1-2,  and  1-2-2  that  are  connected 
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to  the  path  section  1-2,  and  26.1  would  be  included  in  the 
estimate  that  is  made  from  the  selected  terminal  branch  that 
follows  path  section  1-2.  There  are  no  path  apples  other  than 
the  three  that  have  been  mentioned  between  the  base  of  the  tree 
and  the  three  terminal  branches  that  follow  path  section  1-2. 
Therefore,  the  estimated  total  number  of  apples  in  the  event 
any  one  of  the  three  terminal  branches  is  selected  would  be: 


Terminal  Branch 


Estimate 


3 

73 

=  2379 

.  1150 

.03103 

3 

138 

=  2863 

.1150 

.04864 

3 

133 

=  3794 

.  1150 

.03530 

With  the  application  of  either  random-path  method,  and 
assuming  path  fruit  are  recorded  for  each  path  section  that 
is  traversed,  an  estimator  that  includes  path  fruit  can  be 
written  in  generalized  form.  It  appears  to  be  complicated, 
but  an  illustration  follows  that  should  help  clarify  it.  The 
estimator  is: 


y 


y 


ki 


y 


t  i 


Y  =  l£i  +  +  .  *1  +  +  /  tl 

i  Poi  CPoi)---(Pki)  (PoP-'-lPtf 


(3.2) 


where 

i  =  1, 

2, 

.  .  . ,  N  is 

an  index 

of 

the 

terminal  branches 

k  =  0, 

1, 

.  .  .  ,  t  is 

an  index 

of 

the 

path  sections  of 

h 

the  path  to  the  i  terminal  branch  (t  is  not 
constant,  it  depends  on  i)  , 
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■f-  Vi 

yki  =  the  number  of  path  fruit  on  the  k  path  section 
of  the  path  from  the  base  of  the  tree  to  the  1 
terminal  branch, 

Pj^  =  the  conditional  probability  of  selecting  the 

t  l"i  t  h 

k  path  section  of  the  path  to  the  1  terminal 

branch,  given  that  the  preceding  path  section  has 

been  selected. 

When  k  =  0,  the  path  section  referred  to  is  the  part  of 

the  tree  between  ground  level  and  the  bases  of  the  primary  or 

first  stage  branches.  Given  that  the  tree  is  in  the  sample, 

pQ^  =  1  which  means  that  this  path  section  is  a  part  of  the 

path  to  every  terminal  branch.  Generally,  yQ^>  the  number  of 

fruit  on  this  section  of  the  tree,  will  be  zero.  When  k  =  t, 

the  k^*1  path  section  becomes  a  terminal  branch,  thus  y  ^  is 

th 

the  number  of  apples  on  the  1  terminal  branch,  and  p  ^  is 

t  h 

the  conditional  probability  of  selecting  the  i  terminal 
branch  given  that  the  path  section  that  it  is  connected  to 
has  already  been  selected. 
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Suppose  application  of  the  RP-PPS  method  leads  to  terminal 
branch  1-2-1-1.  In  this  case,  k  =  0,  1,  2,  3,  4  and  the  values 

of  yki  and  pki  are: 


Path  section 


No.  of  apples 
on  path  section, 


yki 


Conditional 
probability,  p^.- 


0 


0 


1 


1 


0 


2 


3 


3 


0 


4 


73 


11 , 

.60 

37, 

.  89 

5  , 

.61 

14, 

.94 

4, 

.  13 

5, 

.96 

1 , 

.48 

3. 

.  80 

.  3061 


.  3755 

.6930 


.  3895 


Substituting  these  values  of  y^  and  p^  in  the  estimator, 

(3.2) ,  gives : 

(°)  +  (0)  +  C3)  +  CO) 

ITT  (1) (.3061)  (1) (.3061) (.3755)  (1) (.3061) (. 3755) (.6930) 


+  (73)  _  7379 

(1) (. 3061) (. 3755) (.6930) (. 3895)  ^ 


Exercise  13.2  For  the  RP-EP  and  RP-PPS  methods  use  the 
estimator  (3.2)  to  obtain  estimates  corresponding  to  terminal 
branches  3-1-2,  3-1-4- 1 3  and  3-3.  Your  results  should  agree 

with  the  estimates  presented  in  Table  3.2. 
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There  is  an  alternative  view  of  the  above  method  of  in¬ 


cluding  the  path  fruit  which  leads  to  the  same  answers.  The 
idea  is  to  prorate  the  path  fruit  to  the  terminal  branches  that 
follow  the  sections  of  the  path  where  the  path  fruit  are  found. 

The  prorating  is  done  according  to  the  probabilities  of  selection. 
Under  the  RP-EP  method  the  three  apples  on  the  path  section  1-2 
would  be  prorated  as  follows: 


Prorated  amount 


Terminal  branch 


(i)(|)(3)  =  .75 


1-2-1-1 


(|)(i)(3)  =  .75 


1-2-1-2 


1-2-2 


1.50 

3.00 


Total 


Notice  that  (^) (j) ,  (j) (j) ,  and  ( j )  are  the  conditional 


probabilities  of  selecting  one  of  the  three  terminal  branches, 
given  that  the  path  section  1-2  has  already  been  selected.  The 
conditional  probabilities  add  to  1  which  verifies  that  the  method 
of  prorating  accounts  for  all  of  the  path  fruit. 

If  1-2-1-1,  for  example,  is  the  selected  terminal  branch, 

.75  is  added  to  73,  the  number  of  apples  on  1-2-1-1.  The  estimate 
of  the  total  number  of  apples  on  the  tree  is  then  obtained  by 
dividing  73.75  by  the  probability  of  selecting  1-2-1-1  which  is 

60  ’  Thus>  (60)  (73.7  5)  =  4425.  The  branch  total, 

73.75,  (including  the  prorated  amount)  appears  in  Table  3.1  in 
the  column  titled  EP ,  and  the  expanded  total,  4425,  appears  in 
Table  3.2  in  the  column  titled  Y^- 
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Under  the  RP-PPS  method,  the  system  of  prorating  is  the 
same  but  the  probabilities  are  different.  Thus, 


Terminal  branch 

Prorated  amount 

1  -  2  - 1  - 1 

.81 

1  -  2  - 1  -  2 

<3>  ■ 

1.27 

1-2-2 

rl-83H3) 
l5. 96J  1  J 

.  92 

Total 

3.00 

The 

estimator 

T,  Eq.  (3.2), 

can  be  written  in 

a  form 

that  corresponds  to  the  idea  of  prorating  path  fruit 

to  terminal 

branches . 

Let  p . 

r  l 

(P0i) • • •  (Pt i^  ’ 

which  is  the  probability  of 

selecting 

the  ith 

terminal  branch. 

It  follows  that 

v  -  Ii 

i  Pi 

(3.3) 

where  y. 

1  i 

=  [(PU) 

•••(pti)y0i]+--- 

+  I(p(k+ild  '  '  '  (pti 

) yki ] + • • -+ [y 

t  h 

Thus,  y^  is  the  number  of  fruit  "on"  the  i  terminal  branch  in¬ 
cluding  prorated  amounts  of  path  fruit.  Assuming  the  RP-PPS  method 


and  terminal  branch  1-2-1-1  as  an  example,  the  value  of  y^  is 

(|^|)  (3)  +  73  -  73.81  and  Y^  is  =  2379  which  gives 

the  same  result  that  was  obtained  when  Eq.  (3.2)  was  used. 

Table  3.2,  columns  headed  and  ,  present  estimates  of 

the  total  number  of  apples  on  the  tree  for  the  RP-EP  and  RP-PPS 
methods  and  each  of  the  possible  random  paths.  These  estimates 
were  obtained  by  using  the  technique  of  prorating  path  fruit, 
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Eq.  (3.3).  That  is,  estimates  of  the  total  number  of  apples 
were  obtained  by  dividing  the  values  of  (last  two  columns  of 

Table  3.1)  by  the  appropriate  probabilities  which  are  presented 
in  Table  3.2,  columns  and  . 

For  comparison  of  the  four  methods  we  now  need  to  decide 
how  to  include  the  path  fruit  for  the  DS-EP  and  DS-PPS  methods. 
If  the  amount  of  path  fruit  is  small,  the  best  method  might  be 
to  count  all  path  fruit  at  the  time  the  tree  is  mapped  to  deter¬ 
mine  terminal  branches.  In  this  case,  assuming  a  sample  of  one 
terminal  branch,  the  estimator,  would  be 


! 


(3.4) 


'  i 

where  Y  is  the  number  of  path  fruit,  y^  is  the  number  of  fruit 
t  h 

on  the  i  terminal  branch  and  p  is  the  probability  of  selecting 

i 

t  h 

the  i  terminal  branch.  Alternatives  are  not  considered  in  this 
illustration  because,  from  a  practical  viewpoint,  interest  is  in 
the  random  path  methods.  Thus,  as  a  matter  of  expediency,  the 
estimator  (3.4)  was  used  to  obtain  the  estimates,  Y^  and  Y^  ,  that 

are  presented  in  Table  3.2  for  the  DS-EP  and  DS-PPS  methods. 

Since  only  51  apples  out  of  1901  were  on  path  sections,  the  method 
of  accounting  for  the  apples  on  path  sections  probably  has  a 
very  small  impact  on  the  sampling  variance. 
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Exercise  3.3  For  terminal  branches  3-1-4-1  and  3-3 


calculate  estimates  of  the  total  number  of  apples  on  the  tree 
for  the  DS-EP  and  DS-PPS  methods  using  the  estimator  (3.4). 

Your  answer  should  agree  with  the  estimates  that  are  presented 
in  Table  3.  1  for  these  two  branches. 

For  each  terminal  branch  and  each  of  the  four  estimators 
(methods)  there  is  a  unique  estimate  of  the  total  number  of  apples. 
All  four  estimators  are  unbiased.  By  definition,  an  estimator  is 
unbiased  if  the  expected  (average)  value  of  the  estimates  that 
might  occur  is  equal  to  the  population  value.  To  find  the  expected 
value  of  an  estimator,  each  estimate  must  be  weighted  by  the 
probability  of  its  occurence. 

Exercise  3.4  For  the  RP-EP  and  RP-PPS  methods ,  compute  the 
expected  value  of  the  estimates  presented  in  Table  3.2.  The 
answer ,  except  for  rounding  error 3  should  be  exactly  1901s  which 
is  the  total  number  of  apples  on  the  tree. 

3.5  VARIANCES  OF  THE  ESTIMATORS 


With  reference  to  the  theory  of  expected  values,  the 
variance  of  a  random  variable,  Y,  is  the  average  of  the  squared 
deviations  of  Y  from  its  expected  (average)  value.  To  be  more 
specific,  suppose  Y  is  a  random  variable  that  can  equal  one  of 
a  set  of  values  Y^,  Y^j-.-jY^  with  probabilities  > • • • >  ^ 

where  =  1.  By  definition,  the  average  value  of  Y  is 

Y  =  E  (Y)  =  ZPiYi 
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and  the  variance  of  Y,  which  is  the  average  value  of  (Y-Y)2,  is 

-  2  N  -  2 
E  (Y-Y) z  =  EP.(Y.-Y)Z 

Exevoise  3.5  Show  that  EP.(Y.-Y)2  =  EP-Y.2-Y2 

i  i  li 

Consider  the  estimator  for  the  RP-PPS  method.  It  is  a  random 
variable  that  can  equal  any  one  of  the  set  of  values  in  column 
Y^  of  Table  3.2.  The  set  of  probabilities  is  presented  in  column 
P^.  By  definition,  the  variance  of  the  estimator  (or  estimates) 
is 

(.054  92)  (3 7 51 -1901)  2 +.. .  +  (.0  6163) (8 14  -1901)  2  =  8 00,1 94 
or  using  the  right  hand  side  of  the  equation  in  exercise  3.5, 

(.  05492)  (37S1)  2+. . .  +  (.06163) (814)  2 - (1901)  2  =  800,194 

The  result,  800,194,  is  the  sampling  variance  for  the  RP-PPS 
method  when  only  one  terminal  branch  is  selected.  If  four  ter¬ 
minal  branches  (or  random  paths)  were  selected  with  replacement, 
four  estimates  of  the  tree  total  would  be  computed,  one  for  each 
branch,  and  the  variance  of  the  average  of  the  four  estimates 


wou 


Id  be  800 ^ 194  =  200,048 


The  sampling  variances  (for  a  sample  of  one  branch)  are 
presented  in  Table  3.3  for  each  of  the  four  methods  and  each  of 
the  six  trees.  The  third  tree  is  the  one  that  was  used  above 
as  an  example.  It  is  not  expected  that  the  four  methods  will 
always  rank  in  the  same  order  from  one  tree  to  another.  However, 
the  results  illustrate  some  points  that  are  of  interest  and  impor¬ 
tance  . 
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3.6  DISCUSSION  OF  THE  METHODS 


The  RP-EP  method  requires  considerably  less  time  than  the 
RP-PPS  method,  but  it  has  relatively  high  sampling  variance  be¬ 
cause,  at  any  given  stage  of  branching,  a  large  branch  has  the 
same  probability  of  selection  as  a  small  one.  That  is,  the 
RP-EP  method  is  such  that  the  probability  of  selecting  a  termi¬ 
nal  branch  has  little  or  no  relation  to  the  number  of  fruit  on 
the  branch.  The  result,  as  shown  by  the  sampling  variances  in 
Table  3.3,  is  a  good  illustration  of  a  point  that  was  made 
earlier.  Compared  to  selecting  sampling  units  with  equal  prob¬ 
ability  (as  in  the  DS-EP  method) ,  the  introduction  of  unequal 
probabilities  of  selection  (as  in  the  RP-EP  method)  will  in¬ 
crease  the  sampling  variance  unless  the  selection  probabilities 
are  related  to  the  values  of  the  characteristic  being  measured 
in  a  way  that  will  reduce  sampling  variance. 

Figure  3.1  is  a  dot  chart  with  the  number  of  apples  on  a 
branch  (column  headed  EP  in  Table  3.1)  plotted  against  the  values 
of  P^ •  The  wide  range  in  the  selection  probabilities  and  the 
lack  of  a  relation  explains  the  high  sampling  variance  of  the 
RP-EP  method  compared  with  the  other  methods.  For  comparison, 
Figure  3.2  is  a  dot  chart  for  number  of  apples  and  the  selection 
probabilities  for  the  RP-PPS  method.  Compare  Figure  3.2  with 
Figure  1.2  which  showed  a  dot  chart  where  sampling  with  pps 
would  rank  high. 

After  a  branch  has  been  identified  and  marked,  the  time 
required  to  obtain  its  csa,  with  a  convenient  instrument  that 
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gives  a  reading  directly  in  square  inches  (or  square  centi¬ 
meters),  is  quite  small.  The  use  of  csa  as  an  auxiliary  vari¬ 
able  reduced  sampling  variance  by  a  large  amount.  The  reduction 
in  variance  in  relation  to  cost  is  definitely  advantageous. 
According  to  Table  3.3,  the  sampling  variances  for  DS-PPS  and 
RP-PPS  are  about  the  same  and  much  less  than  the  sampling  vari¬ 
ance  for  the  DS-EP  method.  This  indicates  that  RP-PPS  is  a  good 
choice  because  it  avoids  the  work  of  identifying  all  terminal 
branches  before  sampling.  However,  results  in  Table  3.3  should 
not  be  accepted  as  representative.  The  csa  is  not  always  an 
effective  measure.  Pruning  and  maintenance  practices,  age  of 
trees,  species  or  variety  of  trees,  and  other  factors  have  some 
influence  on  the  relation  between  csa  and  number  of  apples. 

The  purposes  of  an  intensive  investigation  limited  to  a  few 
trees  include  testing  different  procedures  for  counting  apples 
or  measuring  the  size  of  branches,  and  acquiring  ideas  that 
seem  to  be  worth  exploring  as  possibilities  for  large  scale 
application . 

It  is  extremely  important  in  the  processes  of  sampling  to 
understand  the  part  played  by  randomization.  Important  biases 
sometimes  occur  even  when  strict  attention  is  paid  to  details 
in  making  random  selections.  On  the  other  hand,  subjective 
evaluations  or  determinations  in  sampling  are  commonplace.  With 
knowledge  of  how  various  factors  effect  sampling  variance,  the 
exercise  of  good  judgement  can  be  very  effective  in  reducing 
sampling  variance.  But,  there  are  points  in  the  processes  of 


99 


sampling  where  a  determination  should  be  strictly  random.  Some 
design  constraints  may  be  determined  subjectively  but  selections 
of  units  for  a  sample  should  be  in  accord  with  rigorous,  technical 
interpretation  of  randomness.  It  is  generally  preferable  to  have 
random  selections  made  under  competent  supervision  in  an  office, 
but  that  is  not  always  feasible.  Thus,  one  advantage  of  taking 
photographs  of  a  sample  of  bare  trees  (assuming  it  is  feasible) 
is  that  sample  branches  can  be  selected  in  the  office.  The  se¬ 
lected  branches  are  marked  on  photographs  for  enumerators.  In  this 
situation  an  enumerator’s  work  is  subject  to  full  verification. 
Incidentally,  the  economics  of  sample  surveys  suggests  that  larger 
investments  in  sample  design  and  selection  can  often  be  justified 
when  the  same  sample  is  to  be  used  for  several  surveys  rather 
than  one. 


Exercise  3.6  Suppose  the  RP-PPS  method  is  being  applied 
and  in  the  process  you  come  to  the  following  situation: 


has  already  been  selected.  It  divides  into  two  branches  A-l  and 
A- 2  with  csa’s  equal  to  1.4  and  1.6.  With  regard  to  size 3  the  two 
branches ,  A-l  and  A-2 ,  qualify  as  terminal  branches  and  ordinarily 
A-l  and  A-2  would  be  accepted  as  terminal  branches.  But 3  before 
selecting  one  of  the  two ,  you  happen  to  notice  that  A-2  has  no 
apples  on  it  and  that  A-l  appears  to  have  approximately  an  average 
amount.  Consider  the  following  alternatives : 
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(1) 


Accept  A 3  which  includes  A-l  and  A-2  3  as  the  terminal 

branch 3  and  expand  the  count  of  apples  by  \  3  where  PA 

A  A 


was  the  probability  of  selecting  A. 

(2)  Accept  A-l  and  A-2  as  terminal  branches  and  select  one 

1  3.0 

with  pps.  Expand  the  count  on  A-l  or  A-2  by  (p  )  (yyr) 

1  3.0 

or  (p  )  (.yL-r)3  depending  on  whether  A-l  or  A-2  is  selected. 


(3) 


Discard  A-2  since  it  has  no 
as  the  terminal  branch  using 


apples 

Cp  M 


on 


3.0 

1.4 


) 


it  and 
as  the 


take  A-l 
expansion 


factor. 

Discuss  the  alternatives  with  regard  to  bias  and  sampling  variance. 

Exercise  3.7  Refer  to  exercise  3.6  and  as  a  variation  of  the 
situation  assume  that  branch  A-2  has  been  selected  at  random  in 
accord  with  the  instructions  for  the  random  path  method.  The 
enumerator  prepares  to  count  the  apples  on  A-2  but  finds  there  are 
no  apples.  He  recognizes 3  since  a  sample  of  only  one  branch  is  to 
be  selected  for  the  sample  from  this  tree 3  that  the  estimate  of 
the  number  of  apples  on  the  tree  will  be  zero  (assuming  no  path 
fruit  on  thepath  to  A-2).  There  is  obviously  a  large  number  of 
apples  on  the  tree 3  so  he  might  have  a  strong  opinion  that  some¬ 
thing  should  be  done  that  would  give  a  better  sample.  How  would 
you  respond  to  each  of  the  following  possibilities: 

( 1)  Accept  A-2  as  a  terminal  branch3  which  means  using  zero 
as  an  estimate  of  the  number  of  apples  on  the  tree. 
Remember  A-2  has  already  been  selected. 
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(2)  Reject  A-2  as  a  sample.  Start  at  the  beginning  and 
select  another  terminal  branch  to  replace  A-2. 

(3)  Accept  A  which  includes  A-l  and  A-2 3  as  the  terminal 
branch  for  the  sample. 

Discuss  the  three  possibilities  with  regard  to  bias  and  sampling 
variance . 

Exercise  3.8  In  application  of  the  RP-PPS  method  would  it 
be  advisable  to  be  looking  forward 3  as  one  approaches  the  terminal 
branch  stage3  for  branches  that  are  large  enough  to  be  terminal 
branches  but  clearly  have  a  very  small  number  of  apples  on  them. 
With  reference  to  the  diagram  in  exercise  3.6  as  an  example 3  an 
enumerator  looking  forward3  and  considering  what  was  ahead3  could 
have  stopped  when  A  was  selected  and  accepted  A  as  a  terminal 
branch.  Otherwis e 3  he  would  normally  have  followed  the  selection 
procedure  one  stage  further .  In  application  of  the  random  path 
method3  what  is  your  opinion  of  the  feasibility  of  looking  ahead 
and  taking  eye  estimates  of  numbers  of  apples  into  account  in 
determining  the  terminal  branch.  Can  it  be  used  to  reduce  sampling 
error  without  risk  of  introducing  bias?  Think  about  the  matter 
with  regard  to  instructions  that  would  be  given  to  enumerators. 

Exercise  3.9  It  is  not  likely  that  there  would  be  an  interest 
in  estimating  the  average  number  of  terminal  branches  per  tree. 
However 3  as  an  exercise  3  suppose  the  RP-PPS  method  is  applied  to 
the  tree  for  which  data  arepresented  in  Table  3.1.  Assume  that 
the  following  four  terminal  branches  are  selected  as  a  sample: 
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1-1-2,  1-2-1-2,  2-4 ,  and  3-2-1.  From  this  sample,  estimate 

the  number  of  terminal  branches  on  the  tree.  (The  selection 
probabilities  have  already  been  computed,  see  Table  3.2).  The 
parameter  being  estimated  is  26.  Ans .  33.4. 

Exercise  3.10  Suppose  a  sample  of  25  apple  trees  has  been 
selected  and  that  four  enumerators  have  been  trained  in  the  appli¬ 
cation  of  the  RP-PPS  method.  Assume  that  each  enumerator ,  working 
independently  and  using  the  RP-PPS  method,  selects  a  sample  of  one 
terminal  branch  from  each  of  the  25  trees.  It  is  unlikely  that 
enumerators  will  interpret  terminal  branches  in  exactly  the  same 


way.  For  example ,  one  enumerator  might  have  a  tendency  to  follow 
the  random  path  to  terminal  branches  of  the  smallest  permissible 
size,  whereas  another  might  stop  as  soon  as  he  obtains  a  branch 
that  is  small  enough  to  qualify  as  a  terminal  branch.  Or,  a 
branch  along  a  path  might  be  treated  as  a  terminal  branch  by  one 
enumerator  and  as  path  fruit  by  another .  However ,  for  each  enu¬ 
merator  an  estimate  of  the  total  number  of  apples  on  each  tree  is 
made  using  either  (3.2)  or  (3.3)  as  the  estimator.  The  25  esti¬ 
mates  are  added  together  to  obtain  an  estimate  of  the  total 
number  of  apples  on  the  25  trees.  This  gives  four  estimates , 
one  for  each  enumerator,  of  the  total  number  of  apples  on  the 
25  trees. 

(a)  Assume  that  random  selection  is  performed  correctly  at 
each  stage  of  branching  (after  all  branches  at  the  stage  have  been 
completely  identified  and  measured) ,  and  assume  that  apples  have 


103 


been  correctly  counted.  Do  the  four  estimates  of  the  total 
number  of  apples  all  have  the  same  expected  value  and  the  same 
variance ? 


(b)  Suppose  four  estimates 3  one  for  each  enumera 
the  total  number  of  terminal  branches  on  the  25  trees 
Do  these  estimates  have  the  same  expected  value?  Why? 

(c)  Two  enumerators  measuring  the  csa’s  of  any  g 
of  branches  are  not  likely  to  obtain  exactly  the  same 
values.  Is  this  important?  Discuss. 

(d)  The  assumptions  made  in  (a)  are  subject  to  qu 
Try  listing  some  differences  among  enumerators  that  wi 
not3  have  an  effect  on  the  expected  value  of  an  estima 
total  number  of  apples  on  the  25  trees. 

Exercise  3.11  Suppose3  owing  to  pruning  practice 
many  cases  like  the  following  are  found: 


tor  3  of 
are  made. 

iven  set 
numerical 

es tion . 

II 3  and  will 
te  of  the 

s3  that 


Assume  the  instructions  were  to  always  measure  the  csa  at  the 
base  (point  A)  of  a  branch3  Would  you  expect  the  csa  measurements 
under  the  RP-PPS  method  to  be  ineffective}  or  even  incre as e  the 
sampling  variance 3  compared  with  the  DS-EP  method?  In  cases  like 
the  above  drawing 3  perhaps  measuring  the  csa  at  position  B  would 
be  more  effective .  What  is  your  opinion?  Incidentally  3  this  is 
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a  good  example  of  why  it  is  essential  that  a  research  and 
development  staff  should  have  actual  experience  with  practical 
operations  and  decisions  that  must  be  made  by  enumerators.  Do 
not  expect  high  quality  results  when  instructions  are  not  well 
adapted.  Agreement  between  concepts  (the  theoretical  model) 
and  operations  as  actually  performed  is  of  fundamental  importance . 
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Table  3.1 — Data  by  Branches  for  Apple  Tree  No.  3 


Total  number  of  apples  on  terminal  branches  1850 
Total  number  of  apples  on  path  sections  51 
Grand  total 
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1901 


Table  3.2  Probabilities  of  Selection  and  Estimates  of  the  Total 


Number  of  Apples  on  Tree  No.  3 


Terminal 
branch  no . 

DS- 

EP 

RP-EP 

DS- 

PPS 

RP-PPS 

P1 

Y1 

P2 

Y2 

P3 

Y3 

P 

4 

Y4 

1-1-1 

.03846 

5407 

.03333 

6180 

.06095 

3431 

.05492 

3751 

1-1-2 

.03846 

883 

.03333 

960 

.02206 

1502 

.01988 

1610 

1-2-1-1 

.03846 

1949 

.01667 

4425 

.03366 

2220 

.03103 

2379 

1-2-1-2 

.03846 

3639 

.01667 

8325 

.05276 

2667 

.04864 

2863 

1-2-2 

.03846 

3509 

.03333 

4035 

.04162 

3247 

.03530 

3794 

1-3-1 

.03846 

883 

.03333 

960 

.02206 

1502 

.01998 

1602 

1-3-2 

.03846 

831 

.03333 

900 

.02342 

1332 

.02121 

1414 

1-4 

.03846 

753 

.06667 

405 

.03252 

881 

.02930 

921 

1-5 

.03846 

2339 

.06667 

1320 

.05094 

1776 

.04590 

1917 

2-1-1 

.03846 

1143 

.04167 

1026 

.02092 

2059 

.03073 

1384 

2-1-2 

.03846 

2885 

.04167 

2634 

.04526 

2459 

.06647 

1657 

2-2-1 

.03846 

1975 

.04167 

1794 

.03343 

2265 

.04382 

1706 

2-2-2-1 

.03846 

1507 

.02083 

3090 

.03730 

1552 

.05334 

1222 

2-2-2-2 

.03846 

3067 

.02083 

5972 

.03502 

3363 

.05009 

2489 

2-3 

.03846 

3275 

,08333 

1506 

.04526 

2791 

.05757 

2171 

2-4 

.03846 

2105 

.08333 

966 

.04162 

1949 

.05294 

1509 

3-1-1 

.03846 

831 

.02778 

1101 

.03343 

948 

.02528 

1203 

3-1-2 

.03846 

857 

.02778 

1137 

.02752 

1147 

.02081 

1506 

3-1-3 

.03846 

1117 

.02778 

1497 

.04344 

995 

.03284 

1264 

3-1-4-1 

.03846 

467 

.01389 

2001 

.03343 

530 

.03984 

742 

3-1-4-2 

.03846 

649 

.01389 

2505 

.02615 

931 

.03117 

1078 

3-2-1 

.03846 

961 

.02778 

1263 

.03184 

1150 

.02274 

1542 

3-2-2 

.03846 

1637 

.02778 

2199 

.03229 

1940 

.02306 

2648 

3-2-3 

.03846 

3067 

.02778 

4179 

.04003 

2949 

.02858 

4062 

3-2-4 

.03846 

2339 

.02778 

3171 

.07414 

1238 

.05294 

1665 

3-3 

.03846 

1351 

.11111 

453 

.05890 

900 

.06163 

814 

.99996 

1.00001 

.99997 

1.00001 
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Table  3.3  Variances  of  Estimates  of  the  Total  Number  of  Apples 


on  Each  of  Six  Trees  from  a  Sample  of  One  Terminal  Branch 


No.  of 

terminal 

branches 

csa  of 
trunk 

No.  of 
apples 
on  tree 

Variances 

Tree 

DS-EP 

RP-EP 

DS-PPS 

RP-PPS 

1 

13 

7.0 

214 

(000) 

40 

(boo)" 

28 

Coboy" 

24 

(bool 

22 

2 

27 

20.0 

1448 

882 

1383 

674 

478 

3 

26 

23.0 

1901 

1419 

2815 

755 

800 

4 

20 

16.5 

1658 

1148 

1444 

380 

350 

5 

19 

13.5 

403 

82 

263 

65 

79 

6 

30 

19.5 

1575 

894 

4339 

416 

513 

Total 

135 

99.5 

7199 

4465 

10272 

2314 

2242 
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Figure  3.1  Dot  Chart - Number  of  Apples  vs  Selection  Probabilities  for  RP-EP 
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Figure  3.2  Dot  Chart - Number  of  Apples  vs  Selection  Probabilities  for  RP-PPS 
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TWO-STAGE  SAMPLING 


CHAPTER  IV 


4.1  INTRODUCTION 

Most  sampling  plans  for  estimating  or  forecasting  tree-crop 
production  will  involve  three  or  four  stages  of  sampling.  Typi¬ 
cally,  there  will  be  a  sample  of  orchards  (fields),  a  sample  of 
trees  in  selected  orchards,  and  a  sample  of  branches  from  a 
sample  of  trees.  Fruit  on  the  sample  branches  would  be  counted 
and  a  small  sample  of  fruit  on  the  sample  branches  might  be  se¬ 
lected  for  measurements  of  size  of  fruit. 

This  chapter  illustrates  some  alternative  two-stage  sam¬ 
pling  plans  using  data  for  the  six  apple  trees.  Trees  are  the 
psu’s  (primary  sampling  units)  and  terminal  branches  or  "paths" 
are  the  ssu's  (secondary  sampling  units).  The  six  trees  will 
be  treated  as  a  population  to  be  sampled  and  population  vari¬ 
ance  formulas  will  be  used  to  find  the  first  and  second-stage 
components  of  variance.  Incidentally,  the  problem  of  making 
accurate  counts  of  numbers  of  fruit  on  sample  branches  needs 
serious  consideration.  However,  in  the  illustrations  that 
follow,  attention  is  limited  to  matters  of  sampling. 

In  the  application  of  two-stage  sampling,  psu’s  are  often 

selected  with  probabilities  proportional  to  N^ ,  where  N^  is  the 

t  h 

number  of  ssu's  in  the  i  psu.  For  some  surveys,  sampling  with 
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probability  proportional  to  Isk  has  important  advantages.  When 
the  are  not  known,  approximations  of  N.  are  often  used. 

With  regard  to  sampling  trees,  the  Isk  (number  of  branches 
on  trees)  are  not  known  and  it  is  not  feasible  to  determine  the 
for  trees  in  an  orchard.  Some  other  effective  measure  of  size 
must  be  found  or  the  sample  trees  will  need  to  be  selected  with 
equal  probability.  One  possibility  is  to  use  a  double  sampling 
procedure.  For  example,  a  "large"  sample  of  trees  might  be 
selected  with  equal  probabilities.  For  each  tree  in  the  large 
sample  a  measurement  of  size,  that  takes  relatively  little  time, 
might  be  made  and  used  in  the  selection  of  a  small  sample  of  trees 
from  the  large  sample.  Possible  measures  of  size  are  the  csa  of 
the  trunk,  the  sum  of  the  csa's  of  primary  branches,  and  eye 
estimates  of  the  amount  of  fruit.  The  feasibility  of  double 
sampling  would  depend  upon  the  cost  of  obtaining  the  measurements 
of  size  and  the  relation  between  the  measure  of  size  and  the 
amount  of  fruit  on  the  trees.  Stratification  of  trees  within  an 
orchard  also  needs  to  be  considered.  Sometimes  strata  within  an 
orchard  are  readily  recognized;  for  example,  differences  in  age 
or  variety.  Perhaps  a  relation  between  size  of  trunk  and  number 
of  apples  will  be  found  to  be  effective  only  within  strata  com¬ 
prised  of  trees  of  the  same  variety  and  of  a  uniform  condition. 

Stratification,  systematic  sampling,  or  other  techniques 
might  be  applied  at  any  stage  of  sampling.  However,  for  simpli¬ 
city,  the  discussion  will  be  limited  to:  (1)  simple  random  sampling 
of  psu's  (selection  with  equal  probability  and  without  replacement) 
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and  (2)  sampling  the  psu's  with  pps  (sampling  with  unequal  proba¬ 
bilities  of  selection  and  replacement) .  Within  each  selected  psu 
we  will  assume  that  a  simple  random  sample  of  n^  ssu's  is  selected. 

The  number  of  psu's  in  the  sample  is  m  and  the  number  of  ssu's  in 

m 

the  sample  is  n  =  En^. 

Refer  to  Table  4.1  for  an  exposition  of  the  notation  that 
will  be  used  for  representing  data  for  a  population.  Examine  the 
notation  carefully.  Sample  data  are  represented  in  the  same  way 
except  that  lower  case  letters  are  used. 

Since  a  general  mathematical  formulation  of  estimators  and 
their  variances  is  rather  complex  for  two-stage  sampling,  we  will 
procede  from  specific  cases  to  more  general  description.  The 
primary  purpose  of  the  next  section  is  to  present  an  elementary 
view  of  two-stage  sampling. 

4.2  PRIMARY  SAMPLING  UNITS  EQUAL  IN  SIZE 

The  simplest  case  of  two-stage  sampling  is  one  where  all  psu's 
have  the  same  number  of  ssu's,  where  simple  random  is  applied  at 
both  stages,  and  where  the  same  number  of  ssu's  is  selected  from 
each  psu  in  the  sample.  In  this  case,  and  with  reference  to  the 
notation  in  Table  4.1,  the  N^  all  equal  N  and  the  n^  all  equal  n. 

To  summarize,  the  sampling  plan  under  consideration  is  to  select 
a  simple  random  sample  of  m  psu's  from  a  population  of  M  psu's 
and  a  simple  random  sample  of  n  ssu's  from  each  of  the  m  psu's, 
which  gives  a  total  sample  of  n  =  mn  ssu's. 
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For  illustration  a  hypothetical  population  of  4  psu's 

with  5  ssu's  in  each  is  assumed.  The  20  values  of  Y..  are  pre- 

1 J  F 

sented  in  the  top  part  of  Table  4.2.  Deviations  of  Y^  .  from  Y 
are  also  presented.  In  single-stage  sampling,  there  is  one  com¬ 
ponent  of  variance,  namely  the  variance  of  (Y^j  -  Y)  which  in 
the  illustration  is  487.053. 


In  two-stage  sampling,  each  deviation  (Y^.  -  Y)  divides 
into  two  deviations  as  follows: 


(Y.  .  -  Y)  =  (Y.  -  Y)  +  (Y.  •  -  Y.  ) 
ij  }  l  lj  l 

The  values  of  (Y^  -  Y)  are  a  set  of  deviations  which  reflect 
the  variation  among  psu's  and  the  values  of  (Y^ .  -  Y^)  form  the 
other  set  which  reflects  variation  among  ssu's  within  psu's. 

Turn  to  Table  4.2  and  verify  the  deviations  (components)  (Y^  -  Y) 
and  (Yjj  -  Y^) .  Notice  that  the  between  psu  component,  (Y^  -  Y) , 
varies  from  one  psu  to  another  but  is  constant  within  a  psu.  There 
are  only  M  different  values  of  (Y^  -  Y)  and  selecting  a  sample  of 
m  psu's  is  equivalent  to  selecting  a  sample  of  m  values  of  (Y^  -  Y) 
Also,  study  the  values  of  the  within  psu  component,  (Y-jj  “  Y^)  . 

It  varies  from  one  ssu  to  another  within  a  psu,  but  its  average 
value  is  zero  for  each  psu.  Therefore,  these  deviations  reflect 
only  variation  within  psu's.  The  second  stage  of  sampling  is 
equivalent  to  selecting  mn  of  the  deviations,  (Y^ .  -  Y^) 

Now  consider  the  variance  of  y,  the  mean  of  a  two-stage 
sample.  The  difference  between  y  and  Y  may  be  expressed  as  follows 

7  -  Y  =  h  +  32 
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where  d1  is  the  average  value  of  (Y^  -  Y)  for  the  m  psu's  in  the 
sample  and  d£  is  the  average  value  of  (Y^ .  -  Y^)  for  the  mn  ssu's 
in  the  sample. 

Exercise  4.1  With  reference  to  Table  4.23  suppose  that 
psu's  1  and  3  are  selected  at  the  first  stage  and  that  ssu's  1 
and  4  are  selected  within  psu  No.  1  and  ssu's  3  and  5  are  selected 
within  psu  No.  3.  Find  the  values  of  y,  d^  ,  and  d^ .  Verify 
that  y  -  Y  =  d-^  +  •  Ans.  34-43  =  -9.4  +  0.4. 

Since  is  the  average  of  m  random  values  of  (Y^  -  Y)  and 
is  the  average  of  mn  random  values  of  (Y^  -  Y^) ,  it  follows 
that  d-^  and  d2  are  random  variables.  It  happens  that  d^  and  ^ 
are  independent.  Therefore,  the  variance  of  y  is  equal  to  the 
variance  of  d^  plus  the  variance  of  •  From  knowledge  of  the 
variance  of  the  mean  of  a  simple  random  sample,  one  might  antici¬ 
pate  what  the  variances  of  d^  and  d2  are  and  hence  the  formula 
for  the  variance  of  y  which  is: 


V  (y)  = 


M  -  m 
M 


1  +  N  -  n 


m 


N 


2 

mn 


(4.1) 


2  2 
where  S-^  is  the  variance  of  (Y^  -  Y)  and  S2  is  the  variance  of 

=  2 

the  deviations (Y^ .  -  Y^) .  In  this  case,  S2  is  a  simple  average 

2 

of  the  within  psu  variances,  $2^,  which  is  logical  since  the  psu's 


are  equal  in  size  and  are  selected  with  equal  probabilities.  More 

over,  the  within  psu  sample  size  is  constant. 

2  2 

For  the  illustration,  values  of  S-^  and  S2  as  functions  of 
the  deviations  (Y^  -  Y)  and  (Y^j  -  Y^)  are  shown  at  the  bottom 
of  Table  4.2. 
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In  practice,  the  two  sets  of  deviations  (Y^  -  Y)  and 

=  2  2 
(Y..  -  Y.)>  would  not  be  computed.  The  variances,  S,  and  S?, 

1 J  1  1  Lt 

could  be  calculated  as  follows: 


S 


2 

1 


S 


2 

2 


1 

M(N-l) 


(4.2) 


(4.3) 


Exercise  4.2  Use  Eq.  's  4.2  and  4.3  to  find  the  values  of 

2  2.  —2 

and  S2  in  the  numerical  example.  Explain  why  N  appears  as  a 

divisor  in  Eq.  4.2. 

Exercise  4.3  For  m=2  and  n=2  find  the  variance  of  y  using 
Eq.  4.1.  Ans .118.9 

Exercise  4.4  Show  algebraically ,  that  the  right  hand  side 


of  Eq.  4.3  is  equal  to  — 


l  sl. 

■  2i 


2  . 


M 


3  where  S2^  is  the  variance  among 


.  th 


ssu’s  within  the  i  psu. 

One  partial  check  on  a  variance  formula  is  to  determine 
whether  it  reduces  to  known  formulas  for  special  cases.  Two  special 
cases  are  of  interest:  (1)  When  m  =  M,  two-stage  sampling  becomes 
stratified  random  sampling.  That  is,  the  psu’s  become  strata. 
Observe,  when  m  =  M,  that  the  first  term  on  the  right  side  of  Eq. 
4.1  vanishes  and  the  second  term  becomes  the  variance  for  a  stra¬ 
tified  random  sample  of  n  units  from  each  stratum  (psu) .  (2)  When 

n  =  N,  two-stage  sampling  reduces  to  single-stage  cluster  sampling. 
In  this  case  the  last  term  in  Eq.  4.1  vanishes,  leaving  the  first 
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term  which  is  the  variance  for  a  cluster  sample  where  the  clus¬ 
ters  (sampling  units)  are  the  psu's. 

Exercise  4.5  Suppose  m=l  and  n=l.  In  this  case  the  selec¬ 
tion  of  one  psu  at  random  and  the  selection  of  one  ssu  within 
it  is  equivalent  to  a  single-stage  sample  of  one  ssu.  There  fore  3 
the  variance  of  y  given  by  Eq.  4.1  when  m=l  and  n=l  should  be 
equal  to  the  variance  of  y  for  a  single-stage  random  sample 
when  n=l.  Verify  this  using  the  data  in  Table  4.2.  Remember 
the  appropriate  variance  formula  for  the  single-stage  sample  is 
Eq.  1.4. 

It  is  important  to  study  the  structure  of  the  variance 
formula,  Eq .  4.1,  for  the  variance  of  y.  When  the  number  of 
psu's  in  the  sample  is  fixed,  increasing  the  size  of  the  sample 
in  each  psu  reduces  only  the  second  component  of  variance.  As 
n  increases,  a  point  is  reached  where  the  among-psu  variance  is 
the  major  component  and  further  increases  in  n  contributes  very 
little  to  reducing  the  variance  of  y.  Notice  that  increasing  m 
reduces  both  components  when  n^  is  constant  for  all  psu's. 

4.3  PRIMARY  SAMPLING  UNITS  UNEQUAL  IN  SIZE 

Populations  having  psu's  with  equal  numbers  of  ssu's  are 
relatively  infrequent.  In  this  section,  it  is  assumed  that  the 
numbers,  N^,  of  ssu's  vary  and  that  simple  random  sampling  (with¬ 
out  replacement)  is  applied  at  both  stages. 

As  discussed  in  Chapter  I,  Sec.  1.1.2,  "P"  or  "p"  with 
appropriate  subscripts  refer  to  selection  probabilities  on  the 
occasion  of  a  particular  random  draw  and  "f"  with  an  appropriate 
subscript  refers  to  the  probability  that  a  particular  unit  has 
of  being  in  the  sample. 
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A  general  expression  for  the  probability,  f ^ ,  which  any 
given  ssu  has  of  being  included  in  a  two-stage  sample  is: 


.  =  f.f  (j  |  i) 
ij  i  1  J 


(4.4) 


where 


.  th 


f ^  is  the  probability  which  the  i  psu  has  of  being 


in  the  sample,  and 


.  th 


f(i|j)  is  the  conditional  probability  which  the  j  ssu 

t  h 

in  the  i  psu  has  of  being  in  the  sample,  given 
t  h 

that  the  l  psu  is  in  the  sample  of  psu's. 

m  ni 

With  simple  random  sampling  at  both  stages,  f^  =  ,  and  f(j|i)  =  — 


Since  f.  is  constant  for  the  case  under  consideration,  let  f.  = 
i  i 

f^  which  is  the  sampling  fraction  at  the  first  stage.  Also,  let 

f(j\i)  =  f2^  which  is  the  sampling  fraction  at  the  second  stage 


t  h 

within  the  i  psu.  Then  Eq .  4.4  reduces  to: 


f .  .  =  f ,  f  9  • 

ij  1  2i 


(4.5) 


If  the  f2^  (the  sampling  fractions  at  the  second  stage)  are  con¬ 
stant,  f^  is  constant  and  every  ssu  in  the  population  has  the 
same  chance  of  being  in  the  sample.  Then,  Eq .  4.5  becomes: 

f  =  f  f 
12 

where  f2  is  the  constant  second-stage  sampling  fraction.  However, 
in  the  interest  of  generality,  a  requirement  that  f2^  be  constant 
will  not  be  specified  at  this  point  in  the  discussion. 

An  estimator  of  the  population  mean,  Y,  is 


y  - 1  sr)(M)  i 


N .  y . 

i7  i 
m 


(4.6) 
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where 


n . 
z1 


y. 


=  1 


yii 


n 


-  th 


is  the  average  of  ru  ssu's  in  the 


sample  from  the  i  psu  in  the  sample.  Study  the  estimator 
4.7  and  observe  that; 

—  i 

is  an  estimate  of  Y-,  the  total  for  the  i  psu; 


m  N .  y . 

2_ J  i 

Z  — - —  is  an  average  of  the  estimated  totals  for  the  m 
i 

psu's  in  the  sample;  therefore. 


m  N .  y  . 

(M) Z  — —  is  an  estimate  of  the  population  total  and  (^) 

in  Eq.  5.6,  changes  the  estimated  total  to  an 
estimate  of  Y. 


The  variance  of  y  is  given  by: 


V (y)  =  - 
w  J  m 


— 

M 

N2sf  | 

(i-fp5  ♦  1 

£  M(l-f„) 

i  2i  | 

n . 

N 

i 

1  J 

(4.7) 


M  -  2 
£(Y,  -  Y)Z 

,  c2  1  i  1 
where  S,  =  — 

1  NZ 


M 


j —  is  the  variance  among  psu  totals  divided 


—2  2 

by  N  so  will  be  expressed  on  the  basis  of  one  ssu,  and 
N. 

Z1  =2 

2  i  (Yii  '  V 

S2^  =  - — 3 — j—  is  the  variance  among  ssu's  within  the 

i 

.  th 

i  psu . 

12 

The  first  part  of  4.7,  —  (l-f-^)S^,  is  the  variance  of  y  assuming 
all  of  the  m  psu's  are  enumerated  completely.  That  is,  the 
theory  for  single-stage  sampling  applies  to  the  first  stage. 
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The  quantity: 


t1  '  f2i> 


2  2 
NiS2i 

ni 


th 


in  Eq.  4.7  is  recognizable  as  the  variance  of  N^y^  where  y^ 
is  the  mean  of  a  simple  random  sample  of  n^  ssu's  in  the  i 
psu . 


Eq.  4.7  was  written  in  the  above  form  for  comparison  with 
other  variance  formulas  given  later  for  two-stage  sampling.  The 
second  term  within  [  ]  could  be  written  as  follows: 


(4.8) 


because 


M 

¥ 


¥  *• 


Expression  4.8  shows  that  the  variances  of 


Nf^i  are  summed  over  all  psu's  in  the  population  and  the  sum  is 
divided  by  M  giving  an  average  of  such  variances.  The  variances 
of  N^y^  receive  equal  weight  in  the  average  because  the  psu's 
are  selected  with  equal  probabilities.  Since  the  average  variance 

=  _  O 

°f  N^y^  pertains  to  psu  totals,  the  divisor  N  appears  in  4.8  to 
convert  the  variance  to  a  basis  of  one  ssu.  Such  analysis  of 
a  formula  is  helpful  in  determining  whether  one  has  the  right 
formula  for  a  particular  purpose. 


Exercise  4.6  If  the  variance  formulas  (4.1)  and  (4.7) 
are  correct ,  formula  (4.7)  should  reduce  to  (4.1)  when  =  N 
and  n^  =  n.  Show  that  this  is  true. 
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are  constant 


n . 

When  the  second-stage  sampling  fractions 

i 

and  equal  to  f^,  the  estimator,  (4.6),  reduces  to: 

z  Ey . 


y  = 


f  ^mN 


(4.9) 


and  its  variance,  (4.7),  reduces  to: 


y(y)  =  a  -  q)  -f-  *  a 


f2) 

mn 


(4.10) 


,2  . 


where  is  the  same  as  in  4.7, 


7  M  N.  ? 

S9  =  E  -4-  S9. 

2  ^  N  2 1 , 


and 


n  = 


M 

E  n . 
i  1 


Ef?N. 

t  -1  =  f9N 

M  2 


Exercise  4.7.  Show  that  Eq.  's  (4.9)  and  (4.10)  follow 

ni 

from  (4.6)  and  (4.7)  when  f7  =  — — 

Exercise  4,8  Show  that  f2mN,  in  Eq.  4.9 3  is  equal  to 

the  expected  sample  size.  That  is  3  show  that  E(n)  =  f2mN  where 
m 

n  =  E  n . .  In  practice  one  would  probably  use  n,  the  actual 

i 

sample  size 3  in  the  estimator  instead  of  the  expected  size3 
f2inN.  Moreover  3  N  is  not  known  in  most  practical  applications . 
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4.3.1  NUMERICAL  EXAMPLE 


As  a  numerical  example,  the  apple  tree  data  presented  in 
Table  2.1  will  be  treated  as  a  population  to  be  sampled.  The 
psu’s  are  trees  and  ssu's  are  terminal  branches.  The  number 
of  trees  in  an  orchard  is  usually  large  and  in  practice  the 
number  of  sample  trees  selected  from  an  orchard  would  be  rela¬ 
tively  small,  that  is  (1  -  f^)  would  be  nearly  equal  to  1. 
Accordingly,  for  this  illustration,  (1  -  f^)  is  assumed  to  be 
1  even  though  M  =  6  and  (1  -  f-^)  =  pp  is  considerably  less 
than  1 . 


Suppose  we  are  interested  in  knowing  what  the  sampling 
variance  is  for  the  following  three  allocations  of  a  sample  of 
four  terminal  branches  assuming  simple  random  sampling  at  both 
stages : 


Allocation 


No.  of  Trees 
m 


No.  of  Branches 
Selected  from 
Each  Tree 

n .  =  n 

l 


11  4 
2  2  2 
3  4  1 


To  find  the  variances  for  the  three  allocations  we  need 
part  of  the  results  in  Table  2.6.  The  relevant  results,  N.  Y. , 

1  9  1 

2 

and  S^,  from  Table  2.6  are  included  in  Table  4.3  along  with 
some  other  information  that  will  be  used  later. 
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In  each  allocation,  n^  is  constant  (the  same  for  all 

ni 

trees)  which  means  that  is  not  constant  and  the  branches 

1 

do  not  have  equal  probability  of  being  in  the  sample.  Thus,  the 
estimator,  Eq .  4.6,  and  its  variance,  Eq.  4.7,  are  applicable. 
The  variances  for  the  three  allocations  are  presented  in  Table 
4.4. 


Exercise  4.9  Refer  to  the  data  presented  in  Table  4.3 3 

2 

columns  ,  and  S2^  and  perform  the  calculations  that  are 

needed  to  obtain  the  results  presented  in  Table  4.4  for  m  =  2 
and  n^  =  n  =  2.  Assume  that  f^  is  negligable . 

Exercise  4.10  Complete  the  following  table: 


m 


n 


1 

2 

4 

2 

4 

8 

4 

8 

16 


V  (y ) 


1167.0 


Variance  Components 
Among  vsu  rs  Within  vsu  ' s 


306.5 


If  you  understand  the  variance  formula  4.7  and  the  results  in 
Table  4.33  this  table  can  be  completed  very  easily.  Firsts  fill 
in  the  "Among  psu's "  column  by  copying  the  appropriate  numbers 
from  Table  4.4.  Consider  how  to  fill  in  the  "Within  psu's"  column 
by  making  simple  changes  in  the  within  psu  components  in  Table 
4.4.  Study  the  results.  For  a  constant  value  of  n  and  an  increase 
in  m  from  1  to  4  there  is  a  75  percent  reduction  in  the  variance 
of  y;  but s  for  a  constant  m}  increasing  n  from  1  to  4  reduces  the 

variance  of  y  by  less  than  50  percent . 
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Exercise  4.11  One  of  the  numbers  in  Table  4.4  is  the 


sampling  variance  for  m  =  2  and  n.  =  N..  What  is  the  number? 

Is  Is 


Exercise  4.12  Find  the  probability  that  any  given  terminal 
branch  on  tree  No.  1  has  of  being  in  the  sample  when  m  =  2  and 
n.  =  2  for  all  trees.  What  is  the  probability  for  tree  No.  3? 

Is  the  unequal  probability  something  to  be  concerned  about? 

In  what  ways? 

It  is  of  interest  to  compare  the  variance  for  a  simple 
random  (single-stage)  sample  of  4  branches  with  the  variances 

A 

of  y  in  Table  4.4.  The  variance  among  the  135  branches  is  1,762 

(see  Table  2.6).  Hence,  the  variance  of  the  mean  of  a  sample 

1762 

of  4  branches  is  — ^ —  =  440  ,  disregarding  the  fpc.  The  answer, 

A 

440,  is  less  than  the  variances  of  y  in  Table  4.4.  This  is 
expected  with  the  possible  exception  of  the  allocation  m  =  4 
and  n.  =1,  which  has  a  variance  equal  to  583.5.  However,  when 


one  recognizes  in  the  specified  two-stage  plans  that  all  branches 
do  not  have  the  same  probabilities  of  selection,  it  is  resonable 
to  expect  that  the  answer  for  simple  random  sampling  would  be 
less  than  583.5. 

Suppose  we  wish  to  give  every  branch  an  equal  chance  of 
being  in  the  sample.  Considering  samples  of  4  branches  the  over¬ 
all  sampling  fraction  would  be  yys’  If  we  specify  that  m  =  2, 


then  f-^  =  y  and  all 


4 

(or  f 2 )  should  equal  Since  the 


are  small  and  the  n^  must  be  integers,  it  is  not  possible  to  have 
n '  4 

all  exactly  equal  to  This  presents  a  type  of  practical 

i 

problem  that  often  occurs  when  working  with  small  integers.  Ways 
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of  dealing  with  this  problem  will  not  be  discussed  at  this 

ni 

point.  Instead,  we  will  procede  as  though  the  fraction 

1 

is  sufficiently  close  to  ^  to  warrant  use  of  the  unweighted 

average  of  the  sample  data  as  the  estimator  and  the  variance 

formula  4.10.  Assuming  (1  -  f2)  =  1>  for  reasons  explained 

2  2 

above,  and  substituting  the  numerical  values  of  and  in 
4.10,  we  have: 

V(y)  ■  jj  [917.1  +  (1  -  f2)  i^-1  (4.11) 

4 

When  m  =  2  and  f2  =  j^-,  the  value  of  n  is  2  and  the  variance  of 

A 

y  is  769.9.  This  answer  compares  with  797.6  in  Table  4.4. 

Exercise  4.13  Verify  the  numbers ,  917.1  and  1367 3  in  Eq. 

4.  11. 

It  is  often  desirable  to  specify  that  all  ssu's  in  the 

population  have  an  equal  chance  of  being  in  the  sample.  As 

discussed  above,  one  way  of  fulfilling  this  requirement  is  to 

select  psu's  with  equal  probability  and  apply  a  constant  sampling 

fraction  at  the  second  stage  of  sampling.  But,  when  the  sizes 

of  the  pus's  vary  widely,  this  method  often  has  two  important 

disadvantages:  (1)  Variance  associated  with  variation  in  the 

sizes  of  the  psu's  is  included  in  the  variance  of  an  estimate 

2 

unless  such  variation  is  reduced  by  design.  Notice  that  in 
4.7  is  the  variance  among  psu  totals  rather  than  the  variance 
among  psu  means.  Incidentally,  an  auxiliary  variable(s)  might 
be  useful  in  reducing  the  sampling  variance  associated  with  the 
first  stage  of  sampling.  (2)  When  the  second-stage  sampling 


\2S 


fraction,  ,  is  constant,  n^  is  proportional  to  and  the 
1 

workload  varies  from  one  psu  to  another.  For  many  surveys,  it 

n-: 

is  important  for  reasons  of  economy  that  n^,  rather  than  be 

i 

constant.  Selecting  psu's  with  pps  is  often  very  helpful  in 

overcoming  these  disadvantages. 

Exercise  4.14  Under  the  plan  of  applying  a  sampling  frac- 
4 

tion  of  to  each  tree  that  is  selected 3  suppose  that  trees 

numbered  1  and  3  are  selected .  Find  the  values  of  n.  for  these 

4 

two  trees  where  n^  is  Nk  rounded  to  the  nearest  integer .  Also3 
find  n=  En^.  Do  the  same  assuming  trees  numbered  2  and  4  are 
selected .  This  illustrates  that  the  size  of  the  sample ,  n  =  £n^, 

ni 

is  a  random  variable.  Also 3  in  this  case3  cannot  be  exactly 

l 

constant .  One  should  consider  whether  there  is  an  appreciable 
bias  in  the  estimator  (4.9).  Use  (4.6)  instead  of  (4.9)  unless 
there  is  assurance  that  any  bias  in  (4.9) 3  owing  to  unequal  pro¬ 
babilities  of  the  ssu's  being  in  the  sample3  is  negligible. 

4.4  SELECTION  OF  PSU’S  WITH  PPS 

Consider  a  sample  of  m  psu's  selected  with  replacement  and 

with  selection  probabilities  P-^,  •••>  (See  section  1.1.2  in 

Chapter  I).  Let  n^  be  the  size  of  a  simple  random  sample  of  ssu's 

t  h 

that  is  to  be  selected  from  the  i  psu  in  the  event  that  it  is 

t  h 

selected.  If,  by  chance,  the  i  psu  is  selected  a  second  time 
another  sample  of  n.  ssu's  is  selected.  For  a  sample  of  n  psu's 


the  estimator  is: 


1  1  m  Ni>"i 

y  E  -V 

l 


(4.12) 
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Remember  to  interpret  "iM  as  an  index  of  the  psu's  selected 

by  the  m  random  draws.  Notice  that  N.y.  is  an  estimate  of  a 

=  1  1 

N.y. 

psu  total  and  that  — 1  1  is  an  estimate  of  the  population 


total,  Y,  based  on  a  sample  of  one  psu  and  a  simple  random 
sample  of  n^  ssu's  within  it.  Thus,  there  are  m  estimates  of 


1  m  Ni^i 

Y  and  (— )  z  -  is  an  average  of  these  estimates.  The  factor 

m  7  p- 


l 


^  makes  y  an  estimator  of  Y.  The  variance  of  y,  in  Eq.  4.12,  is 


v(y)  -  s 


2  1  “  1  ,  ,  ,  N?Sj. 

1  N  i  i  ^  n. 


(4.13) 


where 


i  M  Y. 

?  pi(pi  '  Y) 
N  l  l 


and 


N. 

l 


-  J 
2i 


J  CYij  -  V 


N.  -  1 

l 


Exercise  4.15  Compare 


m 


in  4.13  with  the  variance  of 


y^  in  Table  1.1 3  using  the  alternative  expression  for  a ^  in  the 
variance  of  .  Change  the  notation  used  in  Chapter  1  to  conform 
to  the  notation  used  for  psu’s.  This  gives; 

,  l  1  m  Yi  ? 

v(y4)  ■  £)(^)  f  pitP7  -  Y> 

Mi  i 
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Why  is  this  expression  for  V(y^)  different  from  the  between 
psu  part  of  the  variance  in  4.13?  In  terms  of  the  notation 

A  _ 

for  two-stage  sampling  is  an  estimate  of  Y  rather  than  Y. 
Change  y^  so  it  will  be  an  estimator  of  Y  and  make  the  cor¬ 
responding  change  in  V(y^).  Your  answer  should  agree  exactly 
a  2 

with  -  in  (4.13). 

m 

2 

Notice  the  correspondence  between  S-^  in  Eq .  4.7  and 

the  variance  of  y^ ,  plan  1,  in  Chapter  I;  also,  notice  the 

2 

correspondence  between  in  Eq.  4.13  and  the  variance  of 

y^  ,  plan  4,  in  Chapter  I.  The  discussion  in  Chapter  I  of  the 
efficiency  of  plan  4  compared  to  plan  1  is  relevant  to  the 

2 

first  stage  of  sampling.  If  is  a  good  measure  of  size, 

2 

will  be  considerably  less  than  S^. 

Compare  the  components  of  variance  in  Eq.  4.7  and  Eq.  4.13 
which  pertain  to  the  second  stage  of  sampling.  The  only  differ¬ 
ence  is  a  reflection  of  the  difference  in  the  probabilities  of 
selection  at  the  first  stage.  When  the  probabilities  are  equal, 

P.  =  i  and  substituting  ^  for  P.  in  4.13  gives  4.7. 

l  M  6  M  l  6 

In  Eq .  4.4,  f^  was  expressed  as  the  probability  that 

any  given  ssu  has  of  being  in  a  sample  assuming  the  sample  at 

both  stages  was  simple  random  sampling  without  replacement. 

This  equation  now  needs  modification  to  be  in  accord  with  sampling 

at  the  first  stage  with  unequal  probability  and  with  replacement. 

An  appropriate  probability  equation  is: 


P.f9. 

l  2i 


(4.14) 
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where  is  the  selection  probability,  at  any  given  random 

t  h 

draw,  for  the  i  psu  in  the  population, 

f 2^  as  defined  before,  is  the  sampling  fraction  within 
t  h 

the  i  psu  of  the  population,  and 

f—  is  the  probability  which  the  j  ^  ssu  in  the  i*'*1  psu 

of  the  population  has  of  being  in  a  sample  obtained  by 

selecting  one  psu  with  pps  and  selecting  a  simple 

random  sample  of  n^  ssu's  within  the  selected  psu. 

It  is  in  the  context  of  the  probability  Eq.  4.14  that  the 

estimator,  4.12  and  its  variance,  4.13,  are  applicable,  assuming 

m  independent  random  selections  of  psu's. 

The  estimator,  Eq.  4.12,  and  its  variance,  Eq.  4.13, 

are  for  any  given  set  of  selection  probabilities  at  the  first 

stage  and  any  given  set  of  sample  sizes,  n^,  at  the  second 

stage.  An  important  special  case  exists  when  f!j  ,  in  Eq .  4.14 

is  held  constant  and  when  the  psu’s  are  selected  with  probabilities 

N. 

proportional  to  ,  that  is,  when  By  letting  f'  be  the 

constant  value  of  f!^  ,  we  obtain  the  following  results  from  Eq . 
4.14: 

n  .  =  f  ’  N  =  n 

l 


and 


f  =  2- 
r2i  N. 


That  is,  the  sample  size  within  a  psu  is  constant,  and,  since 
fj .  is  also  constant,  the  sample  is  self -weighted . 
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The  estimator  and  its  variance  become: 


y  = 


EEy .  . 
n 


(4.15) 


and 


v(f)  =  lL2  +  T  g  2  N.(l  - 

L  n  i 


£2i^S2i 


M 


where  ^  E  N.(Y.  -  Y)  ^ 

1  N  ^  l v  l  J 


(4.16) 


For  computational  purposes  one  might  use: 


and 


,  M  Y. 

—  E  — — 
N  .  N. 
l  i 


(I)2 


M 


Je  N.(l  -  f,.)sL  *  J  EN.  S?. 
N  ^  r  2iJ  2i  N  l  2i 


—  ES^ 

N  lb2i 


Exercise  4.16  Show  that  Eqs.  4.12  and  4.13  reduce 

N. 

to  4.15  and  4.16  when  P.  =  ^  and  n.  =  n. 

l  N  l 

When  the  are  not  known,  estimates  of  or  a  suitable 

measure  of  size  might  be  used  in  place  of  isL  .  In  this  case, 

assuming  f!^  =  f',  the  sampler  would  choose  a  value  of  f'  such 

that  f'N  is  the  desired  average  size  of  sample  from  a  psu.  Since 

the  selection  probabilities  for  psu’s  are  known,  the  second-stage 

f  ’ 

sampling  fraction  f 2 ^  =  p—  would  be  calculated  for  each  selected 

r  i 

psu.  Application  of  these  second-stage  sampling  fractions  gives 
a  self -weighted  sample.  The  n^  will  be  nearly  equal  if  the 
measure  of  size  is  close  to  being  proportional  to  N^.  The  esti¬ 
mator,  Eq .  4.12,  and  its  variance,  Eq .  4.13,  are  applicable. 

They  could  be  modified  by  making  use  of  the  fact  that 
n. 

p.£,.  =  p. 

2i  q  N.  r  • 
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4.4.1  NUMERICAL  EXAMPLE 

Exercise  4.17  With  reference  to  the  apple  tree  example 3 
we  found  for  simple  random  sampling  at  both  stages  that  the 
sampling  variance  was  797.6  when  m=2  and  ru=n=2  ( See  Table  4.4). 
For  comparative  purposes find  the  sampling  variance  for  m=2 
and  n=2  when  the  trees  are  selected  with  probabilities  propor¬ 
tional  to  N^.  The  data  needed,  are  found  in  Table  4.33  columns 

2  2  12 
headed  N^,  Y^,  and  Find  the  values  of  0-^  ,  j^-zN^S2^,  an ^ 

12  - 

NjS2i’  then  compute  the  variance  of  y  for  m=2  and  n=2.  Ans . 
532.6. 

Substituting  results  from  exercise  4.16  in  Eq .  4.16  gives: 

V(y)  =  439-7  +  1367  -  n  (57 . 99  (417 

m  mn 

439  74 

For  m=2  the  between  psu  variance,  — j1 =  219.9,  compares  with 

458.6  (see  Table  4.4)  when  two  psu's  are  selected  with  equal 
probability.  As  indicated  by  this  result,  selecting  psu's  with 
pps  is  often  very  important  in  reducing  the  between  psu  compo¬ 
nent  of  variance.  For  m*2  and  n=2  the  within  psu  component  in 
Eq .  4.17  is  equal  to  312.8  which  compares  with  two  other  results 

that  were  obtained  when  the  trees  are  selected  with  equal  proba- 

n . 

bilities:  339.0  when  n.  =  2,  and  311.4  when  — ^  is  constant  and 

1 

n=2.  The  first  result,  339.0,  was  recorded  in  Table  4.4  and 
the  second,  311.4,  is  readily  obtained  by  Eq.  4.11. 
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N. 

Suppose  that  one  tree  is  selected  with  probability  ~ 
and  that  one  branch  is  selected  from  it  with  equal  probability. 
In  this  case,  m=l,  and  n- 1 ,  and  the  variance  of  y  according  to 
variance  formula  4.17  is  1748.8.  The  probability  of  selecting 


Ni  1 

any  given  branch  in  the  population  is  ( ^— )  ( )  .  This  is  a 


■N- 


special  case  of  two-stage  sampling  that  is  the  same  as  a  single 
stage,  simple  random  sample  of  one  branch.  We  found  earlier 
that  the  variance  among  the  135  branches  was  1762.  The  exact 
variance  for  a  simple  random  sample  of  one  branch  is: 


1261  .  1748 . 8 

4.5  UNEQUAL  PROBABILITY  OF  SELECTION  AT  BOTH  STAGES 


As  a  further  exposition  of  the  theory  for  two-stage 
sampling,  suppose  a  sample  of  trees  is  selected  with  replace¬ 
ment  and  with  selection  probabilities  proportional  to  trunk 
size.  Also,  suppose  that  the  method  of  sampling  at  the  second 
stage  is  the  random-path  method,  RP-PPS,  that  was  discussed 
in  Chapter  III.  You  may  recall  that  the  random-path  method 
was  presented  in  the  context  of  sampling  with  replacement. 

When  the  sampling  at  both  stages  is  with  unequal  proba¬ 
bility,  the  estimator  of  the  population  total  Y  is: 

r  i 


y. 


m  n. 
z  z1 
i  j 


y 


ii 


p. . 

LAll 


n 


(4 


18) 


where 


m 

n  =  Zn. 

.  i 

i 


Pij  '  PiPOlU 
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1"  Vi 

is  the  selection  probability  for  the  1  psu 
in  the  sample,  and 

p(j|i)  is  the  selection  probability  for  the  j  ^ 
ssu  given  that  its  psu  has  been  selected. 


yii 

Consider  the  quantity  — in  the  estimator.  When  the  value 

Pij 

for  a  unit  in  the  sample  (in  this  case,  y^.)  is  divided  by 
its  selection  probability  (in  this  case,  Pjj)  the  quotient  is 
an  estimate  of  the  population  total.  Therefore,  y^  in  Eq .  4.18 
is  an  average  of  n  estimates  of  Y,  one  estimate  from  each  branch 
in  the  sample. 

The  subscript  "t"  was  added  to  y  because  it  is  an  esti¬ 
mator  of  Y  rather  than  Y.  Notice  that  the  estimator  does  not 
contain  N.  In  practice,  one  finds  many  populations  to  be 
sampled  where  N  is  unknown.  An  estimate,  N,  of  N  might  be 

made  from  a  sample  and,  if  needed,  X-  could  be  used  as  an  estimate 

N 

of  Y.  An  estimator  of  N  is  obtained  by  substituting  "1"  for 
y . .  in  4.18. 


Exercise  4.18  Suppose 3  for  m=3,  and  n^=2,  that  appli¬ 
cation  of  the  above  method  to  the  apple  tree  population  gives 
the  following  sample : 

Population  index  Sample  index 


values 

of  i  and  j 

values 

of  i  and  j 

Pi 

P(j  1  i) 

7ij 

Tree 

Path 

Tree 

Path 

1 

2-2-3 

3 

1 

0.07035 

.15996 

59.5 

1 

4-1 

3 

2 

0.07035 

.07779 

7.0 

3 

1-2-1-2 

1 

1 

0.2312 

. 04864 

139.3 

3 

2-3 

1 

2 

0.2312 

.05757 

125.0 

3 

3-1-2 

2 

1 

0 . 2312 

.02081 

31.3 

3 

1  -  2  - 1  -  2 

2 

2 

0.2312 

.04864 

139.3 
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Tree  No.  3  and  path  1-2-1-2  were  selected  twice.  The  selection 
probabilities  were  proportional  to  the  trunk  sizes  which 

are  presented  in  Table  4.3.  Verify  the  values  of  p^.  For  tree 
No.  3  in  the  population,  the  conditional  probabilities ,  P(j  ID, 
are  the  probabilities  in  Column  of  Table  3.2.  Thus,  the  above 
values  of  P(j  U)  and  for  the  branches  in  the  sample  from  this 

tree  were  taken  from  Table  3.2.  The  values  of  P  ( j  1  i)  and  for 

tree  No.  1  are  from  records  not  reproduced  herein.  Using  Eq.  4.18 
as  the  estimator ,  calculate  the  estimate  of  the  total  number  of 
apples.  The  answer  is  7873.,  which  is  an  estimate  of  7199.,  the 
total  number  of  apples  including  " path "  apples  (See  Table  4.3). 


To  find  the  variance  of  y  ,  refer  to  Eq.  4.13  and  make 


two  modifications: 


2  2  2 

(1)  for  the  first  stage  we  want  N  instead  of 

because  y  is  an  estimator  of  Y  rather  than  Y,  and 

(2)  for  the  second  stage,  the  part  of  the  formula  repre¬ 
senting  the  variance  of  an  estimate  of  Y^  for  a  simple 

random  sample  of  n.  needs  to  be  changed.  That  is, 

2  2  1 

NT  • 

(1  -  f £ ^ — —  needs  to  be  replaced  by  the  cor- 


n . 
1 


responding  variance  for  sampling  within  the  i 
with  pps.  Also,  ^  needs  to  be  dropped. 


.  th 


=  i 


N 


M  Y.  ?  M  ,  q2 

Z  P  -  Y)Z  +  Z(±-)  1_ 

i  i  i  i  n 


r  1 


psu 
This  gives: 

(4.19) 
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where 


N. 

1 


Y.  . 


S  .  =  E  P(j  |i)  -  Y.) 

n  j  J  1  P  C  J  1 1 )  1 


The  subscript  r  signifies  random  path 


For  Tree  No.  3  the  values  of  P(j  | i)  and  the  values  of 

Y.  . 

tw*?-i  are  in  columns  P,  and  Y.  of  Table  3.2.  From  these  two 

P(j  I  i)  4  4 

2 

columns  the  value  of  S  ^  for  tree  number  3  can  be  computed. 

2 

The  answer,  800,000  is  recorded  along  with  other  values  of  Sr^ 
in  the  last  column  of  Table  4.3. 

Exercise  4.19  When  n^=n,  the  second  term  in  []  of  Eq. 

,  M  S2.  X. 

r  i 

4.19  becomes  —  Z  -  -f.-  .  In  the  problem  under  consideration ,  P .  =-y- 
n  i  i 

where  is  trunk  size.  From  the  data  in  Table  4.3}  find  the 
2 

S  .  Y.  9 

value  of  E  p-  and  of  E P ^  C p —  '  Y)  .  When  your  results  are  sub- 
r  i  r  i 

stituted  in  4.19 ,  you  should  have: 

V(y  )  -  -  15,322,000  +  11 ,46z  »000  j 
m  l.  n  j 

Exercise  4.20  From  the  sample  data  given  in  Exercise 

4.16  estimate  the  total  number  of  terminal  branches  on  the  six 

trees.  Ans .  122.0. 


When  the  equation  in  Exercise  4.19  is  divided  by  N  we 


obtain : 


■i  1‘ 


v(y)  =  £  292  + 


629 


n 


1 35 


To  summarize,  the  following  variance  equations  have  been 
obtained  for  three  alternative  two-stage  plans  for  sampling 
the  small  population  of  apple  trees: 


(1) 


V (y)  =  - 
7  m 


917.1  + 


(l-f2)  1367 


n 


J 


n . 

for  simple  random  sampling  at  both  stages,  where  — is  constant 
and  equal  to  and  1  -  f ,  was  assumed  to  be  equal  to  1, 


(2) 


V (y)  =  - 
w  7  m 


439.7  + 


1367  - 


n  ( 5  8 . 0  )~j 

n 


for  sampling  trees  with  probability  proportional  to  and  a 
simple  random  sample  of  n  branches  from  each  selected  tree, 
and 


(3) 


V(y)  =  - 
w  7  m 


292  + 


629 


n 


for  sampling  trees  with  probability  proportional  to  (trunk 
size)  and  application  of  the  RP-PPS  method  to  the  sample  trees. 

The  results  are  too  limited  to  provide  a  basis  for  genera¬ 
lization  . 
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1/ 

Table  4.1  Representation  of  Population  Data  for  Two  Stage  Sampling 


psu 

ssu 

psu 

total 

psu 

mean 

variances 

1  . . . 

j  ... 

Nf 

1 

Y  . 

•  Y  .  • 

Y 

Y 

s2  - 

?  ry 

r  u 

^)2 

X11 

lj 

Y1 

T1 

21 

Ni 

-  1 

Y 

•  Y  •  • 

y 

Y 

Y 

s2  = 

•(Y-.  - 
J  ^ 

Y^2 

ll 

iNi 

i 

i 

52i 

N. 

l 

-  1 

M 

Y 

.Y  .  . 

.Y 

Y 

Y 

S2  = 

?  (Y 
r  Mj 

Y  )  2 

MJ 

*M1 

mnm 

M 

M 

b2M 

nm 

-  1 

1/ 


A  single  bar  is  used  for  an  average  of  psu  totals  and  a  double 

bar  "="  indicates  an  average  of  secondary  units.  A  subscript  1  or 
2 

2  affixed  to  S  indicates  first  or  second  stage  variance.  See  defi¬ 
nitions  below. 

"t  h 

Y. .  is  the  value  of  the  characteristic  Y  for  the  j  ssu  in  the  l 
1 1  J 

J  psu, 


N. 

l 


Y^  =  2  Y^  is  the  total  of  Y  for  the  i 


.  th 


psu; 


M  Ni  M 

Y  =  T,  l  Y..  =  IY.  is  the  total  of  Y  for  the  population, 
i  j  1-!  i  1 

M  is  the  number  of  psu's  in  the  population, 

f'L  is  the  population  number  of  ssu's  in  the  i*"*1  psu, 

M 

N  =  z  N.  is  the  number  of  ssu's  in  the  population, 
i 


Y  =  jjj  is  the  population  mean  per  psu, 

=  Y 

Y  =  jq-  is  the  population  mean  per  ssu, 

=  Y'  th 

Y^  =  is  the  average  value  of  Y  per  ssu  in  the  l  psu, 

N  =  j^j  is  the  average  number  of  ssu's  per  psu, 


,  N.  (Y.  -  -  Y.)  t, 

sf.  =  I1  — - i —  is  the  variance  among  ssu's  in  the  l  psu,  and 

Zi  .  Ni 

M 


S1  =  - - MT 


£  (Y-  -  Y)2 


is  the  variance  among  psu's  on  the  basis  of 
one  ssu. 
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Table  4.2  Components  of  Variation  for  a  Hypothetical  Population 


psu 

1 

'  2 

3 

4 

1 

67 

Values 

45 

of  Y. . 

51 

20 

2 

32 

27 

82 

39 

3 

14 

25 

21 

30 

4 

55 

48 

72 

63 

Values 

of  (Y. . 
il 

-  Y) 

1 

24 

2 

8 

-23 

2 

-11 

-16 

39 

-4 

3 

-29 

-18 

-22 

-13 

4 

12 

5 

29 

20 

Values 

of  (Y-.  - 

?) 

1 

0.6 

0.6 

0.6 

0.6 

2 

-3.4 

-3.4 

-3.4 

-3.4 

3 

-19.4 

-19.4 

-19.4 

-19.4 

4 

22.2 

22.2 

22.2 

22.2 

Values 

of  (Y. . 
il 

-  V 

1 

23.4 

1.4 

7.4 

-23.6 

2 

-7.6 

-12.6 

42.4 

-0.6 

3 

-9.6 

1.4 

-2.6 

6.4 

4 

-10.2 

-17.2 

6.8 

-2.2 

5 


Y. 

l 

Y. 

l 

s  . 

2i 

35 

218 

43.6 

308.8 

18 

198 

39.6 

620.3 

28 

118 

23.6 

40.3 

88 

326 

65.2 

242.7 

Y  = 

860  y  = 

43.0 

-8 

-25 

-15 

45 


0.6 

-3.4 

-19.4 

22.2 


-8.6 

-21.6 

4.4 

22.8 


? 

S  =  487.053  is  the  variance  among  the  20  values  of  (Y..  -  Y) 

2  =  1J  = 

=  293.707  is  the  variance  among  the  4  values  of  (Y^-  Y) 

2 

=  303.025  is  the  average  of  the  variances  of  (Y—  -  Y^)  within 
psu’s.  Within  the  first  psu  the  variance  is: 


23. 4* 2 * 4  +  1 . 4  2  +  7 . 4  2  +  ( -  2  3 . 6)  2  +  (-8.6)2 

4 


308.8. 
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Table  4.3 

Summary  Data 

for  Six 

Apple  Trees 

1/ 

No.  of 

No .  of 

Within 

Trunk 

Total  No. 

Within 

Terminal 

Apples  on 

Tree 

Size  in 

of  Apples 

Tree 

Branches 

Terminal 

Variance 

Sq .  In. 

on  Tree 

Variance 

Branches 

DS-EP 

RP-PPS 

Tree 

N. 

l 

.  Yi 

s2.. 

2i 

X. 

i 

Yi 

S2  . 
n 

1 

13 

213 

259 

7.0 

214 

22,000 

2 

27 

1,388 

1,147 

20.0 

1,448 

478,000 

3 

26 

1,850 

2,184 

23.0 

1,901 

800,000 

4 

20 

1,592 

3,106 

16.5 

1,658 

350,000 

5 

19 

402 

241 

13 . 5 

403 

79,000 

6 

30 

1,528 

892 

19.5 

1,575 

513,000 

Total 

135 

6,973 

99 . 5 

7,199 

1/ 

The 

values  of 

N . ,  Y . ,  and 

S?.  are  from 

2i 

Table  2 

.6 

The  values  of  Y^  and 

? 

S-.  are  labeled 

2i 

Yh  and  S*h 

in  Table  2.6. 

"Path 

apples"  are 

not  included 

in 

and  S‘.. 

The  values 

of  Y!  and  S2. 
i  ri 

include 

the  path  apples  and  are 

taken 

from  Table 

3.3.  The 

subscript  "r" 

refers 

to 

random  path. 

DS-EP  and  RP-PPS  refer  to  the  method  of  sampling  a  tree  as  dis¬ 


cussed  in  Chapter  III. 


Table 

4.4 

Variances  for  Alternative  Sample  Allocations 

Components 

(1) 

m 

1 

n 

4 

V(y) 

1225.8 

Among  psu’s— ^ 

917.1 

Within  psu's 

308.7 

(2) 

2 

2 

797.6 

458.6 

339.0 

(3) 

4 

1 

583.5 

229.3 

354.2 

1/ 

Assumes 

f .  is 

negligible . 
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